Changes between Version 25 and Version 26 of Tutorials/Digital Twins/Geo2Sig Map Getting Started

Timestamp:

Jan 27, 2026, 4:10:44 AM (2 weeks ago)

Author:

srivera2

Comment:

—

Tutorials/Digital Twins/Geo2Sig Map Getting Started

@@ -89 +89 @@
    {{{
    #!python
-   # Core imports
-   import numpy as np
-   import mitsuba as mi
-   import matplotlib as mpl
-   import matplotlib.pyplot as plt
-   import sionna
-   from sionna.rt import load_scene, Camera, Transmitter, Receiver, PlanarArray,\
-                         PathSolver, RadioMapSolver, load_mesh, watt_to_dbm, transform_mesh,\
-                         cpx_abs_square  
+# Core imports
+import numpy as np
+import mitsuba as mi
+from sionna.rt import load_scene, Transmitter, Receiver, RadioMapSolver, transform_mesh
    }}}
 
@@ -104 +99 @@
    {{{
    #!python
-   SCENE_FOLDER = "../scenes/COSMOS"
-   scene = load_scene(f"{SCENE_FOLDER}/scene.xml")
+import os
+current_dir = os.getcwd()
+print(current_dir)
+SCENE_FOLDER = "../scene/COSMOS"
+scene = load_scene(f"{SCENE_FOLDER}/scene.xml")
 
-   # Interactive 3D visualization and view of the scene
-   scene.preview();   
+# Interactive 3D visualization and view of the scene
+scene.preview();    
    }}}
 
@@ -115 +113 @@
    {{{
    #!python
-   # Set the operating frequency (n48 band for 5G)
-   scene.frequency = 3.7e9  # 3.7 GHz
+   from sionna.rt import AntennaArray
+from sionna.rt.antenna_pattern import antenna_pattern_registry
 
-   # Define Rx position
-   ue_position = [10.0, 0.0, 0.0]   # Rx position (x, y, z in meters)
+# Set the operating frequency (n48 band for 5G)
+scene.frequency = 3.7e9  # 3.7 GHz
 
-   # Define Tx position
-   gnb_position = [0.0, 0.0, 20.0]
+# Define Rx position
+ue_position = [20.0, 20.0, 0.0]   # Rx position (x, y, z in meters)
 
-   # gNB antenna: 3GPP TR 38.901 pattern (AIRSTRAN D 2200)
-   gnb_pattern_factory = antenna_pattern_registry.get("tr38901")
-   gnb_pattern = gnb_pattern_factory(polarization="V")
-   
-   # Use isotropic antenna pattern
-   ue_pattern_factory = antenna_pattern_registry.get("iso")
-   # Polarization should match the transmitter
-   ue_pattern = ue_pattern_factory(polarization="V")
+# Define Tx position
+gnb_position = [0.0, 0.0, 70.0]
 
-   # SISO: Single antenna element at origin [0, 0, 0] for both TX and RX
-   single_element = np.array([[0.0, 0.0, 0.0]])  # Shape: (1, 3)
+# gNB antenna: 3GPP TR 38.901 pattern (AIRSTRAN D 2200)
+gnb_pattern_factory = antenna_pattern_registry.get("tr38901")
+gnb_pattern = gnb_pattern_factory(polarization="V")
 
-   # Configure antenna arrays
-   scene.tx_array = AntennaArray(
-       antenna_pattern=gnb_pattern,
-       normalized_positions=single_element.T  # Shape: (3, 1)
-   )
+# Use isotropic antenna pattern
+ue_pattern_factory = antenna_pattern_registry.get("iso")
+# Polarization should match the transmitter
+ue_pattern = ue_pattern_factory(polarization="V")
 
-   scene.rx_array = AntennaArray(
-       antenna_pattern=ue_pattern,
-       normalized_positions=single_element.T  # Shape: (3, 1)
-   )  
+# SISO: Single antenna element at origin [0, 0, 0] for both TX and RX
+single_element = np.array([[0.0, 0.0, 0.0]])  # Shape: (1, 3)
 
-   # Create Rx
-   rx = Receiver(name="ue", position=ue_position, display_radius=0.03)
-   scene.add(rx)
+# Configure antenna arrays
+scene.tx_array = AntennaArray(
+    antenna_pattern=gnb_pattern,
+    normalized_positions=single_element.T  # Shape: (3, 1)
+)
 
-   # Create Tx
-   tx = Receiver(name="gnb", position=gnb_position, display_radius=0.03)
-   scene.add(tx)
+scene.rx_array = AntennaArray(
+    antenna_pattern=ue_pattern,
+    normalized_positions=single_element.T  # Shape: (3, 1)
+)  
 
+# Create Rx
+rx = Receiver(name="ue", position=ue_position, display_radius=0.03)
+scene.add(rx)
+
+# Create Tx
+tx = Transmitter(name="gnb", position=gnb_position, display_radius=0.03)
+scene.add(tx)
    }}}
 
@@ -162 +162 @@
    #!python
 rm_solver = RadioMapSolver()
-
 
 # Create a measurement surface by cloning the terrain
@@ -174 +173 @@
     los=True,
     refraction=False,
-    cell_size=(2., 2.),
-    size=[512, 512],          # Size of the radio map in meters. The number of cells in the x and y directions should be [size/cell_size]
+    cell_size=(1.0, 1.0),
+    size=[1000, 1000],          # Size of the radio map in meters. The number of cells in the x and y directions should be [size/cell_size]
     # orientation=[0, 0, 0],    # Orientation of the radio map plane
     measurement_surface=measurement_surface, # TO BE UDPATED