[GRASS-SVN] r40934 - grass-addons/ossim_grass/r.planet

[svn_grass at osgeo.org]

Author: epifanio
Date: 2010-02-11 16:45:12 -0500 (Thu, 11 Feb 2010)
New Revision: 40934

Modified:
   grass-addons/ossim_grass/r.planet/r.planet.py
Log:
add ossim-orthoigen export ability

Modified: grass-addons/ossim_grass/r.planet/r.planet.py
===================================================================
--- grass-addons/ossim_grass/r.planet/r.planet.py	2010-02-11 19:45:40 UTC (rev 40933)
+++ grass-addons/ossim_grass/r.planet/r.planet.py	2010-02-11 21:45:12 UTC (rev 40934)
@@ -64,9 +64,12 @@
 #% key: r
 #% description: Remove raster
 #%END
+#%flag
+#% key: d
+#% description: Orthoigen 
+#%END
 
 
-
 import sys
 import os
 import socket
@@ -76,6 +79,7 @@
 def main():
 	add = flags['a']
 	remove = flags['r']
+	orthoigen = flags['d']
 	host = options['host']
 	dport = options['dport']
 	pport = options['pport']
@@ -123,8 +127,23 @@
 	if remove :
 		removefile(output,host,dport)
 		print 'Removed raster file :', mappa[0]
+	if orthoigen :
+		path = os.path.dirname(output)
+		print path
+		elevdir = os.path.join(path,'elevation',mappa[0])
+		os.makedirs(elevdir)
+		print elevdir
+		elev = mappa[0]+'.tiff'
+		exportiff(output,elev)
+		instr = make3d(10801, elev, elevdir)
+		os.system(instr)
 
 
+def exportiff(infile,outfile):
+	gdal.GetDriverByName('GTiff').CreateCopy(outfile,gdal.Open(infile))
+
+
+
 def addfile(output,host,dport):
 	ossim_data_xml = "<Add target=':idolbridge'><Image groupType='groundTexture'><filename>%s</filename> <id>%s</id><name>%s</name></Image></Add>" % (output,output,output)
 	ossimdata = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
@@ -231,6 +250,56 @@
 	return center
 
 
+def makedir(path):
+	d = os.path.dirname(path)
+	if not os.path.exists(d):
+		os.makedirs(d)
+
+def make3d(tile, elev, outdir):
+	makedir(outdir)
+	kwl = 'elev.kwl'
+	template = 'igen.slave_tile_buffers: 5 \n'
+	template += 'igen.tiling.type: ossimTiling \n'
+	template += 'igen.tiling.tiling_distance: 1 1 \n'
+	template += 'igen.tiling.tiling_distance_type: degrees \n'
+	template += 'igen.tiling.delta: %s %s \n' % (tile,tile)
+	template += 'igen.tiling.delta_type: total_pixels \n'
+	template += 'igen.tiling.padding_size_in_pixels: 0 0 \n'
+	template += 'object1.description: \n'
+	template += 'object1.enabled:  1 \n'
+	template += 'object1.id:  1 \n'
+	template += 'object1.object1.description: \n'  
+	template += 'object1.object1.enabled:  1 \n'
+	template += 'object1.object1.id:  2 \n'
+	template += 'object1.object1.resampler.magnify_type:  bilinear \n'
+	template += 'object1.object1.resampler.minify_type:  bilinear \n'
+	template += 'object1.object1.type:  ossimImageRenderer \n'
+	template += 'object1.object2.type:  ossimCastTileSourceFilter \n'
+	template += 'object1.object2.scalar_type: ossim_sint16 \n'
+	template += 'object1.type:  ossimImageChain \n'
+	template += 'object2.type: ossimGeneralRasterWriter \n'
+	template += 'object2.byte_order: big_endian \n'
+	template += 'object2.create_overview: false \n'
+	template += 'object2.create_histogram: false \n'
+	template += 'object2.create_external_geometry: false \n'
+	template += 'product.projection.type: ossimEquDistCylProjection \n'
+	open(kwl,'w').write(template)
+	instr = 'export DYLD_FRAMEWORK_PATH=/Users/sasha/OssimBuilds/Release/ ; '
+	instr += '/Users/sasha/OssimBuilds/Release/ossim-orthoigen'
+	instr += ' --tiling-template '
+	instr += kwl
+	instr +=' --view-template '
+	instr += kwl
+	instr +=' --writer-template '
+	instr += kwl
+	instr +=' --chain-template '
+	instr += kwl
+	instr += ' %s ' % elev 
+	instr += '%s' % outdir
+	instr +='/%SRTM%'
+	return instr
+
+
 if __name__ == "__main__":
 	options, flags = grass.parser()
 	sys.exit(main())