Home

OSSIM Users Guide - OSGeo Server

1. z0 Undo Crrl Zz Propagate Refresh Redo Ctri Y View Geometry adjustment Image chain Layers Vector Position Information Ej mui Image Info Figure 31 Access Imagelnfo menu option Page 37 of 107 Go to the edit gt image info menu option and select the Image Projection tab and select the ossimBilinearProjection class Non commercial Image Information 2 x Image Layer layer 1 image chein C date dra bettlemn tif ER Image Geometry Image Projection Valid Vertices Projection ossimEquDistOylProjection z Unknown m Deum Sensor Model Units ossimBilinearProjection ossimAlbersProjection Pixel type ossimAzimEquDistProjection assimBonneProjection Zone ossimBngProjection Hemisphere ossimCassiniProjection ___JossimCylEquAreaProjection Origin Latitude gs simEquDistOyiProjection z Central Meridian 000 00 0000000000 E False Easting 0 000000000000000 False Northing P eon000000000000 Tie Point 000 00 00 00000000 E Tie Point Y 00 00 D0 00000000 N Scale X 0 000008883152841 Scale Y 0 000008983152841 Scale Factor STD Parallel 1 00 00 00 00000000 N STD Parallel 2 00 O0 00 00000000 N Save Reset Help Close d Figure 32 Change to a bilinear projection Once selected you will be prompted to edit the projection by using the drop point GUI This projecti
2. Choose a fies to open ixl ban Demm v EY EHE L arian cof nercial Vce Components F 4 iners Fiters Spores Displeys S ler My Documents Filename z pen Flescttyos mages act ect dem DEN ctl di at2 hoi Cancel ES Close Figure 82 Drag the ImageHandler icon onto the canvas and the file browser opens automatically Select the image under test_data sanfran sanfran ccf Page 66 of 107 Z Non Commercial Vce Components Combiners Fiters Exporters Displays ImegeHendler Figure 83 sanfran ccf file loaded manually in VCE Now change the components dialog to the Filters tab Ieee je Fiere Emons Dieve PICorvolutonF ACorrGridRema AlConRemapoe BendavarageFil iter poer r et Figure 84 Filters tab for the VCE components dialog Select the ImageRenderer component with the left mouse button While holding down the left mouse button drag the ImageRenderer icon onto the canvas Now connect the top red triangle on the ImageHandler icon to the left red triangle on the ImageRenderer filter Page 67 of 107 Bk Ed X REL S IP uE Figure 85 Connected renderer Now select the displays tab and drag a display onto the canvas MEE B i Ele Layer Vor tities Window Hep Visual Chain Editor i Sk
3. Sample ossim_preferences file Id ossim preferences template v 1 11 2003 03 27 20 51 20 gpotts Exp Description ossim preferences template This file will be automatically loaded by ossim applications provided the environment variable OSSIM PREFS FILE is set to point to some form of this file Note c comments can be used to comment a line To set the environment variable for automatic preference file loading This assumes a preference file in your home called ossim preferences Typically this would be put in a dot file read at startup of a shell tcsh and csh users setenv OSSIM PREFS FILE ossim preferences bash and sh users export OSSIM PREFS FILE ossim preferences windoze users I ll have to look this up You can also use the P preference file option at application startup where preference file is full path and filename Keyword dted directory Path to DTED elevation cells Can have more than one diretory dted directoryl Data elel DTED3arc dted directory2 Data elel DTEDlkps if Keyword dted_cell Path to a DTED elevation cell Keyword default_elevation_path Default path for the elevation manager popup Add to start at default_elevation_path Data elel Keyword elevation enabled Page 20 of 107 If disabled calls to the elevation manager getHeightAboveMSL and getHeightAboveEllipsoid will return a null height
4. m Open Project Close Project hue ceai o e A A O A e A EENE DaVeProJ cb eeri enn DESCR UU HT ERE Save Project As Display Menuss PERIN LO Zoom and Pan Modes Geometry Adjustment Image Chain Image Info Drop Point GUL c E ii 36 Image Generation isisisi nea a a a a i 40 SWipe MMU C MI Add Layer Remove B D e enne e EUR ERR E REPE ER RO E ERE RENE Eb EE MEE rte ke dee 46 osten teed 46 Page 2 of 107 boerg e aaa eT 46 Enhancements Menu eee eee sees eee eene teneo eoe en senses esses sns ens ens ense essen sens enseosess ees seseessessesses FT Band Selection 4 a te e C ED ce ea eo Port ere eee te vei vel o vu rode e e NE 47 Brightness Contfast 2 eec eren ODD dI Ee E E EO 48 Sosczei o dkeeparidnrec E 48 Histogram Operations eie peer HERR RIETI ERE EEPE E ASEESSA 49 Hue Saturation Intensity Adjustments eese nennen enne enne nennen nnne 50 Propagate Resampler tee eere e eR ee ien ek peres eae aE R teres eT 50 Iced la ee i 3 VGE Components eene ee yep e De doe e EE EU Ie cst M RR NENNEN Le EORR RUE EN ERN Edna 64 Ch in Editor WINGO Wieser ee ero none M a E t nete S PM PR PER ER RE E IRE EE TR Tenere 65 Utilities Fori
5. start sample 100 Will start at sample 100 p or end sample ou can specify the offset in the sample direction on where to stop the extraction By default this is max sample Sample use end sample 100 Will stop at sample 100 create overview Specifies to create an overview for the output ile being generated q or compression quality Specifies the compression quality This is currently only valid for jpeg output t r or res level Specifies the resolution level to extract from he input source By default this is RO or full resolution Resolutions are numbered from 0 o Max 1 t or create thumbnail akes an argument specifying the dimension of he thumbnail Note This just grabs the closest r level that fits that dimension Sample Use create thumbnail 512 Finds the closest r level with closest dimensions o 512 W or tile width Defines the tile width for the handlers that Support tiled output For example iff_tiled_band_separate Page 83 of 107 Examples 1 icp jpeg foo tif myfile jpg Will copy foo tif to myfile jpg where the output type is specified as jpeg 2 icp bands 1 2 tiff tiled foo tif output tif Will copy bands 1 and 2 from foo tif and put them into output tif At the time of this document the standard OSSIM formats supported are tiff strip tiff strip band separate tiff tiled tiff tiled band separate jpeg general raster bip general raster bil general raster bsq genera
6. default true Use true yes y or 1 to enable false no n or 0 to disable elevation enabled true rm Keyword elevation auto_load_dted enabled Enable autoloading of dted elevation as need from above directories default true Use true yes y or 1 to enable false no n or 0 to disable elevation auto load dted enabled true eee Keyword elevation auto_sort enabled Enable sorting of elevation cells on an add by the lowest best mean post spacing default true Use true yes y or 1 to enable fi false no n or 0 to disable elevation auto sort enabled true Please Geoid support GEOID 99 Set keyword to the directory containing the GEOID 99 grids geoid 99 directory data elel geoid geoid99 GEOID EGM 96 Set keyword to the path to the egm96 grd geoid egm 96 grid data elel geoid geoid96 egm96 grd font support font filel usr share fonts default Typel b0180321 pfb font dirl usr share fonts default Typel for MacOSX font filel System Library Fonts Helvetica dfont font dirl System Library Fonts plugin support plugin dirl work gpotts ossim plugin You can also list files individually plugin filel Users Shared Development imglnk lib libossim plugins dylib Page 21 of 107 jhe Keywords for plugin registration plugin dir
7. Combine Feather Again just select the first two lines the source files and press the Apply button Clean up the displays by selecting Window gt Tile Select any of the displays press the Fit button then press the Propagate button and you should end up with something similar to this E Non Commercial unnamedo prj loj x File Layer Vce Utilities Window Help joj x joj x File Edit Swipe Enhancements Help File Edit Swipe Enhancements Help Zoom Pan Zoom Pan Fit Fullres Propagate Refresh nea Fit Fulres Propagate Refresh nea Ele Edit Swipe Enhancements Help Zoom Pan Fit Fullres Propagate Refresh nea 37 49 33 9830 N 122 22 19 2024 W 32 1551 m rgb 1 Figure 13 Three basic mosaics The difference between these three combiners is how it treats pixels in the overlap area A mosaic simply chooses one source over the others the blend averages the values of the pixels and the feather changes the blend gradually as the distance from the seam increases To change the weighting sources or any relevant adjustable parameters select the desired display and Edit gt Layers to bring up the appropriate dialog box Page 15 of 107 3 4 Save your project Main Window File gt Save Project as mapblend prj 3 5 Summary This completes the quickstart tutorial for some of the combiner algorithms available in
8. 0 0 Figure 3 Cascade Select Window gt Minimize All Collapses the displays into the bottom of the Canvas Image chains 35 x Image chain me n x Figure 4 Minimize all Page 9 of 107 Select Window gt Restore All Will restore the minimized images to their original location To place the images side by side use the Tile menu command Select Window gt Tile 8 Non Commercial unnamedo prj loj xl File Layer Vce Utilities window Help Page char BIETET oce an cdo test data fenem TES File Edit Swipe Enhancements Help File Edit Swipe Enhancements Help Zoom Per zoom Panj Fullres Propagate Refresh nea Fulres Propagate Refresh nee 37956 27 0263 N 122 15 25 1861 W 32 016 m rgb 12 Z 37 54 52 3305 N 122926 36 3857 W 32 156 m rgb 0 Figure 5 Tiled Page 10 of 107 3 3 Pan Zoom and Propagate in the Displays The displays have several modes and capabilities that make it easy to navigate and synchronize displays within a project In the center of the Landsat scene is Treasure Island place the cursor over the island and click Note that the Latitude and Longitude of the mouse click is displayed at the bottom of the window For more thorough position information open the Edit gt Position information menu option found on the display window Click on Treasure Island Click on the Full res button You should s
9. Eat Loaders Combiners Fits Exporers Ciepiays EXE HERRERA FP INon corrmerdall vas cariponen HE Figure 86 Final reprojection chain The view on the ossimlmageRenderer can be modified through the Display window Edit gt View See 4 2 5 7 3 2 2 Building Chains from the Layer Manager The VCE canvas supports drag and drop of selected chains from the Layer gt Manager dialog onto the VCE canvas So we have the same data loaded as the tutorial clear the data manager of all information or just re start the ImageLinker application Once the data manager is cleared load images sanfran sanfran ccf and sanfran map highres utm map tif under the test data sub directory Once they are loaded close both displays The data is still in the Page 68 of 107 global Data Manager even though the displays have been closed Next open the Layer gt Manager and select both images and drag with the mouse cursor the chains onto the VCE canvas editor Visual Chain Editor Blxl Ble Edt IS Ei a PI fai E EEEE E Figure 87 Chains dragged from the Layer gt Manager window Select the Displays tab on the Vce Components dialog and select and drag and drop the display onto the canvas six Ele Edit EST HEME ee E eee C 4 H Figure 88 Display added to the canvas Now hold the left mouse button down on the red triangle of the right most filte
10. The final product can be saved by hitting the save image button It will save the image and immediately after saving bring it back up into ImageLinker for further viewing processing If you are using imagelinker to create a very large mosaic and your output product will take a long time to execute then it might make since to output the spec file only by pressing the save spec file only button The file you output can be fed into the igen command line application For example if you saved out a spec file to foo spec then you can execute it on the command line by igen foo spec You must have your path set to the location of the executable or give a full qualified path to the igen executable Page 44 of 107 7 3 Swipe Menu If we have brought in multiple images also know as layers in OSSIM we can swipe the datasets together We first add the layer by accessing the Swipe gt Add layer menu option The second layer is selected through the layer selection dialog box Two geographically coincident data sets can be stacked and rapidly compared by moving the cursor back and forth over the content region of the window The top layer is covered and uncovered as the mouse is moved 7 3 1 Add Layer File Layer Vce Utilities Window Help ta T jojx Image chain C data test_data sanfran sanfran ccf Ele Edit Swipe Enhancements Help Add layer Figure 46 Adding another layer to a display for swiping The add layer menu item
11. Try to get an input projection from the image and replace it with a coarse grid model 2 ogeom2ogeom cg noelev image Will create a coarse grid model and disable elevation Note Elevation is applied but after it is applied it is no longer affected by elevation 10 7 applanix20geom The application gives a quick way to generate OSSIM geometry files given the Applanix parameters The applications can only support ECEF orientations where the platform positions are in ECEF WGS84 points and UTM Orientation with ORTHOMETRIC HEIGHTS and not ELLIPSOIDAL HEIGHTS In CVS the Applanix model support ORTHMETRIC or ELLIPSOIDAL HEIGHTS where the ellipsoidal height must be above the WGS84 ellipsoid Usage applanix2ogeom camera file exterior orientation file lt imageToProcess gt optional output directory Where camera file an ossim formatted camera file for applanix exterior orientation file an Applanix file generated by their POSPAC software imageToProcess Is an applanix image Usually lt id gt extension The id is used to look up its orientation in the exterior orientation file optional output directory Output the result to an output directory If an lt imageToProcess gt is not given then all geometries from the exterior orientation file will be generated with lt id gt geom as the filename We currently only support either UTM or Ecef specified exterior orientation file The camera parame
12. build all overviews for foo tif and produce an output image called foo ovr 2 img2rr a foo tif Will copy the full resolution image into foo ovr and also build overviews into the same file Page 101 of 107 3 img2rr compression quality 70 compression type JPEG foo tif This will produce an overview file call foo ovr and the different resolutions will be jpeg compressed at 70 quality 10 4 create histo This application creates a histogram for the image and writes it to the sme filename but with extension his Usage create histo options image file entry akes an argument that specifies which entry to compute a min max for Sample Use create histo entry 0 foo ntf Will do the first entry of input file and create an external histogram called foo e0 his bins Override for the number of bins Really important for float data and higher bit depths By default it will be 65536 number of bins if can t be determined for the histogram Examples 1 create histo foo tif Will create a histogram for the input image and call it foo his 2 create histo min 100 max 1100 bins 4000 foo ntf Will create a histogram for values between 100 and 1100 and estimate it with 4000 bins 10 5 cmm Page 102 of 107 Compute min max pixel values from within a file Important to run this on radiometries other than 8 bit It will write an external OSSIM Metadata File omd file So if you have an image foo tif and you run cmm foo t
13. into the current directory with name lt id gt geom Page 106 of 107 2 applanix2ogeom camera file exterior orientation output dir Will generate all geometries into the output dir with name output dir id geom 3 applanix2ogeom camaera file exterior orientation imagefile Will generate one geometry file if the image is found in the exterior orientation file The image is usually something of the form lt id gt jpg or lt id gt tif etc It first strips the extension and then uses that for the id to lookup into the exterior orientation file 4 applanix2ogeom camaera file exterior orientation imagefile output dir Will generate one geometry file and then output to the output dir 10 8 igen igen is a powerful command line application that is entirely key wordlist driven It takes a keywordlist as an argument and builds a product from the keywordlist specification One way to use it is to setup a product in imagelnker and save out a spec file This is typically done if you have a product that is very large and might take several hours or maybe a day to run The spec file generated can be used to run the command line application igen that can be ran separate from the imagelinker GUI For example if you saved out the spec file from imagelinker to a file called igen spec then you can do igen igen spec producing a product based on your setting from imagelinker Adding mpirun np np igen igen spec will run the application on a cluste
14. loading All data loaded into ImageLinker will be stored into the global Layer Manager Image chain C data test_data sanfran sanfran cct nl xj Ele Edt Swipe Enhancements Help Fullres Propagate Refresh nearestneighbor ZA Figure 15 San Francisco image opened through the File gt Open Image menu option Image is part of the test data distribution under sanfran sanfran ccf Page 25 of 107 The display windows have a number of capabilities and menus that can manipulate the underlying image chain for the display Each display window contains a title bar that denotes the source Display menus and menu items operate on the underlying image chain for the window The window also has a number of modes for panning and zooming When a display is initially created it determines if there is an underlying geographic context associated with the imagery If none is available through associated metadata or internal geo specific tags a dialog is presented to the user asking if a default geographic projection should be assigned to the image 6 1 2 Open Project The Open command invokes a file navigation dialog allowing the user to navigate the file system and open a prj project file A project file saves the state of your session so that you can retrieve it later OSSIM also supports drag and drop of project files It will first close the current project and open the new project file 6 1 3 Close Project Closes
15. software library It is one of many programs that have been built on top of the functionality in OSSIM Users are encouraged to register online at http www ossim org and join the collaborative community there The site is the focal point for software development user support and related information ImageLinker is an application for viewing and processing remote sensing and Geographical Information System GIS data built on the underlying OSSIM library Capabilities include advanced remote sensing photogrammetry and image processing Underlying design goals and capabilities are covered in more detail in the Design Goals section of this document ImageLinker brings out many of the underlying functions through commands that are accessed through pull down menus and dialogs ImageLinker also allows the user to construct and manipulate complex image flows through the Visual Chain Editor VCE The VCE allows the user to dynamically connect and re order multiple image chains through the visual chain editor canvas Each function is represented by a thumbnail that reflects the state of the imagery at that point Users can double click any function to access custom dialogs and parameters Projects can be saved and retrieved at will and later used as processing templates for batch processing ImageLinker also supports a plugin architecture that allows new functions and processes to be loaded at runtime The design of the plugins allows scientists
16. std parallel 1 36 000000000000000 product projection std parallel 2 37 250000000000000 product projection tie point easting 2000000 000000000000000 product projection tie point northing 500000 000000000000000 product projection type ossimLambertConformalConicProjection Cutting the view out and paste it into a file called view kwl Now run the orthoigen with the view template orthoigen view template view kwl input tile1l jpg lambert jpg Produces an image called lambert jpg Figure 103 Using the view template we produced a lambert projected image See the difference between Figure 100 10 2 2Mosaicing Orthoigen can be used for mosaicing images together Note the layer ordering is as they appear on the command line The first image is the top layer and the last image is the bottom layer Running the following orthoigen will do a default A over B style mosaic Page 91 of 107 orthoigen t 256 utm tilel jpg tile2 jpg output tif i rE vw Figure 104 utm projected 256x256 thumbnail Another variation that allows one to change the mosaic type is to specify another combiner on the command line using the combiner type options orthoigen combiner type ossimClosestToCenterCombiner t 256 utm tilel jpg tile2 jpg output tif orthoigen combiner type ossimFeatherMosaic t 256 utm tilel jpg tile2 jpg output tif E phr Page 92 of 107 Figure 105 utm projected 256x256 thumbnail doing a closest t
17. the target histogram to be the first layer Once set we do not have editor in place to change it For an example there are four files supplied in the test data suite under doqq overlapping directory Load all four images into the system Figure 59 Four overlapping doqqs found under the doqq overlapping directory of the test data suite Make sure you go ahead and propagate the view from one display to all other displays This will guarantee all are on the same view plane Next select the Layer gt Histogram Match menu item and it should bring up a selection box and select all 4 images 8 Non Commercial Choose Layers to Match 32 xl Cancel Figure 60 Histogram matching layer selection When apply is hit there will be 4 new layers added to the layer manager that have two equalization filter with the second filter being the target histogram The default target histogram is the first layer and we currently do not have a GUI to Page 53 of 107 edit the target histogram Once the four new chains are created they are mosaiced automatically together with a feather combiner Histogram match mosaic gt jol xj Eile Edit Swipe Enhancements Help Zoom Pan Full res Propagate Refresh nearestneighbor gt 34 54 03 8765 N 076 48 17 9573 W 39 0946 m rgb 36 28 11 0s File Edit Swipe Enhancements Help Zoom Pan Full res Propagate Refresh nearest neighbor v 34 54 03 8765 N0
18. to perform the operation Page 85 of 107 ELLEN MUN L EE annotate akes an annotation keywordlist chain template Specify a chain template for each image to process combiner type can be ossim class names For example To get actory dump t ossimImageCombiner type Here is a list of combiners tested ossimBlendMosaic ossimMaxMosaic ossimImageMosaic ossimClosestToCenterCombiner ossimBandMergeSource ossimFeatherMosaic cut center ll akes two arguments in degrees specify the latitude followed by the longitude of the cut center For example cut center ll 28 5 118 will specify the cut center at location lt 28 5 118 gt cut radius meters Specify the cut distance in meters This will ake the center cut point and produce a bounding box surrounding the center point with he specified distance enable entry decoding A filename can be separated by a anda number NO space Example a toc O0 will do entry O of a toc file and on unix you might want to use a V since the separator is a pipe sign Example a toc 0 will do entry 0 of a toc Will set the view to a geographic projection with origin at 0 0 The GSD to the input akes a latitude as an argument for purpose of scaling Specifies that no spec file was defined Defaults to a scaled geographic image chain with GSD to the input Specifies a utm projection Specify in normalized percent the low clip and hen the high clip value Note You must
19. 76 48 17 9573 W 39 0946 m rgb 37 30 13 8l x Figure 61 Histogram Match and feathered shown on the left Original mosaic with no matching and feathering shown on the right 7 5 3 Combine Histogram match Mosaic Fusion r Blend Correction Feather Elevation Merge bands Ortho Max Closest to center Factory Figure 62 Combine Sub menu option Page 54 of 107 The combine sub menu holds menu items that invoke various image chain combiner methods The Factory menu item currently does nothing in the current release of OSSIM When selecting any of the combiner types it will bring up the layer chooser window with the title corresponding to the combiner just chosen This chooser was from the Layer gt Combine gt Mosaic menu option A Non Commercial Choose Layers to Mosaic 10 Image chain Cjdataltest data sanfran sanfran ccf 35 Image chain C data test data sanfran map highres utm map tif Cancel Figure 63 Chooser for any combiner selected Brings up all layers in the Layer manager A couple of examples of different mosaics Figure 64 Feather combine left simple combine middle and a blend combine Merge Bands combiner will merge all input bands into a single layer Ortho combiner is used if all input data was originally projected to the same view and all that is required is a simple shift to put them together For example if al
20. 77 N 12207019 1900 W 32 1553m rgb 223 220 211 A 5 Figure 56 View and resampler propagated to all open displays 7 5 Layer Menu Layer Vce Utilities Wind Import image file Histogram match Combine Fusion Correction Elevation Figure 57 Layer menu found on the Main ImageLinker window Page 51 of 107 The Layer menu allows creation of other layers by using the loaded data as input Layers can also be used in other operations to create additional layers All layers are visible throught the Layer gt Manager menu option 7 5 1 Manager A Non Commercial Data Manager A 0 Image chain C data test_data sanfran sanfran ccf i Display 35 Image chain C data test data sanfran map highres utm map tif Delete Duplicate Mosaic Blend Hr ER Feather Figure 58 Layer Manager brought up through the Layer gt Manager menu option Displays all the current layers in the system We have added some quick action buttons on the right hand side indicated by Display Delete Duplicate Mosaic Blend and Feather buttons Display button will open displays for any selected image If multiple layers are selected then a display is opened for each layer Delete button will delete a layer in the system If a display is open for that layer it will automatically close if it is the only layer for that display Page 52 of 107 7 5 2 Histogram Match Currently the matching automatically selects
21. IM project joined the OSGF as one of the founding projects Future efforts will integrate these projects into formidable open source image processing mapping OSGeo org Your Open Source Compass and GIS systems Page 5 of 107 3 Quick Start Tutorial One of the best ways to understand the power of the OSSIM library is to demonstrate it through the ImageLinker application ImageLinker was originally developed as a prototyping tool to test and demonstrate the power of the OSSIM library It has a graphical user interface implemented using the Trolltech Qt library and runs on any platform supported by them Windows Linux Solaris Mac OSX Figure 1 Combiners in the Quickstart Tutorial This quick walk through will demonstrate just a few of the capabilities of OSSIM and ImageLinker It is designed so that the user can follow along with step by step instructions using the available sample data sets The tutorial will demonstrate different ways to combine a map and satellite image OSSIM and ImageLinker provide several ways to combine geo spatial data sets This exercise will demonstrate simple mosaics blends and feathers with a map and image of the San Francisco area The data sets in this example are part of the sample data set that is available on ossim telascience org ossimdata see section 3 2 for download instructions Page 6 of 107 Refer to the Installing sectio
22. NOTE this is only used in igen application and not orthoigen This mask is actually specified Page 95 of 107 on the command line as the output image if tiling template option was specified The mask variables are 9er9e all occurances are replaced by the row number C all occurrences are replaced by the column number i all occurrences are replaced by the tile id or all occurrences are replaced by the origin of the row oc all occurrences are replaced by the origin of the col SRTM does SRTM tiling scheme Nice to use when building a database of SRTM elevation from some other elevation format and projection igen tiling padding size in pixels Specifies a padding to use It will extend the borders out by the specified number of pixels in the horizontal and vertical direction Pixel tiling is used most commonly when you want to guarantee some power of two set of images or the single image is too large and needs to be split apart We will tile up tile1 jpg image in Figure 100 to a 64x64 set of tiles First define the template tiling keyword and assume the name is pixel tiling kwl igen tiling type ossimTiling igen tiling tiling distance 64 64 igen tiling tiling distance type pixels igen tiling delta 1 1 igen tiling delta type per pixel igen tiling padding size in pixels 0 0 Now execute orthoigen utm tiling template pixel tiling kwl input tile1 jpg pixel 96r96 96c96 jpg NOTE we put the file name mask as
23. OSSIM Users Guide Version 1 6 5 Open Source Software Image Map OSSIM pronounced Awesome Document version 1 1 July 10 2006 Written by Garrett Potts Mark Lucas Page 1 of 107 About OS SIM icta eerte ere rhe e ERREUR RNEREETE NE E Sa E tease bee ense rre eub S Quick Start Tutorial eerte oh tore iab Open the source files 8 9 1 eeeseeeessssesssoseescoseeccsosssssossesseeece Managing Windows eee ee esee e eese ee eene te einen etas tne ense tne ta senses eto eene ases etos tns en senses etn sene enae Pan Zoom and Propagate in the Displays eeeeeeeeee eene entente enses ens enessessenessese Save your project 16 Minibar CERAM LO SITET a EE V Sample Data m eeesseseessssesssoseeessoseessossecsssssecsossessoseesssoeeco Source Code for Developers eere eee eee eese eee eese ee eene en sens esses s eneensesseeseeseeseessessessees LD Preferences Setups cssssssscccsssciesscestensesseccesdesssssceasesccoseesseseassestevesancusdedsocsscesssdessesssesneceaicescsescesese LO Elevati n Setup zxismodiefee o d Hentai non IO ORO eere et 19 Sample ossim preferences file Geoid Setup OSSIM Planetary Viewer Setup osgPlanet Plugin setup ImageLInker ioi e RETE ee ie PV Re te A tmd File Menu eee 25 Open Image eesseeeeccssececcoseecoe
24. OSSIM and ImageLinker Using a map and image of San Francisco we demonstrated basic navigation of displays and window management tools in ImageLinker followed by the creation of various types of mosaics and overlap area combiners These types of products are very easy to create in the tool sets OSSIM automatically handles differences in input resolution projections file formats and radiometry types Transformations are performed automatically to produce the product needed by any given display Adjustable parameters can be accessed and modified through the Layer gt Edit menu of any display Covered in more detail later any display can generate a product through the File gt save as command Page 16 of 107 4 Installing This section describes how to download and install the OSSIM software OSSIM is hosted at www ossim org There you will find download links for the source code and and binaries If you are only interested in running the software then download the binary installer for your platform If you are a developer and want to build the software yourself you have the option of downloading the latest stable release version from the download site or checking the latest repository from the CVS tree 4 1 Binaries The easiest way to get up and running is to download the binary installer for your platform double click and follow the instructions The binary package installers are available for many of the common platforms The binary packa
25. al properties can be edited by selecting the edit writer button compression quality compression type pixel type lut file color lut flag Class name assimTiffWriter Description enabled True filename create envi hdr False create external geometry False create fgdc False create jpeg world file False create readme False create tiff world file False create image True create overview True create histogram True V Auto apply Reset Close Figure 44 Writer property list Page 43 of 107 If the GDAL plugin is enabled then you can select other writer types related to GDAL For example gdal HFA will allow you to output an imagine formatted file Fl Non Commercial Property List 3 x Class name lossimGdalWriter bu Description enabled True filename create_envi_hdr False create external geometry False create fgdc False create jpeg world file False create readme False create tiff world file False create image True create overview m True create histogram True gdal overview type none HFA_USE_RRD False HFA_BLOCKSIZE 0 HFA USE SPILL False HFA COMPRESSED False Iv Auto apply Apply Reset Close Figure 45 Writer property list for the gdal_HFA output type To see the supporting formats of the GDAL plugin see Main Menu gt File gt preferences or under MACS see Main Menu gt ImageLinker gt Preferences
26. and developers to quickly add new functionality without the need of learning the internals of the OSSIM software libraries In addition to viewing and manipulating image chains ImageLinker can generate output products by passing the parameters to the igen command line tool called through menu commands ImageLinker and OSSIM are open source software projects that are being widely deployed in various government and commercial organizations Page 24 of 107 6 1 File Menu E Non Commercial unnamedo Ele Layer vce Utlities Window Help Open Image Ctrleo Open Project Close Project New Project CtrieN Save Project Ctl S Save Project As Preferences Ext Figure 14 File menu from Windows version of ImageLinker The project menus contains all of the commands necessary to loading saving closing and creation of new projects Project files can be thought of as sessions Multiple images can be loaded and filters and parameters can be applied to their respective image chains All of these settings including window positions and user settings are captured and saved at will when the user saves a project A project typically has a prj extension 6 1 1 Open Image File gt Open Image will bring up a file browser dialog This is a multi selection file browser allowing one to select more than one image to import Drag and drop is supported by ImageLinker and is another alternative to the file browsing option for data
27. at area of interest and if the shift is used then it will zoom out When in the pan mode clicking in the window will re center the view to the clicked position This is not the case if the image is already fit to the current view window 6 2 2 Fit The Fit button calculates and scales the image so that the entire coverage is visible in the display window as shown above in Figure 15 In this example a Landsat image has been scaled to fit in the display windows 3T 4923 5074 N 12292212 6820 W 32 154 m roh E 76 73 37 43126 9716 N 192020111 9857 W 32 1539 m rgb 209 206 205 Figure 17 Click on treasure island in the left image and then click the full res button The Full res button scales the image to full resolution centered about the last cursor click in the window Each geographic image has a full resolution parameter which defines the native resolution mode of the data set In this example the Landsat image on the left was zoomed to full resolution Page 27 of 107 AT del 37 49128 9716 M 122922119957 W 32 1539 m rgb 208 206 203 3T 4923 8074 N 12202212 6820 W 32 154 m roh 65 76 73 Figure 18 Propagate the view of the left image to all other opened displays The propagate mode syncs the views of coincident geographic data sets to the same scale and resolution as the selected window It transmits a command to the other windows to
28. dee 37 49 28 6189 N 122 22 16 9625 W 32 1553 m rgb Y 37249294760 N 12292215 3360 W 32 1547 m rgb 1 4 Figure 10 Blended layers Page 13 of 107 Tile the windows by pressing the menu option Window gt Tile and in the Landsat image window press the Propagate button The view is propagated to the blend and map windows Non Commercial unnamed0 prj 5 x File Layer Vce Utilities window Help Image chain C data test_d joj lt image chain c data test lolx File Edit Swipe Enhancements Help Ele Edit Swipe Enhancements Help Zoom Pan zem Pan Full res Propagate Refresh it Full res 37 49 34 7901 N 122 22 18 1374 W 32 1547 m rgb D NCIITLICONENENENSSSSS 8 Ele Edit Swipe Enhancements Help Zoom Pan Fullres 37 49 34 7901 N 122 22 18 1374 W 32 1547 m rgb d 3 Figure 11 Blended layers tiled and propagated Minimize the source windows for the map and Landsat image by clicking in their display minimize buttons This should leave only the blended display Now we will create feathered and mosaic displays Layer Combine gt Mosaic and select just the map and Landsat image chains E Non Commercial Choose Layers to Mosaic Figure 12 Mosaic layer chooser Page 14 of 107 Press the Apply button Next we will go ahead and create a feathered display as well In similar fashion select the menu option Layer
29. e object2 prefix and save in the writer template kwl color lut flag 0 compression quality 75 compression type none create envi hdr 0 create external geometry 0 create fgdc 0 create histogram 1 create image 1 create jpeg world file O0 create overview 1 create readme Q create tiff world file O0 description enabled 1 filename data ear1 output docs writer tif id 33 image type tiff tiled band separate input connection1 9 input list fixed 1 Page 99 of 107 number inputs 1 number outputs output geotiff flag 1 output list fixed 0 output tile size x 64 output tile size y 64 overview compression quality 75 overview compression type 1 pixel type point projection units meters type ossimTiffWriter Now execute the template on the input image orthoigen utm writer tempLate writer template kwl tile1 jpg output tif It will use the settings in the writer template kwl to instantiate the writer and output the product 10 2 image info This application displays the general information about an image Usage image info options input file By default all information is output for the input file If multiple entries exists then each is prefixed with image lt number gt where number is the entry number palette 1 fone exists a color palette will be outputted Only image information such as height idth etc is outputted Output the information to the specifi
30. ec 74 Elevation Managetz eni erre eere tre Serene hone ie M Se Ne ene e ee SEVERE E PIER S ee LEM E a 74 Unit CONVE ET i oes RU RE RHET e He MERERETUR NUR RISUS DD MM R 75 buaunn diu M X Q 1O OSSPIONC Rm ati aiii iiei ari ieda itidir ii H OSSIM NITE gt Dynamic Image Chains sccssscsscsssssescsseccsscessscescessssescessssescessesessessscessessssessessssessessesesees 7 7 File formats supported eere eese e esten ee eene enean sensns niesen assieto ii a 70 Projections and Transformations Supported eere ecce eere eee eee senes ees 79 Command line Applications esee eee eese essen seen nn eene senses eene enses ee OO OVCHOIGEM seccisessacceedssnastecnessasonsoccssentededsedsoasacasconeccosstssadssbesassasseasseasosnedwesoseashscasenssveasonsacsbescsoross OD R ptojectiOni ae enenho epp iem D E te 89 Mo SaiCing iss D m 91 Histogram iodin 93 Cut and Resample 55 iei eene HET HE GER Ce Frei e OR SERTE TENS 94 im ge mT PLI ridic L2 Page 3 of 107 senes 104 seese LOS eeessesesssossessosesssoseeesssseessssseccossesssseeesso sesso ogeom2ogeom applanix2ogeom mnm mM mn eeeesesessssssessose esses ee essossessssse
31. ed file Specify the datum you want the ground rect returned Currently it is the datum of the projector Overwrite output geometry file If the input image has any then output the projection information Will also output orner coordinates Examples Page 100 of 107 1 image info file name This will output all information to standard out 2 image info palette file name Will output a palette only if one exists 3 image info i file name Will output image information only 4 image info i p file name Will output image and projection information 10 3 img2rr This is an overview builder application It is important to build overviews for any image file to be processed at different resolutions When executed on any image a multi resolution pyramid is created for the image as an external file with a ovr extension In the case of r multi entry it will ad _e lt number gt ovr to each entry a or include Will copy the full res image into the fullres overview file e or entry If the input image is multi entry then do a specific entry resample type Currently you can do nearest or box Box is the default and is just a 2x2 average Currently you can do nearest or box Box compression quality is the default and is just a 2x2 average compression type The compression type can be NONE JPEG PACKBITS or DEFLATE The default is none and no compression is used Examples 1 img2rr foo tif Will
32. ee a display similar to that shown above Note the tracking cursor lines in the map display track the geographic location of your cursor Tracking cursors will appear in all geographically coincident windows Non commercial unnamedo prj loj x File Layer Vce Utilities Window Help Image chain C data test_data sanman EI EI _ Image ch IES File Edit Swipe Enhancements Help Fie Edit Swipe Enhancements Help Zoom Pan Zoom Pan Fulres Propagate Refresh Fullres Propagate 37 49 32 4832 N 122 22 20 2605 W 32 1557 m rgb amp 37 49 32 4083 N 122 22 18 9834 W 32 1553 m rgb 17 A Figure 6 Full resolution of the Landsat San Francisco Note Depending on the data set you are using your view may differ slightly The test_data tgz data set has been downsampled and therefore will display at a 15m versus 5m resolution when zooming to full resolution Page 11 of 107 In the Image Window Press the Propagate Button The geometric view is updated in all open windows to the same scale and orientation You can toggle the cursor modes between zoom and Pan with the icon buttons at the top of the display for general navigation A Non Commercial unnamedo pri lol x File Layer Vce Utilities Window Help E o moe a jai Ele Edt Swipe Enhancements Help Ele Edit Swipe Ephancements Help Zeom Pan Zoom Pan Fullres Propagate Fit Fu
33. ee section 4 7 2 For using the component window in the Canvas window and connecting individual filters together Page 64 of 107 7 3 2 Chain Editor Window The chain editor Window is a visual canvas for building process flow from the ground up Currently the chain editor is not directly linked o the global data manager but does support dragging and dropping individual chains from the Layer gt Manager window to the VCE canvas window From the main window hit Vce gt New chain menu option and stretch the canvas out to give a little more workspace Ele Edit X Sg amp amp ET 3 a Figure 80 VCE Canvas window from Vce gt New chain menu option 7 3 2 1 Building A Chain From Scratch In this example we build a re projection chain from scratch At the minimum we will need an image loader and an ossimlmageRenderer filter from Vce gt Components dialog First bring up the components dialog and select the loaders tab Page 65 of 107 Non Commercial Voe Components Combiners Fiters Exporters Displays ImegeHendler E Figure 81 VCE component window with the loaders tab selected Select the ImageHandler component with the left mouse button While holding down the left mouse button drag the ImageHandler icon onto the canvas This will immediately cause a file browser native to your platform to open
34. elevation srtm 1arc Load both cells into the system Figure 72 SRTM elevation cells Create an ortho mosaic through the Layer gt Combiner Ortho menu option and select both SRTM layers ossimOrthoImageMosaic Ele Edit Swipe Enhancements Help Figure 73 SRTM Ortho mosaic Page 60 of 107 The hillshading portion is a 2 step process and future releases will merge these into one For now we need to first apply surface normals to the ortho mosaic so it can be used in a shading equation Select Layer gt Elevation gt Normals menu option 2 Non commercial Choose Layers to Apply Normal Calculations 328 Image chain C data test data elevation srtm 1arc N37W123 hgt 364 Image chain C data test data elevation srtm 1arc N38W123 hgt 377 Image Ortho chain Ortho chain data test data elevation srtm 1arc N37W123 hgt 380 Image Ortho chain Ortho chain data test data elevation srtm 1arc N38wW123 hat 403 0ssimOrtholmag ai 436 lmage chain C data test data sanfran map highres utm map tif cm Figure 74 SRTM Normal calculations Select the ortho mosaic chain and hit the apply button Now Select Layer gt Elevation gt Hillshade E Non Commercial Choose Layers to Hill Shade 328 Image chain C data test data elevation srtm 1arc N37W123 hgt 364Image chain C data test data elevation srtm 1arc N38W123 hgt 377 Image Ortho chain Ortho chain da
35. ence between the multi spectral and the high resolution grey scale image The larger the scale difference the larger the kernel needs to be The second parameter applies a gain to the high resolution pan Currently the gain is not adaptive Figure 70 5 zoomed to 1 meter image using cubic resampling left Final SFIM sharpened image with blurring kernel set to 5 using the 1 mter grey scale image for detail preservation The local correlation fusion has the following parameters that can be adjusted E Non Commercial Parameter Adjustments Adjustment Selection 0 bd Adjustment description Regression clamp High pass gain Blurring kernel width JV Auto apply Keep Copy Delete Reset Close Figure 71 Fusion parameter edits for the Local Correlation Algorithm The attenuator parameter is a global attenuation on how much information is added back into the multi spectral dataset The Regression clamp clamps the slope of the regression equation Allows one to clamp large slopes Page 59 of 107 The High Pass Gain runs a high pass over the high resolution image before applying the fusion The Blurring kernel width is increased based on the scale change The larger the scale difference is between the multi spectral and the grey scale image then the larger the blurring kernel value should be 7 5 5 Hill shading In the test data suite we have supplied wo SRTM elevation cells under
36. er 11 37 92626 122 37761 cut radius meters 50 meters 5 input tilel jpg sinc jpg Page 94 of 107 Figure 107 Shows a cut and resample differences These are tiles cut at lat 37 92626 and lon 122 37761 with cut radius of 50 meters and gsd resolution of 5 meter The original was a 5 meter image giving a 10x zoom Left is nearest neighbor resampling middle is cubic and the right side is sinc See options table for other resampler types 10 2 5Tiling Orthoigen supports geographic meters and pixel tiling options The tiling specification is currently through an external keyword list and is loaded using the tiling template command option The format of the template is igen tiling type ossimTiling igen tiling tiling distance 10000 10000 igen tiling tiling distance type meters igen tiling delta 100 100 igen tiling delta type per pixel igen tiling tile name mask tile 6r 6 96c96 igen tiling padding size in pixels 0 0 where igen tiling tiling distance specifies the unit of the distance type igen tiling tiling distance type can be pixel meters degrees minutes or seconds igen tiling delta specifies the delta distance If delta type is per pixel then it is in the unit of the distance type If the delta type is total pixels hen the delta is the total number of pixels within the tiling distance igen tiling delta type can be total pixels or per pixel igen tiling tile name mask specifes a mask for the filename
37. ersion tool is accessed through the Utilities gt Unit converter menu option Page 75 of 107 Meters NEN Feet 228083989500 U S survey feet 320083333335 Miles fo o00621371192237 Nautical miles 0 000538989170472 Radians 1 56785594289e 007 Degrees f8 9831528412e 006 Minutes 0 000538989170472 Seconds 0 0323393502283 O Dms 00000000323 Close Figure 98 Unit converter dialog interface To use the converter type a value in any of the fields above and then tab out of the field 7 5 Window Menu Window Help Cascade Tile Minimize all Restore all Maximize Refresh all Close all Figure 99 Window management menu The window menu allows for managing the layouts and group closing of all displays in the project Window gt Cascade will cascade all open or minimized displays Window gt Tile will grid the and fit all displays onto the main window Window gt Minimize all menu option minimizes all shown displays Window gt Restore all menu will restore all minimized windows Page 76 of 107 Window gt Maximize menu option causes the current display to be maximized Window gt Refresh all menu option causes all open displays to be refreshed Window gt Close all Closes all open displays Note The data is still loaded in OSSIM and can be accessed in the Layer gt Manager menu option To delete a layer from memory this can only occur in the layer manager 8 osgPlanet see t
38. filter from the Components dialog s filters tab onto the canvas Drag an EdgeFilter to the canvas Page 71 of 107 8 Non Commercial vce components xl Loaders Combiners Filters Exporters Displays 3x3ConvolutionF AtCorrGridRema AtCorRemappe BandAverageFil ilter pper r ter BandClipFilter BandSelector BrightnessContr BrightnessMatch astSource CacheTileSourc CastTileSource ConvolutionFilter ConvolutionSour e Filter 1D ce EastingNorthing EsriShapeFileC EsriShapeFileFil FiFilter Cutter utter ter GeoAnnotationS GeoPolyCutter GridRemapSour HistogramEqual ource ce ization HistogramMatch HistogramRema HistogramThres HsiRemapper Filter pper hhaldFilter HsiToRgbSourc HsvToRgbSour ImageGaussian ImageRenderer e ce Filter ImageSharpenF ImageToPlane IndexToRgbLut JpegYCbCrToR gt Figure 92 EdgeFilter under the Filters tab P Non Commercial unnamedo prj Ele Layer Vce Utilities Window Help Visual Chain Editor Bl x eere al Os ail A Figure 93 EdgeFilter dragged from the component filters window to the canvas Now cut the link from the left most filter and the second filter and connect them to the new edge filter Page 72 of 107 Non commercial unnamedo prj Ele Layer Vce Utiles Window Help Visual Chain Editor laxi Ele Edt Ra E113 m al v kMenepepepee 7 Figure 94 Final edge filter
39. future use or for use by external OSSIM applications By default it should fill the correct lt filename gt geom file in the file browser dialog The close button will close the dialog box Page 30 of 107 7 1 2 Image Chain In the image window selecting the menu option Edit gt Image Chain will bring up the image chain editor m Non Commercial Image Chain Editor 2 xl Image Layer TENTS RIG Eam TERUP Ie s re er MT nr onem E Filter Action Available Filters Image Chain ccf 3x3ConvolutionFilter Insert Before AtCorrGridRemapper Band Selector AtCorRemapper Insert After Histogram Remapper AtbController Tile Cache BandAverageFilter Delete Image Renderer BandClipFilter Tile Cache Edit Properties close gt BandMergeSource BandSelector Rlanrh4neaic didi Figure 21 Image window image chain The left are filters that you can add to the current chain and on the right is the current processing chain for the window OSSIM supports dynamic image chains This will become more apparent later in the demonstration with the visual chain editor see section 4 7 An image chain typically begins with a geo spatial file on disk and is connected through a series of band selectors caches filters and resamplers This dialog shows the current image chain right side associated with the particular display The left hand list contains all the registered filters
40. gers will typically install the OSSIM library in the proper system location command line utilities and the ImageLinker and lview applications Most users will want to start with the ImageLinker application to get a better feel for the capabilities of the OSSIM package 4 2 Sample Data You can download the test data at ossim telascience org ossimdata demodata tgz full rez large map and image test data tgz test data as shown below Some of the directories you should see listed in test data tgz are doqq doqq_overlapping drg elevation geoid03 geoid96 geoid99 katrina orthoigen sanfran sanfran map These test sets are used throughout the users guide Each directory is explained Page 17 of 107 doqq directory is used to demonstrate the fusion algorithms see section 4 6 4 Fusion below You should see a 5 meter color image and a 1 meter grey scale image of the same area This is used to simulate a 1 meter high resolution color product by using the 1 meter grey scale image for detail enhancement doqq overlapping directory has 4 tiles that overlap each other This is used to show different mosaics drg directory holds an image without geometry This is used to show how to drop ground points to create an ossim projection see section 4 2 9 elevation directory holds a couple 30 meter cells SRTM cells geoid03 directory holding NGS distributed geoid 2003 grids geoid99 directory holding NGS distributed geoid 1999 grids ge
41. h field can be edited and changed The slider is there for convenience and if you need to have more precision then you need to hand edit the field in the Page 29 of 107 param column The bottom of the dialog we have the buttons apply keep copy delete reset save and close The apply button will apply the current adjustments to the image This is not needed if the auto apply button is enabled If the auto apply is enabled then the image changes immediately when any parameter changes This is true by default The keep button will make a copy and re center the adjustments So if you move the top slider over and hit keep button it will make a copy of the adjustment adjust the center to that value reset the param back to 0 center center param sigma param 0 You can go back to the previous adjustment by selecting the adjustment number in the adjustment selection at the top of the dialog The copy button simply makes a copy of the current adjustment The previous adjustment can be retrieved through the adjustment selection list at the top of the dialog The delete button will delete the current adjustment If the current adjustment is the first adjustment it will simply reset the adjustment to the initial condition The reset button will reset the current adjustment to the sensors initial adjustment values The save button will write the adjustments to disk This is important if you want your adjustments to stick around for
42. hat can take any projection understood by OSSIM and produce a coarse grid or RPC projection It is important to note that the application will currently overwrite the passed in geometry if it was originally an ossim geometry file with a geom extension Usage ogeom2ogeom options lt geometrylimagefile gt The argument to the applications can take either an image file and try to extract the geometry and convert it or an ossim geometry file Whatever file is given it will produce that file with a changed extension to geom Some of the options are Options Description Create a coarse grid projection This will estimate the input projection with a grid pplies to cg option This will disable the creation of adjustable parameters in the coarse grid The coarse grid generates a layer for every adjustable parameter taking up a lot of Space This can decrease the amount of space used by the coarse grid If the output is RPC then all heights are O For coarse grid it just turns off applying heights to the model akes 4 arguments ulx uly width height If you give the application a geometry file only then you must specify the image coordinates for the geometry rpc Generate an rpc output model Page 104 of 107 Used as an error tolerance for the estimate his is mainly used by the coarse grid for splitting to higher estimates if it doesn t fall within tolerance Examples 1 ogeom2ogeom cg image Will
43. have ran create histo on all images before applying his operation input proj he output projection will default to the input projections Page 86 of 107 meters Change the GSD to the passed in meters per So if you want 50 meters per pixel GSD hen make this meters 50 resample type Specify what type of resampling you would like alues can be nearest bilinear cubic gaussian hanning hamming lanczos mitchell catrom blackman sinc qessel quadratic hermite bspline slave buffers Specifies number of tiles to buffer when doing MPI applications When passing the tile from he slave node back to master we can continue processing the next tile if the buffers is larger han 1 For performance you can play with this value If not present then by default it will be 2 buffers per node tiling template Specify an external file that contains tiling information A template can be passed to orthoigen that enables tiling support Tiling ican be done in pixels meters decimal degrees For meters and decimal degrees tiling it will depend on the view projection If you are geographic then use decimal degrees else use meters view template Is an OSSIM keywordlist the describes the view projection and parameters writer template By default a writer is determined based on the extension of the output file A writer keywordlist can be used a template to orthoigen if you need more settings than the default action taken based on
44. have 3 one band input images then you will have a single 3 band output image Since OSSIM is factory driven for a complete list you can execute the application factory dump t ossimlmageCombiner and it will list all ossimlmageCombiner type objects 10 2 3Histogram Histogram operations can be applied to any image that has histograms built Histograms can be built by using the create_histo application Currently we have only exposed a linear stretch given a penetration percent The clip percents are normalized between 0 and 1 orthoigen hist stretch 1 1 utm tilel jpg histogram tif Page 93 of 107 Penetration is given as the first number corresponding to the left side and the second number cooresponding to the right side Figure 106 Histogram stretched utm image 1096 1 penetration on both the left and right side 10 2 4Cut and Resample The filter kernel can be changed by setting the filter type using the resample type We will zoom tile1 jpg to 5 meter The original is a 5 meter image and then cut it by 50 meters Here we show 3 command line executions with nearest cubic and sinc respectively orthoigen resample type nearest utm cut center ll 37 92626 122 37761 cut radius meters 50 meters 5 input tilel jpg nerest jpg orthoigen resample type cubic utm cut center ll 37 92626 122 37761 cut radius meters 50 meters 5 input tilel jpg cubic jpg orthoigen resample type sinc utm cut cent
45. have large NULL areas since the image is rotated The valid vertices can specify the valid corner points This is currently for visual purposes only and cannot be edited here Page 35 of 107 7 1 4 Drop Point GUI If you would like to follow the steps here please remove the file test data drg battlmn geom If this is not the first time running this part and you want to see how to do a drop point projection then remove the file test data drg battlemn geom Next load the image test data drg battlmn tif and say yes to create a default projection Figure 28 Load image test data drg battlemn tif The loaded image should look like this Ini xj Eie Layer Vce Utilities Window Help Image chain C data test_data drg battlemn tif Ble Edit Swipe Enhancements Help Figure 29 Loaded image Page 36 of 107 Zoom to full resolution to the upper left corner of the image Non Commercial unnamedO prj Ble Layer Vce Utilities Window Help Image chain C data test_data drg battlemn tif Ele Edit Swipe Enhancements Hep Zeom Pan Fk Fullres Propagate Refresh nearestneighbor BATTLE MO mor oe R4AE ELK T36N I m xil 00700 08 4921 S 000 00 16 6305 E 17 1608 m rgb 255 255 255 m Non Commercial unnamedO prj ijs x File Layer Vce Utilities Window Help Image chain C data test data drg battlemn tif He Edt Swipe Enhancements Help
46. he geometry information into the output file Typing icp return will display the help screen followed by a list of supported output types If the gdal plugin is enabled then you will see a number of gdal supported output types prefixed with gdal At the minimum you must specify and output format and input file and an output file Usage icp options output format input file output file a or use scalar remapper Will use a scalar remapper to transform to an 8 bit image b or bands akes a quoted comma separated list of bands o select For example bands 1 2 4 Will create a 3 band output where 1 is the first followed by 2 as the second and followed by 4 as the third c or compression type akes an argument that describes the compression type Currently this argument is only valid for tiff output types The argument value can be jpeg packbits deflate or zip akes an argument that specifies which entry to extract Tis argument should only be used for multi entry input sources Sample Use entry 0 Page 82 of 107 ou can specify an offset in the line direction hat specifies where to start the extraction By default this is 0 Sample Use start line 100 Will start at line 100 ou can specify an offset in the line direction on where to stop the extraction By default this is ou can specify the offset in the sample direction on where to start the extraction By default this is 0 Sample use
47. he output of the first row Figure 109 These are the fist row of tiles produced from the orthoigen meter tiling example Let s specify the tiling in seconds and put it in a file called seconds tiling Kwl will create a tiling for the 5 meter UTM file in Figure 100 Now since we are a geographic tiled output we will need to convert the input image into a geographic projected image using the geo option for the output view The tiling keywords are Page 97 of 107 igen tiling type ossimTiling igen tiling tiling distance 10 10 igen tiling tiling distance type seconds igen tiling delta 16 igen tiling delta type per pixel igen tiling padding size in pixels 0 0 This says to do a 10x10 seconds image where each pixel size is 16 seconds Now execute the command line orthoigen geo tiling template seconds tiling kwl input tilel jpg seconds_ r _ c b jpg and we get tiles named seconds 0 O0 jpg seconds 0 1 jpg etc Using Figure 100 as the input we get Figure 110 These are the fist row of tiles produced from the orthoigen meter tiling example 10 2 6Writer template Writer templates are currently available with CVS builds since the version 1 6 4 release If running from the CVS source code you can apply the writer template to specify specific output attributes and change to other writers The easiest way to create a template is to use imagelinker to create an igen spec file much like was done for the view templa
48. he separate tutorial on osgPlanet osgPlanet is a 3D planetary viewer that is built on top of ossim and OpenSceneGraph 9 OSSIM Internals Open Source Software Image Map OSSIM is a cross platform open source software distribution of remote sensing image processing and geo spatial mapping technologies The OSSIM distribution is centered around a C object oriented software library Key functionality is implemented in the library for use in command line tools GUI applications and web based services The software distribution includes utility programs image viewers and sophisticated custom production applications As improvements are made to the core library all applications and services that use the library benefit OSSIM can load and process a wide range of geo spatial and image formats It supports rigorous sensor models map projected products and Residual Polynomial Coefficient RPC mechanisms Most government and commercial formats are supported by OSSIM When geospatial data is loaded the associated meta data is typically processed to correctly map project and provide precision terrain correction over elevation data where appropriate OSSIM understands and handles map projection and datum transformations and resampling to arbitrary resolutions 9 1 Dynamic Image Chains EEEESH aes SESE SS Basic to OSSIM is the support of Dynamic Image Chains The user can dynamically connect loaders combiners filters and outputs within a
49. if then the file foo omd is written updated to the new values Usage cmm image file Under unix type shells you should be able to wildcard it since they get expanded cmm tif Executing the command within a shell that supports wildcarding like a bash shell would compute min max for all files with a tif extension This program will force a re compute of Min Max Pixel Values Options Description min Min pixel override The min value is not computed rom the image but uses this value as the min value or all bands max Max pixel override The max value is not computed rom the image but uses this value as the max alue for all bands nul an override the null value e or entry akes an argument that specifies which entry to ompute a min max for Sample Use mm entry 0 foo ntf ill do the first entry of input file l or list entries ill list the id indices of each entry Note Since ome handlers have multiple le entries where some entries can t be rendered either they aren t renderable or we just don t support it yet you ould get non sequential numbering For example ou might have an image with 10 entries but only 3 Page 103 of 107 Examples 1 cmm foo tif Min and max values are computed 2 cmm min 1 max 200 foo tif The file is not scanned since both min and max are overridden and will be Written to the foo omd 10 6 ogeom20geom ogeom2ogeom is a geometry conversion application t
50. in OSSIM Using the San Francisco image in Figure 15 we can add an edge detector to the chain by first selecting the EdgeFilter in the available filters list and then select the ccf filter on the image chain list Page 31 of 107 m Non Commercial Image Chain Editor Image Layer layer 1 image chain C data ccf_utm foo ccf Available Filters ConvolutionFilterl D ConvolutionSource EastingNorthingCutter mj ElevimaqeSource EquationCombiner EsriShapeFileCutter EsriShapeFileFilter Y gt m Filter Action Insert Before Insert After Delete Edit Properties close didi Image Chain 3 ccf Band Selector Histogram Remapper Tile Cache Image Renderer Tile Cache Figure 22 Chain editor Select EdgeFilter left and select the image handler ccf on the right Top is the input source followed by the rest of the processing filters and then do an insert after m Non Commercial Image Chain Editor Image Layer layer 1 image chain C data ccf_utm foo ccf Available Filters ConvolutionFilterl D ConvolutionSource EastingNorthingCutter m ElevimaqeSource EquationCombiner EsriShapeFileCutter EsriShapeFileFilter X gt m Filter Action Insert Before Insert After Delete Edit Properties close didi L Image Ch amp in 23 ccf ossimEdgeFilter Band Selecto
51. insertion into the chain Bring up the property editor by double clicking on the filter Change the Edge type option to see different edge detectors Initially it holds the value of Sobel Select the field and it will have a drop down combo box for selecting other edge detection options 8 Non Commercial Property List 2 x Class name lossimEdgeFilter H Description enabled True Edge type Sobel M Auto apply Reset Close Figure 95 EdgeFilter properties Page 73 of 107 7 4 Utilities Menu Utilities Window Help Elevation manager Unit converter Figure 96 Utilities menu for editing the elevation manager and unit converter Utilities menu gives quick access to some utility dialogs Currently we have utilities for manipulating the elevation manager and doing unit conversions 7 4 1 Elevation Manager The elevation manager can be accessed through the Utilities gt Elevation manager menu option By default any elevation directories registered with the System is auto loaded Any current elevation cell will appear in the list 8 Non Commercial Elevation Manat x Elevation cels 34 Add Add directory Remove Remove all Move up Move down Move to top Move to bottom ERR Mosaic v elevation enabled MV autoload V auto sorting Close Figure 97 Elevation manager dialog interface Add button will allow the user to add individual elevati
52. invokes a layer chooser dialog Select the layer that you wish to place on top of the current display window and press the apply button E Non Commercial Layer Chooser 32 Image chain C data test data sanfran sanfran ccf 57 lmage chain C data test data sanfran map highres utm map tif Cancel Figure 47 Layer chooser dialog box Select the second layer for swiping Page 45 of 107 The resulting swipe window is finished aA Non Commercial unnamedo prj Ele Layer Vce Utiiies Window Help Image chain C data test_data sanfran sanfran ccf Ele Edit Swipe Enhancements Help Zoom Pan Full res Propagate Figure 48 Result of adding a layer for swiping 7 3 2 Remove Layer Once an additional layer has been added it can be removed using this command Simply select the image chain to be removed and click the Apply button 7 3 3 Horizontal Swipe modes are ideal for change detection or rapidly comparing two geographic coincident data sets As the cursor is moved back and forth the top layer covers and uncovers the base layer The horizontal swipe is the default swipe when adding a new layer If you are not in horizontal swiping it can be selected by the Swipe gt Horizontal menu option 7 3 4 Vertical The vertical mode is selected through the Swipe Vertical menu option This give another direction for detecting differences Page 46 of 107 7 4 Enhancements Menu E Non Commercial
53. l lt directory where plugins are gt you can also list by individual file names plugin filel lt full path and file name gt Some keywords for the OSSIM GUI application igen spec files output This is the directory location you would like the igen export GUI to default to when outputing spec files ossim igen spec output directory work gpotts igen test this is the igen executable Make sure you put the full path ossim igen executable work gpotts ossim bin igen END KEywords for OSSIM GUI 5 1 1 Geoid Setup Current version supports the 96 99 and 2003 grids The 99 and 2003 grids are distribution from NGS and the 96 grid was a distribution from Geotrans Assuming you have downloaded the test data and have it saved to c ossim then you should have a geoid96 geoid99 and a geoid03 under the c ossim test_data directory For 96 geoid support use geoid egm 96 grid c NossimNtest dataNgeoid96Negm96 grd The 1999 grids and 2003 grids are NGS distributions and use the keyword Geoid ngs directory c NossimNtest dataNgeoid99 geoid ngs directory byte order big endian Or Geoid ngs directory c NossimNtest dataNgeoid03 geoid ngs directory byte order big endian Page 22 of 107 In future releases we will support multiple geoid directories and searching for correct geoids to use For now you can only use one or the other 5 1 2 OSSIM Planeta
54. l input data was say UTM zone 10 5 meters then they could be put together with an Ortho combiner and remove the resampler out of each chain and put it after the combiner This saves on a lot of processing Page 55 of 107 Max combiner takes all input layers and finds the max pixel for the output of each pixel value Closest to center combiner will replace the output pixels with the closes valid pixel to the center of the image See Figure 105 for closest to center example 7 5 4 Fusion Manager Import image file Histogram match Combine Color normalized Correction id Local Correlation Elevation SFIM Figure 65 Feather combine left simple combine middle and a blend combine The fusions implemented in OSSIM take a low resolution multi band image and a high resolution grey pan band as input to the fusion process Each algorithm tries to enhance the lower resolution image with the higher resolution details of the pan image NOTE The color inputs and the pan inputs could have been from a mosaic layer and you are not restricted to single image layers Typical fusion may be Landsat Multi spectral with the co registered pan band Color Normalized Fusion is a slightly modified Brovey fusion It is important to note that this fusion only works on 3 band inputs Local Correlation uses local linear regression on a 5x5 window to correlate the pan to the same area on the multi spectral data The fusion is only applied a
55. l raster bip envi general raster bil envi general raster bsq envi nitf block band separate nitf block band sequential If the GDAL plugin is supported it will depend on the GDAL build on what output formats are enabled For my machine I have gdal VRT gdal GTiff gdal NITF gdal HFA gdal ELAS gdal AAIGrid gdal DTED gdal PNG gdal JPEG gdal MEM gdal GIF gdal XPM gdal BMP Page 84 of 107 gdal PCIDSK gdal PCRaster gdal_ILWIS gdal_PNM gdal_ENVI gdal_EHdr gdal_PAux gdal_MFF gdal_MFF2 gdal_BT gdal_IDA gdal_FIT gdal_RMF gdal USGSDEM The plugin output types are all prefixed with the gdal 10 2 orthoigen orthoigen application is a powerful product generation tool Orthoigen generates an igen spec internally and calls the same code the igen application uses Command line options are explained followed by before and after screen shot samples Usage orthoigen options input files output file Note The last file on the command line is assumed to be the output filename The output writer type is determined by that extension unless a writer template is given orthoigen is also an MPI aware applications If MPI support was included during the build stage then you can use MPI for the processing For example If you have LAM MPI and OSSIM was compiled with MPI turned on then you could execute mpirun np number processors orthoigen options lt input files output file and would use number processors
56. long correlated areas This fusion will probably not work well in cross sensor fusions This work was derived from the paper SFIM or Smoothing Filter based Intesity Modulation The Local Correlation and the SFIM fusion are band separate and will fuse all bands whereas Color Normalized expects the input to be 3 bands only When the scale resolution is more than 2x you will probably want to adjust the fusion parameters through the display window Enhancements gt fusion menu option For an example in the directory doqq in the test data suite are two images doqq5 meter tif and doqq_1meter_grey tif This simulates a 5x difference in resolution and we would like to enhance the doqq_5meter color image to the 1 meter grey and create a new high res color product Examples will be shown for each fusion Note the color normalized fusion has no adjustments available but Page 56 of 107 the Local Correlation and SFIM fusion do have parameters to tweak and is a necessity in this particular example First bring both of the images into imagelinker and select the grey scale image and hit the full resolution button You should be at a 1 meter view Now propagate the view to the color window by hitting the Propagate button You will notice a blurred image since the color image was a 5 meter product and is now zoomed to 1 meter Change the resampler type to be something smoother For this example we chose cubic Fig
57. lres Propagate Refresh nea 37949 28 6189 N 122 22 16 9625 W 32 1553 m rgb 1 37 49 29 4760 N 122 22 15 3360 W 32 1547 m rgb 1 4 Figure 7 Full resolution of the Landsat San Francisco scene propagated to all other displays Layer Vce Utlites Window Help With the displays zoomed to treasure island as Manas shown above blend the views In the main menu Hd peat uten select Blend under the Combine submenu under Histogram match the Layer menu Combine Mosaic Fusion Correction Feather Layer gt Combine gt Blend Elevation Merge bands ortho Max P f Closest to center Figure 8 Combine gt Blend menu option Factory A dialog box will appear with a list of image chains that can be combined Select both the map and the image and press the Apply button Page 12 of 107 FA Non commercial Choose Layers to Blend p highres_utm_map tif Cancel Figure 9 Blend layer chooser A new display appears that is blended press the Fit button on the new display to get an overview You will notice that the data sets are blended together in the overlap area 8 Non Commercial unnamedo prj inl xi File Layer Vce Utilities Window Help ossimBlendMosaic ESTES ETE File Edit Swipe Enhancements Help Help Zoom Pan Fulres Propagate Refresh nearestneighbor agate Refresh nea VA
58. mmercial unnamedo prj Ble Layer Voe Uttities Window Help Non Commercial image oix outputimage outputimage type E ift tiled_band_separate x editwriter geometry lines samples n IE use scene rect sevespecfleonly ssveimage close Figure 39 Scale changed to 5 meters Region of interest box changes dimensions Notice the rectangle box changes size to match the scale estimate Now change the dimensions of the box by entering 1024 in the lines field and 1024 in the samples Make sure you hit return in each field or the change will not take affect E Non Commercial unnamedo prj Ble Layer Vce Uttities Window Hep Image chain data test_data sanfran sanfran ect Ele Edt Swipe Enhancements Help Fullres Propagate Refresh nearestneighbor xl Non Commercial image generi BEE outputimage output image type tft tiled_band_separate x editwriter lines geometry samples 1024 1024 use scene rect sevespecfleonly ssveimage close E Figure 40 Lines and samples changed to 1024x1024 Region of interest box changes dimensions You can also move the box with your middle mouse button Move the mouse over the box and press and hold down the middle mouse button and drag the box The b
59. n 3 of this document for instructions on how to get and install the software for your platform For this example we will use the Landsat image and raster map of San Francisco Uncompressing the demodata tgz file on ossim telascience org ossimdata you should find the following sanfran sanfran ccf sanfran map highres utm map tif These files have been prepared as a test set for the quickstart guide to the ImageLinker program ImageLinker will take advantage of auxiliary files if they are present There are several OSSIM command utilities to create these files for data setsThe initial files were sanfran tif and highres utm map tif both geotiff files GeoTiffs are tiffs that have geospatial tags embedded in the file The respective ovr files are overview files containing reduced resolutions sets These files can be generated with the img2rr image to reduced resolution command line utility The geom files contain geometry information and were created with the create geom command Histogram files can be created with create histo or can be created when initially ingested into ImageLinker The readme txt files are human readable files about the data sets These two files will be used to demonstrate some of the basics of ImageLinker So start up ImageLinker and follow the following steps in the rest of this section Page 7 of 107 3 1 Open the source files Open the data files for the map and UTM projected color fused Landsat image of Sa
60. n Francisco by navigating the file browser to your test data installation and open them The following command assumes they are located in install location test data File gt Open navigate to test data sanfran san fran tif File gt Open navigate to test data sanfran map highres utm map tif The ImageLinker canvas should reflect that the map and the Landsat image have been loaded 37 56220263 N 122 1525 1861 W 32 016 m rgb 12 10 18 37954 52 3305 N 122926362867 W 32 156 m rgb 0 0 0 Figure 2 Initial opening of the UTM projected San Francisco Landsat left and the San Francisco map right Page 8 of 107 3 2 Managing Windows In this introductory tutorial we will get a feel for some of the window management and display commands Each project in ImageLinker lives in its own canvas Standard maximize minimize and collapse menu options work on the project canvas within the desktop Additionally there are a couple of useful commands for organizing the windows Lets try a couple Select Window gt Cascade The open displays are cascaded within the project canvas A Non Commercial unnamed0 prj n x Ele Layer Vce Utilities Window Help Image chai a test data sanfran s Image chair jata test_data sanfran_Mas Je E Ble Edit Swipe Enhancements Help Zeom Pan Fit Fullres Propagate Refresh nearestr 37954 52 3305 N 122 26 36 3857 W 32 156 m rgb 0
61. n setup Sample Use executable disable elev his disables a specified notification ALL WARN NOTICE INFO FATAL DEBUG Sample Use executable disable notify ALL his option will disable loading of all plugins Specify an individual keyword name value pair to be passed land added to the preferences It takes one string argument Sample Use executable K name value Specify the classes to trace Takes one argument that is a quoted string in a regular expression format The string indicates which classes to output debug messages for Sample Use executable T ossim Will trace all classes that start with the word ossim executable T ossimElev ossimImage Will trace all classes that start with ossimElev and ossimImage akes a log file as an argument All output messages will be lappended to the given log file name Sample Use executable ossim logfile mypath logfile txt Important to note that disable notify and K can appear more than once on the command line For example executable disable notify WARN disable notify NOTICE will disable WARN and NOTICE messages For K option we can have K namei value1 K Page 81 of 107 name2 value2 etc for all name value pairs you wish to update add to the global preferences 10 1 icp Image copy or icp can translate one image format into another If the image is map projected it will try to carry over t
62. ne implementation of the MPI standard The location of LAM MPI s main web site at the time of this document is at www lam mpi org This is a robust MPI implementation and is currently tested under OSSIM v1 6 4 within the UNIX flavor environments NOTE The Windows version does not come with LAM MPI support pre built into the OSSIM libraries Orthoigen igen and pixelflip have LAM MPI support and have achieved near linear speedups The speedups become apparent when the processing time increases and requires a process flow that is more complex than just an equivalent image copy Adding reprojections and other area analysis filters would benefit from a LAM MPI environment Future releases of OSSIM will be migrated to support LAM s new implementation called OpenMPI This should be no more than a recompile of the OSSIM core libraries to the new OpenMP libraries 10 Command line Applications All command line applications have the following common options D ee h or help Should display some help for the command line application It will display help for the common options and any additional options required by the command line application being ran Page 80 of 107 disable notify his disables elevation lookups It takes no arguments If you are manipulating already map projected data it is more efficient to disable elevation overhead This can give you a substantial increase in performance depending on how you have your elevatio
63. ns Ea Non Commercial histogram x stretch mode linear_one_piece 0 0000 low clip percent display channel master 0 0000 high clip percent Zz low clip value 255 high clip value fo mid point iim output min value 255 output max value histogram file C date ccf utm foo his set histogram file enabled apply reset clase Figure 52 Histogram operation editor ImageLinker contains a dialog for viewing and manipulating image histograms The current version of ImageLinks needs to attach to an external histogram file If one is not present it will ask if you would like to create one for the input file An alternative is to use the command line application create_histo The dialog allows the setting and application of clip points and has selectable standard deviation values from the stretch mode list Page 49 of 107 7 4 5 Hue Saturation Intensity Adjustments hue offset TS B 0 hue low range E ae J 0 hue high range _ _ _ n hue blend range 0 saturation offset lt a 0 00 intensity offset 7 _ _ 000 low intensity clip 7 000 bihihtensiy Clm Sj 4 X8 MW 2 1 00 white object clip 1 00 T a G enable reset all close Figure 53 Hue Saturation and Intensity editor This dialog provides color adjustments via hue saturation and intensi
64. o center combiner left A feather combiner right There are several other combiner types not shown ossimBlendMosaic Blends all layers into a single average value Areas that have only one image are just copied ossimMaxMosaic Takes the maximum value of all layers ossimColorNormalizedFusion Takes 2 input images where one is a color image and the other is a pan grey scale image Creates a pan sharpened product ossimLocalCorrelationFusion Takes 2 inputs where one is a color image and the other is a pan grey scale image It s an adaptive fusion that uses the local correlation between the two images to sharpen the color image Note This currently should only be used in the Imagelinker GUI since do not automatically determine the scale difference between the pan grey scale image and color image These parameters will need to be hand tweaked in the editor ossimSFIMFusion Acronym stands for Smoothing Filter based Intesity Modulation and takes 2 inputs where one is a color image and the other is a pan grey scale image Note This currently should only be used in the Imagelinker GUI since do not automatically determine the scale difference between the pan grey scale image and color image These parameters will need to be hand tweaked in the editor This algorithm appears to work the best at preserving spectral information ossimBandMergeSource Takes all input files and merges the bands together and produces a merged output So if you
65. oid96 directory holding the geoid 96 grid distributed by geotrans katrina directory holding some sample aerial data for hurricane Katrina Example of the Applanix sensor model orthoigen directory holds the sample keywordlists and tile data files to run the samples in the orthoigen section 7 2 sanfran directory holds a 5 meter color product originally created from a color Landsat and Spot fusion sanfran map directory that contains a UTM map of the San Francisco area 4 3 Source Code for Developers The main site for accessing the source download link is at www ossim org Selecting the downloads will take you to the sourceforge project tree The released source code is available for download and also contains dependency source code For CVS access we host and own the following modules ossim ossim qt ossim plugins libwms osgPlanet osgplanet qt Page 18 of 107 4 4 Preferences Setup OSSIM preference file are global preferences for all command line a GUI applications Preferences are enabled either by passing the file as command line using the P option or bey specifying it as an environment variable OSSIM PREFS FILE Assume you have a preference file stored at home myaccount ossim preferences then you can set your environment variable OSSIM PREFS FILEz home myaccount ossim preferences or pass as command line options P home myaccount ossim preferences 4 4 1 Elevation Setup OSSIM supports loading of SRTM and DTED elevati
66. on cells Currently only SRTM and DTED formatted cells are loadable Page 74 of 107 Add director button will find all cells in the specified directory and add them to the elevation manager Remove button will remove the currently selected cell from the elevation manager Remove all button will remove all elevation cells Move up button will move the currently selected cell up one location Move down button will move the currently selected cell down on location Move to top button will move the currently selected cell to the top of the list Move to bottom button will move the selected cell to the bottom of the list Mosaic button will mosaic all selected cells together and add it to the data manager Note no overviews are created in this process If your cells don t have overviews and you want to be able to zoom efficiently then you will need to create overviews externally from ImageLinker with the command line application call img2rr elevation enabled checkbox enables disables the elevation lookup values in the elevation manager auto load checkbox will enable disable automatic loading of cells If a height for a lat lon position is not found in the current cells then it will find the cell in the manager and load it If the autoloading is disabled i e not checked then it will not automatically try to find the cell to load auto sorting button will sort the loaded cells based on accuracy 7 4 2 Unit Converter The Unit Conv
67. on data Let s assume you placed the test data under c ossimNest data elevation sourced1 filename c NossimNtest dataNelevationNsrtmN1arc elevation sourcedl type srtm directory Note elevation source1 type is optional and if not present it will use auto detection to determine the type The other supported type is dted directory srtm cell and dted cell For an up to date keyword list template please look in source distribution ossim etc templates ossim preferences template for a complete list of supported detectors If you have other source directories for elevation you just increment the number and add another one elevation source2 filename elevation source n filename You can mix SRTM and DTED and OSSIM will handle it Note Elevation lookups will be in the order you specified So if you have SRTM 30 meter and DTED 30 meter and you have a cell that overlaps a lat lon value then it will use the first one in the list and if the height is NULL Invalid then it will go to the next manager in the list until a valid height is found Currently they are assumed to be relative to a geoid grid Mean Sea Level Other keywords that can be controlled are elevation enabled true elevation auto load dted enabled true elevation auto_sort enabled true elevation compute_statistics enabled true Page 19 of 107 5 if
68. on information 2 Non Commercial Image Information Image Layer IJ xl zi Image Geometry Image Projection Valid vertices Projection ossimUtmProjection z Datum WGE World Geodetic System 1984 z Units Meters z Pixel type Center of pixel paint z Zone 10 z Hemisphere Northern z Origin Latitude 00 00 00 00000000 N Central Meridian 123 00 00 00000000 W False 500000 000000000000000 False Northing 1000000000000000 Tie Point 538185 000000000000000 Tie Point Y 4208785 000000000000000 Scale X 5 000000000000000 Scale Y 5 000000000000000 Scale Factor STD Parallel 1 00 D0 00 00000000 N 00 00 0000000000 N Reset Close Page 34 of 107 Figure 26 Image information window Image Projection tab Image Projection information is currently only given for map projected data Finally valid vertices are displayed by selecting the Valid Vertices tab Image Geometry Image Projection Valid Vertices g Valid Vertices Vertex 0 Vertex 1 Sample Line Sample m Line Vertex 2 Vertex 3 Line 1 Sample Line Sample Scan for valid vertices Help Close Figure 27 Image information window Valid Vertices tab The Valid Vertices tab displays the images valid vertices Many time such as Landsat scenes the images
69. on requires 4 points 2 This projection has a custom dialog that requires you to drop geogrphic points Would you like to drop points Figure 33 Question box when changing to a bilinear projection Select Yes and the bilinear projection drop point dialog will come up Page 38 of 107 FA Non Commercial quad projection dialog 7 3 x Select row use left button on image to drop point Edit cells as needed C roaming enabled Figure 34 Drop point GUI Selecting in the window will add an adjustment to the dialog Select the in the corner of the cross section where 41 degrees North and 117 degrees west meet x 2 Adjust point for row 1 wj Figure 35 Popup that comes up when selecting with the left mouse button in the image window Selecting yes will add the line and sample to the selected row Type in the lat lon location EA Non Commercial quad projection dialog j 2 x Select row use left button on image to drop point Edit cells as needed 275 41 542 16 41 00 00 0000 N 117 00 00 0000 W 0 000 228 58 8821 54 41 00 00 0000 N 115 00 00 0000 Ww 5689 79 888466 40 30 00 0000 N 116 00 00 0000 W 0 000 a 573866 5442040 30 00 0000 N 117 QU 00 0000 W 17 167 Datum WGE World Geodetic System 1984 E roaming enabled delete row save reset close Figure 36 Drop point GUI Page 39 of 107 Repeat f
70. or all 4 corners of the image and then hit save the geometry to test data drg battlemn geom You now have a projection in OSSIM and you can output the product to a geotiff 7 2 Image Generation Image chain C Ele Edit Swipe Enhar Save AS Print Ctrl P Figure 37 File menu under the Image display The File gt Save As allows the user to set parameters for product output save those specifications to an external spec file and or save the image to disk When the File gt Save As menu option is hit the product output dialog box comes up 2 Non Commercial unnamed prj 3 lol x Ble Layer vce Utlties Window Hep Image chain C data test data sanfran sanframecf Ele Edit Swipe Enhancements Hep Zoom Fen Fulres Propagate Retres neerestneighbor x 2 Non Commercial image oulputimage output image type tittled_band_separate v editwriter geometry scele lines samples Bo49773755656077 aat n use scene rect save specfile only saveimage close Figure 38 Image generation dialog Once this dialog is up it will link to the current display Initially you should se a white bounding rect covering the entire scene There are several was to change the dimensions of the rectangle or redoing of interest Try changing the scale field in the Image generator dialog box to 5 Page 40 of 107 2 Non Co
71. ox should change to a green color Once you have the box moved to the desired location release the middle mouse button In this example we moved the box over Treasure Island Page 41 of 107 2 Non Commercial unnamedoprj Ble Layer Vce Uttities Window Help Ele Edit Swipe Enhancements Hep Zoom Fen Fulres Propagate Refresh nearestneighbor P Non Commercial image generi oix outputimage output image type ift tiled_band_separate v editwriter geometry scale lines samples s000000000000000 1024 1024 use scene rect sevespecfleonly saveimage close Figure 41 Box moved over Treasure Island in the San Francisco image Using the left mouse button can create a custom box Press and hold the left button down while dragging A box should be created When you have reached the desired size release the mouse button E Non Commercial unnamedo prj Ele Layer Vce Utiities Window Help Image chain data test_data sanfran sanfran cch Ele Edit Swipe Enhancements Hep P Non Commercial image Zoom Fan Fullres Propagate Retresn nearestneighbar x lolx outputimage output image type Eil i ied_band_separcte S geometry scale lines samples s000000000000000 4573 fasts use scene
72. present themselves in the same view as the window where the button was pressed In this example the full resolution view of the Landsat image on the left was propagated to the map on the right Both images are presented in the same scale and area of coverage 7 View Accessing Edit gt View will bring up the view dialog box FA Non Commercial View Geometry 2 x Transformation C State Plane Datum WGE World Geodetic System 1984 x Projection ossimUtmProjection x Meters per pixel 5 000000000000000 000000000000000 OriginLetion 0 000000000000000 123 000000000000000 Zone INMEEEEN Hemisphere Noth Parallels 1 amp 2 p t Scale factor e 1 Help Propagate Apply Close A Figure 19 View dialog box Page 28 of 107 This dialog displays and modifies the view projection of the displayed image Currently only map products are supported Transforming to another projection datum is as simple as selecting the desired output from the pull down menus Parameters for resolution origin etc are also available from this dialog The OSSIM application is projection and resolution independent All transformations and resamplings are performed automatically to meet the needs of any display To apply changes you can either hit the return key or click on the apply button The propagate button takes the same action as the propagate button on the displa
73. r Histogram Remapper Tile Cache Image Renderer Tile Cache Figure 23 Chain editor EdgeFilter inserted after the selected ccf image handler Page 32 of 107 and finally the result data ccf_utm foo ccf Bc nearest neighbor v File Edit Swipe Enhancements Help Zoom Pen Fit Full res Propagate Refresh Figure 24 Result of inserting the EdgeFilter 7 1 3 Image Info Selecting the edit gt image info menu option will bring up the image information window for the current image in the display Page 33 of 107 Non Commercial Image Information 2x Image Layer Image Geometry Image Projection Valid vertices Image Type Radiometry Byte Order cct ossim uint amp Lines Samples Bands R Levels 8000 goo g e Selected Band min mex null values Min 1 0000 Bend n s Max Null 255 0000 0 0000 Scan For Min Max Help Close Figure 25 Image information window Brought up through the edit gt image info menu option on the Image display window A This dialog contains geometry projection and vertex information on the displayed image The scan for min max is currently disabled and can only be ran by the external application cmm which stands for compute min max Each section is available through the tabbed interfaces on the dialog The Image Projection tab shows the images input projecti
74. r in the top chain and drag the mouse over the top left red triangle on the display icon A display window should appear in the Main window and both arrows should go to green showing that they are linked Page 69 of 107 Figure 89 Display connected to the chain Double clicking on any filter in the chain will either bring up a property list editor for each editable property that has been exposed by OSSIM or a custom dialog editor The property list editor is a fallback editor if no custom editor has been written Double click on the left most filter of the top chain m Non Commercial Property List Class name ossimCcfTileSource Description enabled v gt entry 0 Auto apply Apply Reset Close Figure 90 Property editor for the image handler Page 70 of 107 Try toggling the enabled button off and on When it is disabled data does not flow and a NULL empty image is sent through the chain indicated by the Black output 2 Non Commercial unnamed0 prj Ele Layer vce Utilities Window Help Visual Chain Editor Bee Ble Edit amp 1 3 I9 al 4 m gt Sm AE Property Value Class name imCcfTileSource Description enabled F zj eny 0 iL P F Auto apply Reset Clase Figure 91 Toggling of the Image Handler s enabled property Filters can be added to the chain by drag and dropping a
75. r of machines for faster product generation Note I have not done windows builds of mpi Page 107 of 107
76. rect 37 43 39 9998 N 122 13 45 6546 W 32 0268 m rgb B1 79 66 save spectile only saveimage clase Zi Figure 42 Custom bounding box drawn using the left mouse button The use scene rect button located in the Image Generator dialog box will stretch the region of interest to cover the entire scene The scale is in meters and corresponds to the meters per pixel product output The output Image field has a file browse button indicated by the three dots When pressed it will open up a file browser to identify the location of your product If you edit this field manually without the file browser you must hit return in order for the change to take affect Also note we currently do not automatically set the correct extension based on Page 42 of 107 writer type and you must specify it The output writer type will have a list of output types to choose E Non Commercial image generator output image output image type x ift tiled bend seperate x edit writer tiff tled band separate 4 Jpeg Jeome general raster bip scale lines general raster bil i0 497737556561077 441 general_raster_bsq se scene rect general raster bip envi general raster bil envi save spec file only _savg general_raster_bsq_envi nit block band separate Figure 43 Image generation dialog writer type selection Once the writer type is selected individu
77. running program This building block approach allows complex image processing flows to be interactively constructed and modified Each object or image unit in the chain may have its own controls and adjustable parameters The entire state of the chain including adjusted parameters can be easily saved and retrieved for OSSIM enabled programs Page 77 of 107 9 2 File formats supported mu E E E WrwaBeww p y E Wm 77 E E E WSewmexme 0p o oo E Wemes 0p E E pe EC EE GeneralRaster Band Interleave by pixel BIP Band Interleave by line BIL bem p E E mu 770 E RE Enemy t E pew ww 0 FR E bas OR m0 EE E b E RE In the table P represents the GDAL Plugin For file formats that don t support geo encoding an external geometry is supported Jpeg format and the general raster formats are examples of a supported external georeference GDAL has a host of other exporting capabilities and a complete list can be found at www gdal org The ones listed in the table are the only ones we have tested Page 78 of 107 9 3 Projections and Transformations Supported OSSIM supports RPC Landsat7 and Landsat5 Spot5 Applanix ECEF and UTM orientations and Coarse grid models and a host of map projections and state plane codes and datum shifting transformations Much of the Map projections and datum shifting code was brought over from Geotrans Geo
78. ry Viewer Setup osgPlanet The windows binary installer will automatically setup the preference file with a default low resolution background reference image This will be used as the base texture map for the earth For other platforms the keywords may have to be added manually to your preference file ossimosgplanet background type local ossimosgplanet background file0 C Program Files Common Files ossim images earth tif ossimosgplanet background transparent color flag 0 ossimosgplanet background transparent color 0 0 0 ossimosgplanet background opacity 255 The type currently should be local and the files can be listed with first keyword starting at fileO and going to file lt n gt You can have any number of images as a background reference and can be any image loadable by the OSSIM core library 5 1 3 Plugin setup At the time of this document a Geospatial Data Access Library GDAL was the only supported plugin t gives access to a host of other data formats Under windows this is automatically setup for you Here is a sample of the keywords for unix plugin file1 Users Shared Development ossim plugins lib libossimgdal plugin dylib If you are working from the latest snapshot there is a registration plugin started plugin file2 Users Shared Development ossim plugins lib libossimreg plugin dylib Page 23 of 107 6 ImageLinker ImageLinker is a GUI based application demonstrating the power of the OSSIM open source
79. ssssssessosesessssecssossesssssssssse sesso esee igen Page 4 of 107 2 About OSSIM This document provides an overview of the Open Source Software Image Map OSSIM project http www ossim org OSSIM is a contrived acronym pronounced Awesome The core of OSSIM is a C software library that provides advanced remote sensing image processing and geo spatial functionality A quick summary of OSSIM functionality includes ortho rectification precision terrain correction rigorous sensor models very large mosaics and cross sensor fusions a wide range of map projections and datums and a large range of commercial and government data formats The architecture of the library supports parallel processing with mpi a dynamic plugin architecture and dynamically connectable objects allowing rapid prototyping of custom image processing chains Around the core ossim library the software distribution includes a large number of command line utilities that can be easily scripted for batch production systems and higher level GUI applications imagelinker and iview Additionally bindings have been generated for other languages Also included in the distribution are the osgPlanet and osgPlanet qt modules These modules build on top of OSSIM and OpenSceneGraph to provide geo spatially accurate 3D visualization capabilities On 4 Feb 2006 the Open Source Geospatial Foundation http www osgeo org was formed in Chicago IL The OSS
80. ta test data elevation srtm 1arc N37wW123 hgt 380 Image Ortho chain Ortho chain data test_data elevation srtm 1 arc N38W123 hgt sanfran map highres utm map tif Figure 75 Choose layers for hill shading Page 61 of 107 Resulting in the image ossimBumpShadeTileSource Propagate Refresh 37950015504 N 12230265504 W 32 294 m rgb 0 0 0 Figure 76 Hillshaded image Page 62 of 107 The hillshade can be edited through the Edit gt Layers menu option 8 Non Commercial Combiner Editor 2 xl Available Inputs 329 Image chain C data test de 364 Image chain C data test de 377 Image Ortho chain Ortho che 390 Image Ortho chain Ortho che 7E 455 ossimBumpShadeTileSourc gt Up Down Elevation angle 45 000000000000000 Azimuth angle 45 000000000000000 Smoothness factor 1 000000 409 ossimOrtholmageMosaic 436 Image chain C data test dc Iv Auto apply Apply Beset Close Figure 77 Adjustable parameters for the Hill Shading algorithm Shading of the image can be adjusted by changing the sun s elevation and azimuth angles Azimuth s O starts pointing in the direction up the screen so an azimuth of 0 says the sun is at the top pointing down In this example we have the azimuth and elevation set to 45 degrees which states the sunc is coming from the direction top right and is elevated by 45 degrees above the plane The smoothness factor is used to change
81. te Open an image and then in the file save as set your writer definitions and then output all your settings to the igen spec file Edit the spec file and copy out the writer text located with object2 prefix value Move this to a file called writer template kwl Although the current orthoigen will take it as valid you really should strip the object2 prefix from the writer keywordlist For example I get the following from when extracting the writer from the igen spec file created from the imagelinker application object2 color lut flag O0 object2 compression quality 75 object2 compression type none object2 create envi hdr O0 object2 create external geometry 0 object2 create fgdc O0 object2 create histogram 1 Page 98 of 107 object2 create image 1 object2 create jpeg world file 0 object2 create overview 1 object2 create readme 0 object2 create tiff world file 0 object2 description object2 enabled 1 object2 filename data ear1 output docs writer tif object2 id 33 object2 image type tiff tiled band separate object2 input_connection1 9 object2 input_list_fixed 1 object2 number_inputs 1 object2 number outputs 0 object2 output geotiff flag 1 object2 output list fixed 0 object2 output tile size x 64 object2 output tile size y 64 object2 overview compression quality 75 object2 overview compression type 1 object2 pixel type point object2 projection units meters object2 type ossimTliffWriter Now strip th
82. ters must be specified in an OSSIM keywordlist The Applanix generates a calibration report in the form of a pdf Copy those parameters to the associated keywords A sample OSSIM Applanix camera file is provided sensor sn0016 focal length 54 909 principal point 032 331 Page 105 of 107 ima ge_size 4077 0 4092 0 pixel_size 0 009 do dl d2 d3 d4 d5 d6 d7 d8 d9 d10 dll d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 1 0 02 2 0 2 3 0 66 4 1 57 5 3 05 6 5 25 7 8 3 8 12 32 9 17 44 10 23 76 11 31 39 12 40 42 13 50 93 14 62 99 15 76 64 16 91 94 17 108 89 18 127 52 19 147 8 20 169 71 21 193 19 22 218 16 23 244 53 24 272 17 25 300 94 26 330 64 distortion_units microns where sensor specifies the sensor id focal length specifies the focal length in millimeters principal point specified in millimeters image size specified in pixels Typically 4077x4092 The parameter is width followed by height separated by a space pixel size pixel size on the camera in millimeters d0 d25 are distortion parameters Is a tuple that is the distance from the center followed by the distortion distortion units The units of the distortion Usually they are given in microns They will be converted internally to millimeters Examples l applanix2ogeom camera file exterior orientation Will generate all geometries
83. the active project and asks if you would like to save the project file 6 1 4 New Project Creates a new project and associated project window 6 1 5 Save Project Saves the current project and invokes a dialog box if the project hasn t been assigned a name 6 1 6 Save Project As Invokes a Save dialog box and allows the user to assign a new name to the project 6 2 Display Menus When the File 2Open Image menu is invoked the image comes up in a display window Several menu options are available for advanced image manipulation Image chain C data test data sanfran sanframccf SE E File Edit Swipe Enhancements Help Zoom Pan Fit Full res Propagate Refresh cubic gt Figure 16 Image display menu options and toolbar buttons Page 26 of 107 On the top row we have the basic file edit swipe enhancements and help menu options On the second row we have a toolbar with quick access to zooming panning fitting full resolution propagate refresh and specification of the resampling type Each are explained in further details 6 2 1 Zoom and Pan Modes The window palette contains two icons for selecting zooming and panning modes for the display Select the mode then control by clicking with the left mouse button in the content of a display When in the zoom mode clicking in the window will do a power of 2 zoom in Shift click will cause a power of 2 zoom out Dragging a rectangular area will zoom the display to th
84. the file extension Note this didn t make it into 1 6 4 release but is available from current Cvs t or thumbnail Will manipulate the view so the output product corresponds to the thumbnail dimension for example thumbnail 512 will produce a 512x512 close to it thumbnail image of the product being generated Page 87 of 107 For all examples we will use the following overlapping input images Figure 100 This will correspond to filename tile1 jpg and is a 256x256 UTM input file Page 88 of 107 Figure 101 This will correspond to filename tile2 jpg and is a 256x256 UTM image This has been brightened to show visual differences in the algorithms One of the images has been brightened to show visual differences of certain operations These images are also located under the test data orthoigen directory on the test data distribution 10 2 1Reprojection A view must be specified for the orthoigen application to execute properly First operation will reproject the UTM tile1 jpg image to a geographic output tif orthoigen geo tilel jpg output tif Page 89 of 107 Figure 102 Geographic projected output Tile1 jpg is a UTM 5 meter product and is reprojected to a 5 meter geographic projection See the difference between Figure 100 Two other quick options are utm and input proj The utm will reproject to utm and the input proj will default the output view to the input images projections For more complex
85. the last argument Make sure that the name generated is unique Supplying only a column id or just row id will generate same named files but putting both a row and column will produce a unique name The above command line produces a 4x4 grid of tiles with a total of sixteen tiles named pixel O O jpg pixel O 1 jpg pixel O 2 jpg etc Only the first row is shown Page 96 of 107 Figure 108 These are the fist row of tiles produced from the orthoigen pixel tiling example The mask could also have been pixel_ i jpg and would have generated pixel O jpg pixel 1 jpg pixel 15 jpg Let s specify the tiling in meters and put it in a file called it meter tiling kwl will create a tiling for the 5 meter UTM file in Figure 100 igen tiling type ossimTiling igen tiling tiling distance 320 320 igen tiling tiling distance type meters igen tiling delta 5 5 igen tiling delta type per pixel igen tiling padding size in pixels 0 0 This will generate a tile that is 320 meters wide by 320 meters high and will generate the tile based on a 5 meters per pixel setting This should give something along the lines of a 64x64 pixel tile where the gsd is 5 meters Now execute orthoigen utm tiling template meter tiling kwl input tile1 jpg meter 96r96 96c96 jpg Currently the tiles are snapped to an even tiling distance You might get some tiles that have a sliver of valid information while the rest of it is black or null Here is t
86. the roughness of the scene The smaller the value the rougher the scene and the larger the value the smoother it will look 7 2 Vce Menu Vce Utlites Windo New chain Components 7 3 Figure 78 Utilities menu for editing the elevation manager and unit converter Page 63 of 107 The Visual Chain Editor VCE is a start of allowing one to build your own chain flows The New chain menu option will bring up a chain editor window The Components menu option displays some of the OSSIM image processing factory contents to allow for visual drag and drop of processing algorithms onto the Canvas window brought up by the New chain menu option 7 3 1 VCE Components The Vce gt Components menu option will bring up a dialog box that allows for drag and dropping of any processing filter onto the Chain editor canvas The initial dialog comes up with the loaders tab selected 8 Non Commercial Vce Components 2 x Combiners Filters Exporters Displays ImageHandler Figure 79 Loader for the VCE chain editor Each tab has a list of objects that can be dragged onto the canvas Each tab has a list of objects that can be dragged onto the canvas window Note currently the exporters tab should not be used since the VCE does not support execution of the exporters Instead of using the exporters tab selection to output products we should use the Displays tab to bring up a display and then do a File gt Save as S
87. trans is a NGA certified geospatial coordinate transformation library Here is a listing of all sensor based ossimProjection classes ossimApplanixEcefModel ossimApplanixUtmModel ossimCoarseGridModel ossimRpcModel ossimQuickbirdRpcModel ossimNitfRpcModel ossimLandSatModel ossimNi tfMapModel ossimFcsiModel ossimRpcProjection ossimSpot5Model Note the ossimRpcProjection is a simplified version of the ossimRpcModel Here is a listing of all Map projection based models ossimAlbersProjection ossimAzimEquDistProjection ossimBonneProjection ossimBngProjection ossimCassiniProjection ossimCylEquAreaProjection ossimEquDistCylProjection ossimEckert4Projection ossimEckert6Projection ossimGnomoni cProjection ossimLambertConformalConicProjection ossimL1LxyProjection ossimMercatorProjection ossimMillerProjection ossimMollweidProjection ossimNewZealandMapGridProjection ossimObliqueMercatorProjection ossimOrthoGraphicProjection ossimPolarStereoProjection Page 79 of 107 ossimPolyconicProjection ossimSinusoidalProjection ossimStereographi cProjection ossimTransCylEquAreaProjection ossimTransMercatorProjection ossimUpsProjection ossimUtmProjection ossimVanDerGrintenProjection ossimSpaceObliqueMercatorProjection Here is a listing of some miscellaneous projections ossimQuadProjection ossimWarpProjection ossimAffineProjection ossimBilinearProjection 9 4 LAM MPI Local Area Multicomputer Message Passing Interface or LAM MPI is o
88. ty sliders The changes can be applied to individual or multiple channels by selecting the radio buttons on the top Clip points and white object clip points can also be directly adjusted The enable checkbox makes the changes active on the parent display 7 4 6 Propagate Resampler This command propagates the current windows resampler selection to all other display windows in the project For demonstration bring in the San Fracisco image and zoom to full resolution into the Treasure Island location and zoom a few times past full resolution so you get a block image and change the resampler type to cubic 3754953609 N 122022120 9142 W 32 1553 m rg 0 0 0 37548959410 N 122922209911 Vj 32 1553 m rgb 224 221 216 Figure 54 San Francisco color left zoomed to 1 meter Page 50 of 107 Select the resampler type in the drop down list box s Image chain 6 date test_data santra B Edt Swipe Echarcements Hep Ek Ent Swpe Ennancements Leb TERES Popa Rees AE 3754955609 N 122922120 9142 W 32 1553 m rgb 0 0 0 57548959410 N 1229221209911 Vj 321553 m rgb 224 221 216 Figure 55 Cubic resampler selected Now make sure the view is the same on both sides and do propagate view Now lets watch the resampler propagate to all other open displays by hitting the menu option Enhancements gt propagate Resampler e 4 d 37 48 92 75
89. unnamedoO prj Ele Layer Vce Utilities Window Help Image cha ata test data sanfran sd File Edit Swipe Enhancements Help Zoom Pan Band selection Brightness contrast Topographic correction Histograrn operations Hsi adjustments Fusion Propagate resampler Figure 49 Enhancements menu 7 4 1 Band Selection Access the band selector through the Image display menu option Enhancements gt Band selection x C Grey C N Band available bands selected bands 1 1 d 2 3 3 Figure 50 Band Selection Dialog Page 47 of 107 The band selection dialog is used for selecting which source bands are to be mapped to designated output bands Modes exist for grey scale select one input band True select three bands for red green and blue channels and N band select and arbitrary number of bands The bands can be reordered as well 7 4 2 Brightness Contrast This brings up a brightness contrast GUI window A brightness slider and a contrast slider is present to manipulate the current images brightness contrast levels brightness contrast IV enable Figure 51 Brightness contrast editor 7 4 3 Correction gt Topographic Used for providing topographic illumination correction to imagery based on the normals generated from the underlying elevation surface Used as an error removal process in preparation for material classification Page 48 of 107 7 4 4 Histogram Operatio
90. ure 66 5 meter color image left zoomed to the 1 meter resolution of the grey scale image right Cubic resampler set for the zoomed color image Now select Layer gt Fusion gt SFIM bringing up the layer chooser for the SFIM fusion algorithm A Non Commercial Choose Layers to Fuse Using Local Correlation Ft er_grey tif meter tif Figure 67 SFIM layer chooser Hitting the apply button will initiate the SFIM fusion and bring up a fusion display window Page 57 of 107 I 5 3 ossimSFIMFusion chain Ble Edt Swipe Enhancements Hep Figure 68 Initial SFIM fusion Needs editing for 5x scale change through the Enahancements gt fusion The fusion currently doesn t look that good since the scale changes are not auto detected and assumes 2x by default In this example we have a 5x scale change To make the fusion look better we can tweak individual parameters 88 Non Commercial Parameter Adjustments x Adjustment Selection Adjustment description sigma Prem Peremadjust veue High pass gain 1 1 1 2 Bluning kernel w 7 57142857142857 M Auto apply Keep J Copy Delete Reset Close Figure 69 SFIM fusion Parameters Page 58 of 107 Increasing the blurring kernel width parameter 2 to the value of 5 by editing the value column or by moving the slider to the right gives a better quality fusion This parameter corresponds to the scale change differ
91. view definitions you can use the view template View templates are generally reserved for scripting where other applications drive the execution of the orthoigen application A detailed keyword list for the view template is beyond the scope of this document The easiest way to generate a keywordlist for the view is to use the imagelinker GUI and save out an igen spec file and cut the product view template out of the spec file and save it to a template file that can be used in orthoigen For example opened up tile1 jpg into imagelinker and changed to a California state plane projection using the edit view Next do a file save as and do a spec file only save of the image Open the spec file and cut and paste the product information into a file called view kwl You should get something that looks like the following product projection central meridian 119 000000000000000 product projection datum NAR C product projection ellipse_code RF product projection ellipse_name GRS 80 product projection false easting 2000000 000000000000000 product projection false northing 500000 000000000000000 product projection major axis 6378137 000000000000000 product projection meters per pixel x 5 000000000000000 product projection meters per pixel y 5 000000000000000 Page 90 of 107 product projection minor axis 6356752 314100000075996 product projection origin latitude 35 333333333333336 product projection pcs code 26944 product projection
92. y window It will take the current view setting and propagate to all other displays 7 1 1 Geometry Adjustment Geometry adjustments allow for manual tweaking of mis registered images If the input is coming from a sensor supported by OSSIM you can apply geometry adjustments to the image Supplied under the katrina directory are two Applanix images that are loaded through the Applanix sensor model and can be fine tuned with the geometry adjustment interface B Non Commercial Parameter Adjustments 2 x dusenSebcion E Adjustment description Name Sigma Param Param adjust Value x offset 100 0 y_offset 100 0 orientation x 1 o orientation y 1 0 orientation z 1 a Altitude delta 50 0 slolslolala BERASE IV Auto apply Keep Copy Delete Reset Seve Close j Figure 20 Geometry adjustments for the katrina image The Adjustment selection allows you to switch between different adjustments Each adjustment can have a description by editing the Adjustment description field This dialog is available for sources that are derived with sensor models OSSIM has support for rigorous sensor models and adjustable parameters Each parameter has 5 attributes Name Sigma adjustable param param adjustment slider and the calculated value Not shown is the parameters center The calculation of each parameter is computed as value center param sigma Eac

OSSIM Users Guide - OSGeo Server

Contents

Download Pdf Manuals

Related Search

Related Contents