Hello - I'm using two tools, the "Create Orthocorrected Raster" and also the "Geometric" raster function to apply the RPC file for a WorldView 3 image and a high resolution DEM to orthorectify my image. The raster functions work great, I get an offset of ~1-3 pixels but they don't use nearest neighbor (NN) resampling which is imperative for my work flow. The "Create Orthocorrected Raster" tool lets me choose NN resampling, but the image is 200m-1200m off. Has anyone else figured this one out? Thank you!
can you give us more information? Are you using ArcGIS Pro 2.3? Can you describe how you are applying the raster function?
You should be able to:
Does this not work?
Alternatively if you add your WV image into a mosaic dataset, you can apply the function to the mosaic dataset.
I'm testing https://pro.arcgis.com/en/pro-app/tool-reference/data-management/create-ortho-corrected-raster-datas... now but I have a pre-release build of Pro 2.4. If you're getting a large horizontal error, I'd have to confirm if we had a bug in a prior (or current) version.
I'm sorry I didn't give you the specifics. I'm using vs 10.6.1, ArcCatalog and ArcMap depending on which tool I'm using. When using the raster functions:
1. I'm "adding" a NTF file (using the .IMD) to ArcMap.
2. I then use the geometric function with my DEM (10m USGS DEM) and the RPC file is used.
3. I tried using the resampling function set to NN and it didn't work, but maybe I had it in the wrong place in the function chain?
I've tried the ArcGIS Desktop vs 10.6.1 Create Orthocorrected Raster tool, which I would prefer to use but for some reason it doesn't geocorrect the image, it warps it a bit, but registration is seriously off, like 200m-1200m.
Hi Cody, It seems that I can add the NN resampling function to the raster function chain and it's working. It looks like I can use the "Generate Raster from Raster function" tool in ArcMap and ArcCatalog to create an output raster in geoTIF format. I need to be able to code this in arcpy, so it seems like this might all work. I'm still trying to figure out how to use the raster function template in my arcpy code. I'll post again once I'm sure this all works. Thank you! Leanne
I'm curious how many satellite scenes you have, and how you'll use them after orthorectification. I generally recommend against processing an image to create a new output file, unless you have no alternative. Doing this will double your disk space, and if you find anything that you need to change, you need to take more time and disk space to do it over again. The alternative - if you're doing all of your work in ArcGIS - is to use the mosaic dataset and generate the required products when accessed, as virtual datasets.
If you're not using the mosaic dataset, I would encourage you to explore if it will meet your requirements more easily/automatically - e.g. it can load multiple raw scenes into a single data structure, and orthorectify them on-the-fly. You can view and work with any individual scene, or work with a mosaic of multiple scenes. If you need more information on this, let me know and I'll show you where to find it in the ArcGIS Imagery Workflows as well as our online Help.
I just tried to use the "Create ortho-corrected dataset" tool on a mosaic that I created from multiple worldview-2 rasters, and I received the error message "This tool does not support mosaic datasets". I am using ArcGIS 2.4. If this is correct, not sure about the proper workflow for creating mosaic of worldview images.
the help for the tool you mention states that it is meant to be used with single raster datasets only. The good news is, as Cody stated above, you can apply a function in any MosaicDataset that does the DEM integration on the fly. To do so - using your worldview-2 datasets, proceed like this:
"Behind the scenes, every item will now have a DEM properly added to the items "Geometry" function (See other example used for screenshot 3),
which you can see if you open the Raster Item explorer and the function chain for any item
Hope this answers it
Thanks for the detailed instructions. Unfortunately, when I followed this workflow, the resulting dataset was not orthorectified and only 4 bands were included. Eight bands should have been available. Not sure how to proceed. It seems like the best option may be to process the individual rasters in a loop and then add to a mosaic. I know this isn't ideal, but seems to be the most stable.
there are different "flavours" of WorldView. For general introductions, please check this: Managing High-Resolution Satellite Imagery—Imagery Workflows | ArcGIS - which also includes a link to access additional python scripts that can help you with some WoldView "standards"
To check yourselves, you have several options:
1) Remember WorldView is a "Raster Type" - so you can go to Catalog, scroll down the appropriate folder to the WorlView Metadata-File
When you see the Sat-Icon, just dropdown and see what this data supports
2) If you know what you have and want - go back to "Raster Type Properties" (Screenshot 2 above previous answer), this time to the "general" Tab:
There is a list of ProductTypesw, depending on the WorldView version you have. Choose from this
Hope this helps