Le portail francophone de la géomatique

Jared Morgan, Write the Docs podcast host

forbsdates[2020-01-01:2020-07-31]

forbsdates[2020-02-01,2020-03-11,2020-04-21]

forbsdates[2010-01-01:2010-06-31, 2010-08-01:2010-08-31, 2010-10-01:2010-12-31]

forbandsets[1:3]

"bandset#b1"

"bandset1b1"

max("bandset#b*")

min("bandset1b*")

mean("bandset*b1")

median("bandset{2019-01-01,2019-07-31}b1")

( "#NIR#" - "#RED#" ) / ( "#NIR#" + "#RED#" )
forbandsets[1:3]
"bandset#b1" @#BANDSET#

"bandset1b1" @calc1
"bandset2b1" @calc2
"bandset3b1" @calc3
"raster1" * 2 @first_calculation
"raster1" * 2 @/home/user@first_calculation

"raster1" * 2 @temp@first_calculation
"raster1" @first_calculation.vrt

<?xml version="1.0"?>
<StyledLayerDescriptor xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/sld [schemas.opengis.net] version="1.0.0">
  <NamedLayer>
    <Name>InputLayer</Name>
    <UserStyle>
      <Title>Time comparison</Title>
      <FeatureTypeStyle>
        <Transformation>
          <ogc:Function name="ras:Jiffle">
            <ogc:Function name="parameter">
              <ogc:Literal>coverage</ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>script</ogc:Literal>
              <ogc:Literal>
              before = src[0];
              after = src[1];
              dest[0] = before;
              dest[1] = after;
              dest[2] = after - before;
            </ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>bandNames</ogc:Literal>
              <ogc:Literal>before,after,difference</ogc:Literal>
            </ogc:Function>
          </ogc:Function>
        </Transformation>
        <Rule>
          <RasterSymbolizer>
            <Opacity>1.0</Opacity>
            <ChannelSelection>
              <GrayChannel>
                <SourceChannelName>3</SourceChannelName>
              </GrayChannel>
            </ChannelSelection>
            <ColorMap>
              <ColorMapEntry color="#800000" quantity="-1"/>
              <ColorMapEntry color="#008000" quantity="1"/>
            </ColorMap>
          </RasterSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor>

{ "version": 8,
  "name": "landmarks", "layers": [
      {   "id": "Expression example",
          "type": "circle", "paint": {
              "circle-color": "#dd9933", "circle-radius":
              ["/",["get","LAND"],15],
          }
      }
  ]
}

table {
border-collapse: collapse;
}

table, th, td {
border: 1px solid black;
}
gdal_translate inputfile.tif cogfile.tif -co TILED=YES -co COPY_SRC_OVERVIEWS=YES -co COMPRESS=DEFLATE
$ pip install python-openstackclient
$ source myproject-openrc.sh
Please enter your OpenStack Password for project myproject as user myuser:
**********
$ openstack token issue                                 
+------------+-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+
| Field      | Value                                                                                                                                                                                   |
+------------+-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+
| expires    | 2020-07-21T08:15:12+0000                                                                                                                                                                |
| id         | xxxx_my_token_xxxx
| project_id | 97e2e750f1904b41b76f80a50dabde0a                                                                                                                                                        |
| user_id    | 18f7ccaf1a2d4344a4e35f0d84eb065e                                                                                                                                                        |
+------------+-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------+
openstack object create MyContainer cogfile.tif --name cogfile.tif
$ export SWIFT_AUTH_TOKEN=xxxx_my_token_xxxx
$ export SWIFT_STORAGE_URL=https://storage.sbg.cloud.ovh.net/v1/AUTH_$OS_PROJECT_ID
<?xml version="1.0" encoding="UTF-8"?>
<StyledLayerDescriptor xmlns="http://www.opengis.net/sld"
 xmlns:ogc="http://www.opengis.net/ogc"
 xmlns:xlink="http://www.w3.org/1999/xlink"
 xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
 xsi:schemaLocation="http://www.opengis.net/sld [schemas.opengis.net]  version="1.0.0">
  <NamedLayer>
    <Name>InputLayer</Name>
    <UserStyle>
      <Title>Time comparison</Title>
      <FeatureTypeStyle>
        <Transformation>
          <ogc:Function name="ras:Jiffle">
            <ogc:Function name="parameter">
              <ogc:Literal>coverage</ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>script</ogc:Literal>
              <ogc:Literal>
                before = src[0];
                after = src[1];
                dest[0] = before;
                dest[1] = after;
                dest[2] = after - before;
              </ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>bandNames</ogc:Literal>
              <ogc:Literal>before,after,difference</ogc:Literal>
            </ogc:Function>
          </ogc:Function>
        </Transformation>
        <Rule>
          <RasterSymbolizer>
            <Opacity>1.0</Opacity>
            <ChannelSelection>
              <GrayChannel>
                <SourceChannelName>3</SourceChannelName>
              </GrayChannel>
            </ChannelSelection>
            <ColorMap>
              <ColorMapEntry color="#800000" quantity="-1"/>
              <ColorMapEntry color="#008000" quantity="1"/>
            </ColorMap>
          </RasterSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor>
{ "version": 8,
  "name": "landmarks",
  "layers": [
      {   "id": "Expression example",
          "type": "circle",
          "paint": {
              "circle-color": "#dd9933",
              "circle-radius": ["/",["get","LAND"],15],
          }
      }
  ]
}
<code><?xml version="1.0" encoding="UTF-8"?>
<StyledLayerDescriptor xmlns="http://www.opengis.net/sld"
 xmlns:ogc="http://www.opengis.net/ogc"
 xmlns:xlink="http://www.w3.org/1999/xlink"
 xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
 xsi:schemaLocation="http://www.opengis.net/sld [schemas.opengis.net]  version="1.0.0">
  <NamedLayer>
    <Name>InputLayer</Name>
    <UserStyle>
      <Title>Time comparison</Title>
      <FeatureTypeStyle>
        <Transformation>
          <ogc:Function name="ras:Jiffle">
            <ogc:Function name="parameter">
              <ogc:Literal>coverage</ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>script</ogc:Literal>
              <ogc:Literal>
                before = src[0];
                after = src[1];
                dest[0] = before;
                dest[1] = after;
                dest[2] = after - before;
              </ogc:Literal>
            </ogc:Function>
            <ogc:Function name="parameter">
              <ogc:Literal>bandNames</ogc:Literal>
              <ogc:Literal>before,after,difference</ogc:Literal>
            </ogc:Function>
          </ogc:Function>
        </Transformation>
        <Rule>
          <RasterSymbolizer>
            <Opacity>1.0</Opacity>
            <ChannelSelection>
              <GrayChannel>
                <SourceChannelName>3</SourceChannelName>
              </GrayChannel>
            </ChannelSelection>
            <ColorMap>
              <ColorMapEntry color="#800000" quantity="-1"/>
              <ColorMapEntry color="#008000" quantity="1"/>
            </ColorMap>
          </RasterSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor></code>

<code>{ "version": 8,
  "name": "landmarks",
  "layers": [
      {   "id": "Expression example",
          "type": "circle",
          "paint": {
              "circle-color": "#dd9933",
              "circle-radius": ["/",["get","LAND"],15],
          }
      }
  ]
}</code>

import ee
ee.Authenticate()
ee.Initialize()

area_estudo = ee.Geometry.Rectangle([-49.7, -28.3, -49.3, -28.7])
img_landsat = ee.Image('LANDSAT/LC08/C01/T1/LC08_220080_20200420')\
    .select(['B2','B3','B4']).clip(area_estudo)

print(img_landsat.reduceRegion(reducer = ee.Reducer.minMax()).getInfo())

print(img_landsat.getThumbURL({'min':5000, 'max': 23000}))

task = ee.batch.Export.image.toDrive(**{
    'image': img_landsat,
    'description': 'img_B2E',
    'folder': 'EE_Output',
    'scale': 30,
    'region': area_estudo.getInfo()['coordinates']
    })
task.start()

import time
while task.active():
    print('Salvando a imagem (id: {}).'.format(task.id))
    time.sleep(5)

download_url = img_landsat.getDownloadURL({
    'image': img_landsat,
    'region': area_estudo,
    'name': 'img_B2E_down',
})

print(download_url)

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Carte des membres Sites amis	Annuaire Agenda Biblio Netvibes Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Carte des CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Feeds

Autres planets