Odd error when using mapproject

[Johanna Alen-Bella]

Hello everyone,

I'm running ASP version 3.3.30 in a recent Ubuntu install and getting the following error when trying to mapproject a mosaic .tif of rock glaciers in the Sierra Nevadas onto a DEM (Copernicus, converted to WGS84):

mapproject -t rpc --tr 2.0 dem-adj.tif *_WVCORRECT.r100.tif *_WVCORRECT.r100.xml 016031408010_01_MP.tif

mapproject_single --query-projection dem-adj.tif 016031408010_01_P001_WVCORRECT.r100.tif 016031408010_01_P001_WVCORRECT.r100.xml 016031408010_01_MP.tif -t rpc --tr 2.0
        --> Setting number of processing threads to: 4
Using session: rpc
Loading camera model: 016031408010_01_P001_WVCORRECT.r100.tif 016031408010_01_P001_WVCORRECT.r100.xml
Computed image to DEM bbox: (Origin: (241.381, 37.0085) width: 0.217001 height: 0.18502)
Output pixel size: 2
Projected space bounding box: (Origin: (121, 19) width: 0 height: 0)
Image box: (Origin: (2147483647, 2147483647) width: 0 height: 0)
Output image size:
(width: -2147483648 height: -2147483648)
Query finished, exiting mapproject tool.

Output image size is -2147483648 by -2147483648 pixels.
Splitting into -419430 by -419430 tiles.
parallel --will-cite --workdir /mnt/g/Programs/Stereo-Pipeline/ThesisImages/RockGlaciers/1020010023E6E300_016031408_10_0/016031408010_01_003/016031408010_01_003/016031408010_01 -u --env PATH --env PYTHONPATH --env ISISROOT --env ASP_LIBRARY_PATH --env ISISDATA -a ./016031408010_01_MP_tif_tiles/argumentList.txt -P 16 --colsep \t /mnt/g/programs/stereo-pipeline/StereoPipeline-3.3.0-Linux/bin/python /mnt/g/programs/stereo-pipeline/StereoPipeline-3.3.0-Linux/libexec/mapproject --pixelStartX {1} --pixelStartY {2} --pixelStopX {3} --pixelStopY {4} --work-dir ./016031408010_01_MP_tif_tiles/ dem-adj.tif 016031408010_01_P001_WVCORRECT.r100.tif 016031408010_01_P001_WVCORRECT.r100.xml 016031408010_01_MP.tif --threads 8 -t rpc --tr 2.0
Removing: ./016031408010_01_MP_tif_tiles/
Traceback (most recent call last):
  File "/mnt/g/programs/stereo-pipeline/StereoPipeline-3.3.0-Linux/libexec/mapproject", line 557, in <module>
    sys.exit(main(sys.argv[1:]))
  File "/mnt/g/programs/stereo-pipeline/StereoPipeline-3.3.0-Linux/libexec/mapproject", line 504, in main
    raise Exception("No mapprojected tif tiles were generated")
Exception: No mapprojected tif tiles were generated

Any advice would be welcome!

Thank you,
Johanna Alen-Bella

[Oleg Alexandrov]

I think your DEM is in degrees, and you are trying to mapproject in meters. You should set a local projection with --t_srs. I will look later into having some better error message.

--
You received this message because you are subscribed to the Google Groups "Ames Stereo Pipeline Support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ames-stereo-pipeline...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ames-stereo-pipeline-support/302c2044-c256-4449-8d4f-9a0c347cb174n%40googlegroups.com.