merida

examples/downloadTest.py

@@ -5,15 +5,18 @@ sys.path.append('../')
 from geosentinel import APISentinel
 from datetime import date 
 sentinel = APISentinel.APISentinel('asalazarg',  'geo135asg')
-area = "POLYGON((-89.79030210118333 21.122657323983717,-89.77308220413153 21.122657323983717,-89.77308220413153 21.140540053466538,-89.79030210118333 21.140540053466538,-89.79030210118333 21.122657323983717))"
-footprint = "Intersects(POLYGON((-89.79030210118333 21.122657323983717,-89.77308220413153 21.122657323983717,-89.77308220413153 21.140540053466538,-89.79030210118333 21.140540053466538,-89.79030210118333 21.122657323983717)))"
-products = sentinel.getProducts(area, ('20150101', '20180517'), {"platformname":"Sentinel-2", "cloudcoverpercentage":"[0 TO 10]"})
-print len(products)
-for p in products:
-	print products[p]['filename']
-print len(products)
+#area = "POLYGON((-89.79030210118333 21.122657323983717,-89.77308220413153 21.122657323983717,-89.77308220413153 21.140540053466538,-89.79030210118333 21.140540053466538,-89.79030210118333 21.122657323983717))"
+
+area = "POLYGON((-89.5561050716 20.9757315378,-89.5590410642 20.9683665889,-89.5598293382 20.9618947599,-89.5644248796 20.9568453837,-89.5707737727 20.9488129469,-89.5803545363 20.94091954,-89.5912310095 20.9328149479,-89.594380882 20.9236165089,-89.5952350711 20.915838738,-89.5975202078 20.9086333128,-89.5959439943 20.9006865194,-89.5875063006 20.8895372658,-89.5712810763 20.8843713932,-89.5538128258 20.8836099452,-89.5517958006 20.8655171969,-89.5695029128 20.8357088015,-89.5735638178 20.8120013809,-89.5795207653 20.7567549541,-89.5873128094 20.7287324626,-89.6269598778 20.7286922801,-89.6488731594 20.7063652726,-89.6723514384 20.7240677105,-89.6939676021 20.7716497589,-89.6955100309 20.7981175448,-89.6858152774 20.81129115,-89.6958557672 20.8107672715,-89.6955988131 20.8218317231,-89.6759209821 20.8306601785,-89.67204547 20.8388460306,-89.664849651 20.8413880842,-89.6598956805 20.8485923295,-89.6579042394 20.8561926932,-89.6704337521 20.8605823042,-89.6780978375 20.8779319538,-89.6791756295 20.8826596253,-89.69348393 20.8890852547,-89.6888048894 20.9019279937,-89.6798250025 20.9094271566,-89.6858002563 20.9202587142,-89.6989811754 20.9256666572,-89.7103870631 20.9189464115,-89.7110304142 20.9208494505,-89.712676599 20.9226026348,-89.7147148496 20.9133478584,-89.7313020243 20.9356925305,-89.7546868889 20.9829179179,-89.7381090163 20.9986647279,-89.729487739 21.0190758895,-89.7214339689 21.0347214094,-89.7641638638 21.1008705868,-89.7301215117 21.1736304389,-89.6789257827 21.1526816617,-89.6590806942 21.1454855962,-89.6541281932 21.1541566094,-89.6473347737 21.1686560553,-89.6420220094 21.1820061868,-89.6050631916 21.185129499,-89.5581848021 21.1591233553,-89.5504918058 21.1318386163,-89.5436526047 21.1157931405,-89.5603294053 21.0873491505,-89.5626131743 21.0691481578,-89.5486807542 21.0608497442,-89.5359415494 21.051064318,-89.5145052616 21.0441816521,-89.50660196 21.0265824363,-89.4963971858 21.032536721,-89.4671928337 21.031868919,-89.4731264674 21.0199914001,-89.491727691 21.0052102358,-89.4923715383 20.9736082179,-89.4719914774 20.96565571,-89.4857281705 20.9622817036,-89.5046478479 20.9763668184,-89.5237870952 20.981804659,-89.5434177791 20.9808817858,-89.5521430612 20.9811364779,-89.5581868011 20.9810389744,-89.5583976058 20.9777274612,-89.5561050716 20.9757315378))"
+
+#footprint = "Intersects(POLYGON((-89.79030210118333 21.122657323983717,-89.77308220413153 21.122657323983717,-89.77308220413153 21.140540053466538,-89.79030210118333 21.140540053466538,-89.79030210118333 21.122657323983717)))"
+products = sentinel.getProducts(area, ('20150101', '20180517'), {"platformname":"Sentinel-2", "cloudcoverpercentage":"[0 TO 30]"})
+#print len(products)
+#for p in products:
+#	print products[p]['filename']
+#print len(products)
 products=sentinel.filterProducts(products)
-for p in products:
-        print products[p]['filename']
+#for p in products:
+#        print products[p]['filename']
 print len(products)
 #sentinel.downloadProducts(products,dir)

tools/producListToNDVI.sh

+#!/bin/sh
+# Get Normalized Vegetation Index (NIR-RED)/(NIR+RED) for all products in PRODCUTSDIR
+
+
+

tools/productListToIndex.sh

 #!/bin/sh
-# Get Normalized index (B1-B2)/(B1+B2) for all products in PRODCUTSDIR
+# Get Normalized Index (B1-B2)/(B1+B2) for all products in PRODCUTSDIR
 PRODCUTSDIR=$1	#Products directory
 JP2DIR=$2	#JP2 Directory
 WINDOW=$3	#Crop Window