Forum

Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - Midkhat

Pages: [1]
1
Python and Java API / Re: Compound CS transformation
« on: March 29, 2022, 02:14:24 AM »
@Paulo
Thank you for your reply, I have sorted that out then and got busy with some other stuff and missed your reply, but I have recently encountered similar issue and have found your reply to my own post through search :)
The issue I have now - I get the Missing Vertical Datum error when import shape as Outerboundary (from .kml) and try to export raster or point cloud:
1. Export Raster (DEM or Ortho). I get error when Z component is present in .kml. As the kml is created in third party software - it is possible that it is indeed not right value. I don't really need Z in this shape, I just remove Z and this sorts it - DEM and Ortho come out clipped.
2. Export Points. In this case removing Z doesn't help as I presume it is required for point cloud clipping. At the same time I can open the project later and export clipped cloud from GUI...

Would you be able to advise or may be there is a work around for point cloud clipping?

PS Thanks again, I really appreciate your help and have learnt some new stuff from the prnscreens you shared (and other your comments in this forum)

2
Bug Reports / failed: error 1286
« on: January 24, 2022, 03:31:57 AM »
I get bunch of following two errors since switched to 1.8:

2022-01-24 11:13:01 OpenGL Vendor: Apple
2022-01-24 11:13:01 OpenGL Renderer: Apple M1
2022-01-24 11:13:01 OpenGL Version: 2.1 Metal - 76.3
2022-01-24 11:13:01 Maximum Texture Size: 16384
2022-01-24 11:13:01 Quad Buffered Stereo: not enabled
2022-01-24 11:13:01 ARB_vertex_buffer_object: supported
2022-01-24 11:13:01 ARB_texture_non_power_of_two: supported
2022-01-24 11:13:01 Agisoft Metashape Professional Version: 1.8.0 build 13794 (64 bit)
2022-01-24 11:13:01 Platform: Mac OS
2022-01-24 11:13:01 RAM: 16.0 GB
2022-01-24 11:13:02 LoadProject: path =
2022-01-24 11:13:02 loaded project in 0.008163 sec
2022-01-24 11:13:02 Finished processing in 0.008331 sec (exit code 1)
2022-01-24 11:13:02 renderer_qgl.cpp line 659: blitFramebuffer(0, 0, width, height, 0, 0, width, height, 0x00000100, 0x2600) failed: error 1286
2022-01-24 11:13:02 renderer_qgl.cpp line 674: glBlitFramebuffer(0, 0, width, height, 0, 0, width, height, 0x00000100, 0x2600) failed: error 1286
....
2022-01-24 11:13:05 renderer_qgl.cpp line 659: blitFramebuffer(0, 0, width, height, 0, 0, width, height, 0x00000100, 0x2600) failed: error 1286
2022-01-24 11:13:05 renderer_qgl.cpp line 674: glBlitFramebuffer(0, 0, width, height, 0, 0, width, height, 0x00000100, 0x2600) failed: error 1286

They only appear when I create a shape, or import a kml (in both cases - from GUI or API), save and reopen the project.
The errors keep streaming in even when I only move a mouse (or Pan, Zoom, click) in Model view. They stop coming in when switch to Ortho view.

3
Python and Java API / Compound CS transformation
« on: December 15, 2021, 02:00:40 AM »
Hi

I am involved in a drone data processing. The photos may come in GDA94 or GDA2020 and the output has to be in MGA94 with AHD geoid height.
I have figured out how to do that in Agisoft UI with the help of the manual and it works fine, but I am having trouble automating that process in API.
This is the WKT I am using when transforming GDA2020 to MGA94:
Code: [Select]
COMPD_CS["GDA94 / MGA zone 55 + AHD",
    PROJCS["GDA94 / MGA zone 55 + AHD",
        GEOGCS["GDA94",
            DATUM["Geocentric Datum of Australia 1994",
                SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],
                TOWGS84[0.06155,-0.01087,-0.04019,0.0394924,0.0327221,0.0328979,-0.009993999999999999],AUTHORITY["EPSG","9688"]],
                PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
                UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9102"]],
            AUTHORITY["EPSG","4283"]
            ],
            PROJECTION["Transverse_Mercator",AUTHORITY["EPSG","9807"]],
            PARAMETER["latitude_of_origin",0],
            PARAMETER["central_meridian",147],
            PARAMETER["scale_factor",0.9996],
            PARAMETER["false_easting",500000],
            PARAMETER["false_northing",10000000],
            UNIT["metre",1,AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","28355"]
    ],
    VERT_CS["AHD height",
        VERT_DATUM["Australian Height Datum",2005,AUTHORITY["EPSG","5111"]],
        UNIT["metre",1,AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","5711"]
    ]
]

When I try to apply this transformation to all cameras I get an error:
RuntimeError: Vertical datum missing

I have added the geoid tiff to the "C:/Program Files/Agisoft/Metashape Pro/geoids" folder, trying to add it to the Metashape.CoordinateSystem as well. Have also tried to add it to the script folder.

This is related piece of code:
Code: [Select]
   
wgs = Metashape.CoordinateSystem("EPSG::4326")

geoid_file = "C:/Program Files/Agisoft/Metashape Pro/geoids/au_ga_AUSGeoid09_V1.01.tif"
if not os.path.exists(geoid_file):           
    print(f"file {geoid_file} doesn't exist!!!!!")
mga_file = "./mine-sites/MGA55_AHD.prj"   
with open(mga_file) as f:
    wkt = f.read()
mga = Metashape.CoordinateSystem(wkt)
mga.addGeoid(geoid_file)   

for camera in cameras:
    camera.reference.location = Metashape.CoordinateSystem.transform(camera.reference.location, wgs, mga)
chunk.crs = mga
doc.save()

Best regards

Pages: [1]