def zoom_to_tiles(self):
if self.layer_to_edit.navigation.is_valid:
# compute the wrapper extent of all tiles
extent_tiles = QgsRectangle()
[
extent_tiles.combineExtentWith(tile.extent)
for tile in self.layer_to_edit.navigation.tiles
]
extent_tiles = extent_tiles.buffered(
(extent_tiles.yMaximum() - extent_tiles.yMinimum()) * 0.02)
self.render_widget.update_canvas_to(extent_tiles)
self.render_widget.canvas.refreshAllLayers()
else:
self.render_widget.update_canvas_to(self.layer_to_edit.extent())
def UpdateLayers(packet, parent=None, mosaic=False, group=None):
''' Update Layers Values '''
gv.setGroupName(group)
groupName = group
# frameCenterLat = packet.FrameCenterLatitude
# frameCenterLon = packet.FrameCenterLongitude
gv.setFrameCenterElevation(packet.FrameCenterElevation)
gv.setSensorLatitude(packet.SensorLatitude)
gv.setSensorLongitude(packet.SensorLongitude)
sensorTrueAltitude = packet.SensorTrueAltitude
gv.setSensorTrueAltitude(sensorTrueAltitude)
sensorRelativeElevationAngle = packet.SensorRelativeElevationAngle
slantRange = packet.SlantRange
OffsetLat1 = packet.OffsetCornerLatitudePoint1
LatitudePoint1Full = packet.CornerLatitudePoint1Full
UpdatePlatformData(packet, hasElevationModel())
UpdateTrajectoryData(packet, hasElevationModel())
frameCenterPoint = [
packet.FrameCenterLatitude, packet.FrameCenterLongitude,
packet.FrameCenterElevation
]
# If no framcenter (f.i. horizontal target) don't comptute footprint,
# beams and frame center
if (frameCenterPoint[0] is None and frameCenterPoint[1] is None):
gv.setTransform(None)
return True
# No framecenter altitude
if (frameCenterPoint[2] is None):
if (sensorRelativeElevationAngle is not None
and slantRange is not None):
frameCenterPoint[2] = sensorTrueAltitude - \
sin(sensorRelativeElevationAngle) * slantRange
else:
frameCenterPoint[2] = 0.0
# qgsu.showUserAndLogMessage("", "FC Alt:"+str(frameCenterPoint[2]), onlyLog=True)
if OffsetLat1 is not None and LatitudePoint1Full is None:
if hasElevationModel():
frameCenterPoint = GetLine3DIntersectionWithDEM(
GetSensor(), frameCenterPoint)
CornerEstimationWithOffsets(packet)
if mosaic:
georeferencingVideo(parent)
elif OffsetLat1 is None and LatitudePoint1Full is None:
if hasElevationModel():
frameCenterPoint = GetLine3DIntersectionWithDEM(
GetSensor(), frameCenterPoint)
CornerEstimationWithoutOffsets(packet)
if mosaic:
georeferencingVideo(parent)
else:
cornerPointUL = [
packet.CornerLatitudePoint1Full, packet.CornerLongitudePoint1Full
]
if None in cornerPointUL:
return False
cornerPointUR = [
packet.CornerLatitudePoint2Full, packet.CornerLongitudePoint2Full
]
if None in cornerPointUR:
return False
cornerPointLR = [
packet.CornerLatitudePoint3Full, packet.CornerLongitudePoint3Full
]
if None in cornerPointLR:
return False
cornerPointLL = [
packet.CornerLatitudePoint4Full, packet.CornerLongitudePoint4Full
]
if None in cornerPointLL:
return False
UpdateFootPrintData(packet, cornerPointUL, cornerPointUR,
cornerPointLR, cornerPointLL, hasElevationModel())
UpdateBeamsData(packet, cornerPointUL, cornerPointUR, cornerPointLR,
cornerPointLL, hasElevationModel())
SetGCPsToGeoTransform(cornerPointUL, cornerPointUR, cornerPointLR,
cornerPointLL, frameCenterPoint[1],
frameCenterPoint[0], hasElevationModel())
if mosaic:
georeferencingVideo(parent)
UpdateFrameCenterData(frameCenterPoint, hasElevationModel())
UpdateFrameAxisData(packet.ImageSourceSensor, GetSensor(),
frameCenterPoint, hasElevationModel())
# detect if we need a recenter or not. If Footprint and Platform fits in
# 80% of the map, do not trigger recenter.
f_lyr = qgsu.selectLayerByName(Footprint_lyr, groupName)
p_lyr = qgsu.selectLayerByName(Platform_lyr, groupName)
t_lyr = qgsu.selectLayerByName(FrameCenter_lyr, groupName)
iface = gv.getIface()
centerMode = gv.getCenterMode()
if f_lyr is not None and p_lyr is not None and t_lyr is not None:
f_lyr_out_extent = f_lyr.extent()
p_lyr_out_extent = p_lyr.extent()
t_lyr_out_extent = t_lyr.extent()
# Default EPSG is 4326, f_lyr.crs().authid ()
# Disable transform if we have the same projection wit layers anf
# canvas
epsg4326 = "EPSG:4326"
curAuthId = iface.mapCanvas().mapSettings().destinationCrs().authid()
if (curAuthId != epsg4326):
xform = QgsCoordinateTransform(
QgsCoordinateReferenceSystem(epsg4326),
QgsCoordinateReferenceSystem(curAuthId),
QgsProject().instance())
transP = xform.transform(
QgsPointXY(
list(p_lyr.getFeatures())[0].geometry().asPoint().x(),
list(p_lyr.getFeatures())[0].geometry().asPoint().y()))
transT = xform.transform(
QgsPointXY(
list(t_lyr.getFeatures())[0].geometry().asPoint().x(),
list(t_lyr.getFeatures())[0].geometry().asPoint().y()))
rect = list(f_lyr.getFeatures())[0].geometry().boundingBox()
rectLL = xform.transform(
QgsPointXY(rect.xMinimum(), rect.yMinimum()))
rectUR = xform.transform(
QgsPointXY(rect.xMaximum(), rect.yMaximum()))
f_lyr_out_extent = QgsRectangle(rectLL, rectUR)
t_lyr_out_extent = QgsRectangle(transT.x(), transT.y(), transT.x(),
transT.y())
p_lyr_out_extent = QgsRectangle(transP.x(), transP.y(), transP.x(),
transP.y())
bValue = iface.mapCanvas().extent().xMaximum() - iface.mapCanvas(
).center().x()
# create a detection buffer
map_detec_buffer = iface.mapCanvas().extent().buffered(bValue * -0.7)
# recenter map on platform
if not map_detec_buffer.contains(p_lyr_out_extent) and centerMode == 1:
# recenter map on platform
iface.mapCanvas().setExtent(p_lyr_out_extent)
# recenter map on footprint
elif not map_detec_buffer.contains(
f_lyr_out_extent) and centerMode == 2:
# zoom a bit wider than the footprint itself
iface.mapCanvas().setExtent(
f_lyr_out_extent.buffered(f_lyr_out_extent.width() * 0.5))
# recenter map on target
elif not map_detec_buffer.contains(
t_lyr_out_extent) and centerMode == 3:
iface.mapCanvas().setExtent(t_lyr_out_extent)
# Refresh Canvas
iface.mapCanvas().refresh()
return True
