Ortho Color balancing for large\u00a0UAV\u00a0ortho mosaics can be a challenging task. Usually you start to get into troubles when illumination conditions are quickly changing – at sub-optimal conditions with clouds \u00a0and low illumination condition in spring or autumn. \u00a0color balancing in Agisoft does a good job but sometime we need more – here a hot spot correction and global adaption is sometimes needed. PCIs Orthoengine does a great job here.<\/p>\n
This is a small python script that uses the ortho corrected geotiff files from Agisoft and applies a local and global correction.<\/p>\n
<\/p>\n
1. Head section of a Python PCI script<\/strong><\/p>\n Here we \u00a0load some PCI libraries.<\/p>\n 2. The xml\u00a0section of the script that defines the processing parameters<\/strong><\/p>\n where xml files are defined\u00a0– you usually want that from the user and not hard coded into your script – so f.e. using the following expression the user would have to start the script by typing:<\/p>\n python pythonscriptname \/dir\/where\/the\/files\/are<\/p>\n to start the script and to include the “InFolder” string\u00a0(in this example). We later use “InFolder” to place the string into pathname.<\/p>\n <\/p>\n","protected":false},"excerpt":{"rendered":" Ortho Color balancing for large\u00a0UAV\u00a0ortho mosaics can be a challenging task. Usually you start to get into troubles when illumination conditions are quickly changing – at sub-optimal conditions with clouds \u00a0and low illumination condition in spring or autumn. \u00a0color balancing in Agisoft does a good job but sometime we need more – here a hot spot correction and global adaption...<\/p>\n","protected":false},"author":1,"featured_media":0,"parent":18,"menu_order":8,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_uag_custom_page_level_css":"","footnotes":""},"class_list":["post-847","page","type-page","status-publish","hentry"],"uagb_featured_image_src":{"full":false,"thumbnail":false,"medium":false,"medium_large":false,"large":false,"1536x1536":false,"2048x2048":false,"post-thumbnail":false,"cd-small":false,"cd-medium":false,"cd-standard":false},"uagb_author_info":{"display_name":"S\u00f6ren Hese","author_link":"https:\/\/jenacopterlabs.de\/?author=1"},"uagb_comment_info":0,"uagb_excerpt":"Ortho Color balancing for large\u00a0UAV\u00a0ortho mosaics can be a challenging task. Usually you start to get into troubles when illumination conditions are quickly changing – at sub-optimal conditions with clouds \u00a0and low illumination condition in spring or autumn. \u00a0color balancing in Agisoft does a good job but sometime we need more – here a hot…","_links":{"self":[{"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/pages\/847","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=847"}],"version-history":[{"count":5,"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/pages\/847\/revisions"}],"predecessor-version":[{"id":852,"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/pages\/847\/revisions\/852"}],"up":[{"embeddable":true,"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=\/wp\/v2\/pages\/18"}],"wp:attachment":[{"href":"https:\/\/jenacopterlabs.de\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=847"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}\r\nprint ""\r\nprint "***************************************************************"\r\nprint ""\r\nprint " -= SHESE 4 hainich ortho processing =-"\r\nprint " workflow"\r\nprint " Last modified 24112018"\r\nprint ""\r\nprint "***************************************************************"\r\n# This script runs the following PPFs in a top-down sequence and generates a \r\n# ortho mosaic based on the input ortho images exported from AGISOFT Photoscan.\r\n# All orthos must have the same georef and resolution. \r\n# 1. mosrep\r\n# 2. mosdef\r\n# 3. mosrun\r\nimport sys \r\nimport os\r\nimport pci\r\nimport fnmatch\r\nfrom pci.fun import *\r\nfrom pci.lut import *\r\nfrom pci.fimport import fimport\r\nfrom pci.pcimod import pcimod\r\nfrom pci.model import model\r\nfrom pci.exceptions import *\r\nfrom pci.ortho import ortho\r\nfrom pci.mosprep import mosprep\r\nfrom pci.mosdef import mosdef\r\nfrom pci.mosrun import mosrun\r\nfrom pci.exceptions import PCIException\r\n<\/pre>\n
The comments are defined using a hash sign - you can also put these comments directly behind any coded line<\/pre>\n
\r\nmfile = r"W:\\work-hainich\\Hainich08-2018\\exported-orthofiles"\r\nsilfile = r"W:\\work-hainich\\Hainich08-2018\\exported-orthofiles\\mosaic.xml"\r\nnodatval = [0]\r\nnormaliz = "HotSpot" #HotSpot method often useful for airphoto imagery\r\nbalspec = "BUNDLE"\r\ncutmthd = "minsqdiff"\r\n\r\ntry:\r\n mosprep(mfile, silfile, nodatval, "", normaliz,\r\n balspec, "", "", [], cutmthd)\r\nexcept PCIException, e:\r\n print e\r\nexcept Exception, e:\r\n print e\r\n\r\n## Create Mosaic Definition XML file to\r\nsilfile = r"W:\\work-hainich\\Hainich08-2018\\exported-orthofiles\\mosaic.xml"\r\nmdfile = r"W:\\work-hainich\\Hainich08-2018\\exported-orthofiles\\mosaicdef.xml"\r\ndbic = [1,2,3]\r\ntispec = ""\r\ntipostrn = ""\r\nmapunits= "UTM 32 D000"\r\npxszout= [0.03]\r\nblend = [5]\r\nnodata = [0]\r\nftype = "PIX"\r\nfoptions = ""\r\ntry:\r\n mosdef(silfile, mdfile, dbic, tispec, tipostrn, mapunits, pxszout, blend, nodata, ftype, foptions)\r\nexcept PCIException, e:\r\n print e\r\nexcept Exception, e:\r\n print e\r\n\r\n## Finally create the full resolution mosaic\r\noutdir = r"W:\\work-hainich\\Hainich08-2018\\exported-orthofiles"\r\ncrsrcmap = ""\r\n\r\ntry:\r\n mosrun(silfile, mdfile, outdir, "", crsrcmap,\r\n "", "", "", "BILIN")\r\nexcept PCIException, e:\r\n print e\r\nexcept Exception, e:\r\n print e\r\n<\/pre>\n