Changeset 9206


Ignore:
Timestamp:
Jun 19, 2014, 8:00:56 PM (11 years ago)
Author:
steve
Message:

Adding in the flood extent figure

Location:
trunk/anuga_core/demos/cairns
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/anuga_core/demos/cairns/cairns_fixed_wave.qgs

    r9204 r9206  
    2424        </destinationsrs>
    2525    </mapcanvas>
    26     <legend updateDrawingOrder="true" activeLayer="cairns_fixed_wave_elevation_max20140618233339039">
     26    <legend updateDrawingOrder="true" activeLayer="cairns_fixed_wave_stage_max20140618233338237">
    2727        <legendlayer drawingOrder="-1" open="true" checked="Qt::Checked" name="coast_contour" showFeatureCount="0">
    2828            <filegroup open="true" hidden="false">
    … …  
    3030            </filegroup>
    3131        </legendlayer>
    32         <legendlayer drawingOrder="-1" open="true" checked="Qt::Unchecked" name="flood" showFeatureCount="0">
     32        <legendlayer drawingOrder="-1" open="true" checked="Qt::Checked" name="flood" showFeatureCount="0">
    3333            <filegroup open="true" hidden="false">
    34                 <legendlayerfile isInOverview="0" layerid="flood20140618233428952" visible="0"/>
     34                <legendlayerfile isInOverview="0" layerid="flood20140618233428952" visible="1"/>
    3535            </filegroup>
    3636        </legendlayer>
    37         <legendlayer drawingOrder="-1" open="true" checked="Qt::Checked" name="cairns_fixed_wave_velocity_max" showFeatureCount="0">
     37        <legendlayer drawingOrder="-1" open="true" checked="Qt::Checked" name="max_velocity" showFeatureCount="0">
    3838            <filegroup open="true" hidden="false">
    3939                <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_velocity_max20140618233225048" visible="1"/>
    4040            </filegroup>
    4141        </legendlayer>
    42         <legendlayer drawingOrder="-1" open="true" checked="Qt::Unchecked" name="cairns_fixed_wave_elevation_max" showFeatureCount="0">
    43             <filegroup open="true" hidden="false">
    44                 <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_elevation_max20140618233339039" visible="0"/>
     42        <legendlayer drawingOrder="-1" open="false" checked="Qt::Checked" name="elevation" showFeatureCount="0">
     43            <filegroup open="false" hidden="false">
     44                <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_elevation_max20140618233339039" visible="1"/>
    4545            </filegroup>
    4646        </legendlayer>
    47         <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="cairns_fixed_wave_depth_max" showFeatureCount="0">
     47        <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="max_wave_depth" showFeatureCount="0">
    4848            <filegroup open="false" hidden="false">
    4949                <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_depth_max20140618233338702" visible="0"/>
    5050            </filegroup>
    5151        </legendlayer>
    52         <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="cairns_fixed_wave_stage_max" showFeatureCount="0">
     52        <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="max_stage" showFeatureCount="0">
    5353            <filegroup open="false" hidden="false">
    5454                <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_stage_max20140618233338237" visible="0"/>
    5555            </filegroup>
    5656        </legendlayer>
    57         <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="cairns_fixed_wave_stage_0.0" showFeatureCount="0">
     57        <legendlayer drawingOrder="-1" open="false" checked="Qt::Unchecked" name="initial stage" showFeatureCount="0">
    5858            <filegroup open="false" hidden="false">
    5959                <legendlayerfile isInOverview="0" layerid="cairns_fixed_wave_stage_0_020140618233337613" visible="0"/>
    … …  
    6161        </legendlayer>
    6262    </legend>
    63     <Composer title="cairns_elevation" visible="1">
     63    <Composer title="cairns_elevation" visible="0">
    6464        <Composition printResolution="300" paperWidth="297" snapGridOffsetX="0" snapGridResolution="10" snapGridOffsetY="0" alignmentSnap="1" snapping="0" numPages="1" paperHeight="210" alignmentSnapTolerance="2" printAsRaster="0">
    6565            <ComposerLegend symbolWidth="7" title="Legend" boxSpace="2" splitLayer="0" equalColumnWidth="0" map="0" symbolHeight="4" wrapChar="" fontColor="#000000" columnSpace="2" columnCount="1">
    … …  
    8585                    </LayerItem>
    8686                </Model>
    87                 <ComposerItem width="74.8531" x="3.00815" transparency="0" y="116.988" rotation="0" height="79.4" positionMode="0" frame="false" outlineWidth="0.3" zValue="4" lastValidViewScaleFactor="2.32701" uuid="{9ca87187-a045-47f6-a771-c4d02b008bc2}" id="" background="true" blendMode="0" positionLock="false">
     87                <ComposerItem width="74.8531" x="3.00815" transparency="0" y="116.988" rotation="0" height="79.4" positionMode="0" frame="false" outlineWidth="0.3" zValue="4" lastValidViewScaleFactor="2.78673" uuid="{9ca87187-a045-47f6-a771-c4d02b008bc2}" id="" background="true" blendMode="0" positionLock="false">
    8888                    <FrameColor alpha="255" red="0" blue="0" green="0"/>
    8989                    <BackgroundColor alpha="255" red="255" blue="255" green="255"/>
    … …  
    137137        </Composition>
    138138    </Composer>
    139     <Composer title="cairns_flood" visible="1">
     139    <Composer title="cairns_flood" visible="0">
    140140        <Composition printResolution="300" paperWidth="297" snapGridOffsetX="0" snapGridResolution="10" snapGridOffsetY="0" alignmentSnap="1" snapping="0" numPages="1" paperHeight="210" alignmentSnapTolerance="2" printAsRaster="0">
    141141            <ComposerLegend symbolWidth="7" title="Legend" boxSpace="2" splitLayer="0" equalColumnWidth="0" map="0" symbolHeight="4" wrapChar="" fontColor="#000000" columnSpace="2" columnCount="1">
    … …  
    210210        </Composition>
    211211    </Composer>
    212     <Composer title="cairns_velocity" visible="1">
     212    <Composer title="cairns_velocity" visible="0">
    213213        <Composition printResolution="300" paperWidth="297" snapGridOffsetX="0" snapGridResolution="10" snapGridOffsetY="0" alignmentSnap="1" snapping="0" numPages="1" paperHeight="210" alignmentSnapTolerance="2" printAsRaster="0">
    214214            <ComposerLegend symbolWidth="7" title="Legend" boxSpace="2" splitLayer="0" equalColumnWidth="0" map="0" symbolHeight="4" wrapChar="" fontColor="#000000" columnSpace="2" columnCount="1">
    … …  
    234234                </ComposerItem>
    235235            </ComposerLegend>
    236             <ComposerMap keepLayerSet="true" overviewFrameMap="-1" id="0" previewMode="Cache" drawCanvasItems="false">
     236            <ComposerMap keepLayerSet="false" overviewFrameMap="-1" id="0" previewMode="Cache" drawCanvasItems="true">
    237237                <overviewFrame overviewBlendMode="0" overviewFrameMap="-1" overviewInverted="false">
    238238                    <symbol alpha="0.3" type="fill" name="">
    … …  
    250250                </overviewFrame>
    251251                <Extent ymin="7.97918e+06" xmin="291454" ymax="8.30513e+06" xmax="721654"/>
    252                 <LayerSet>
    253                     <Layer>coast_contour20140619121652371</Layer>
    254                     <Layer>cairns_fixed_wave_velocity_max20140618233225048</Layer>
    255                 </LayerSet>
     252                <LayerSet/>
    256253                <Grid gridStyle="0" crossLength="3" intervalX="0" offsetX="0" intervalY="0" offsetY="0" gridFrameStyle="0" gridFrameWidth="2" gridBlendMode="0" show="1">
    257254                    <symbol alpha="1" type="line" name="">
    … …  
    282279    </Composer>
    283280    <projectlayers layercount="7">
    284         <maplayer minimumScale="0" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
     281        <maplayer minimumScale="-4.65661e-10" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
    285282            <id>cairns_fixed_wave_depth_max20140618233338702</id>
    286283            <datasource>./cairns_fixed_wave_depth_max.tif</datasource>
    … …  
    290287                <value></value>
    291288            </keywordList>
    292             <layername>cairns_fixed_wave_depth_max</layername>
     289            <layername>max_wave_depth</layername>
    293290            <srs>
    294291                <spatialrefsys>
    … …  
    308305            <provider>gdal</provider>
    309306            <pipe>
    310                 <rasterrenderer gradient="BlackToWhite" opacity="1" alphaBand="-1" type="singlebandgray" grayBand="1">
     307                <rasterrenderer gradient="BlackToWhite" opacity="1" alphaBand="0" type="singlebandgray" grayBand="1">
    311308                    <rasterTransparency/>
    312309                    <contrastEnhancement>
    … …  
    322319            <blendMode>0</blendMode>
    323320        </maplayer>
    324         <maplayer minimumScale="-4.65661e-10" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
     321        <maplayer minimumScale="0" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
    325322            <id>cairns_fixed_wave_elevation_max20140618233339039</id>
    326323            <datasource>./cairns_fixed_wave_elevation_max.tif</datasource>
    … …  
    330327                <value></value>
    331328            </keywordList>
    332             <layername>cairns_fixed_wave_elevation_max</layername>
     329            <layername>elevation</layername>
    333330            <srs>
    334331                <spatialrefsys>
    … …  
    348345            <provider>gdal</provider>
    349346            <pipe>
    350                 <rasterrenderer opacity="1" alphaBand="0" classificationMax="100" classificationMinMaxOrigin="User" band="1" classificationMin="-100" type="singlebandpseudocolor">
     347                <rasterrenderer opacity="1" alphaBand="-1" classificationMax="100" classificationMinMaxOrigin="User" band="1" classificationMin="-100" type="singlebandpseudocolor">
    351348                    <rasterTransparency/>
    352349                    <rastershader>
    … …  
    371368            <blendMode>0</blendMode>
    372369        </maplayer>
    373         <maplayer minimumScale="0" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
     370        <maplayer minimumScale="-4.65661e-10" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
    374371            <id>cairns_fixed_wave_stage_0_020140618233337613</id>
    375372            <datasource>./cairns_fixed_wave_stage_0.0.tif</datasource>
    … …  
    379376                <value></value>
    380377            </keywordList>
    381             <layername>cairns_fixed_wave_stage_0.0</layername>
     378            <layername>initial stage</layername>
    382379            <srs>
    383380                <spatialrefsys>
    … …  
    397394            <provider>gdal</provider>
    398395            <pipe>
    399                 <rasterrenderer gradient="BlackToWhite" opacity="1" alphaBand="-1" type="singlebandgray" grayBand="1">
     396                <rasterrenderer gradient="BlackToWhite" opacity="1" alphaBand="0" type="singlebandgray" grayBand="1">
    400397                    <rasterTransparency/>
    401398                    <contrastEnhancement>
    … …  
    419416                <value></value>
    420417            </keywordList>
    421             <layername>cairns_fixed_wave_stage_max</layername>
     418            <layername>max_stage</layername>
    422419            <srs>
    423420                <spatialrefsys>
    … …  
    451448            <blendMode>0</blendMode>
    452449        </maplayer>
    453         <maplayer minimumScale="-4.65661e-10" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
     450        <maplayer minimumScale="0" maximumScale="1e+08" type="raster" hasScaleBasedVisibilityFlag="0">
    454451            <id>cairns_fixed_wave_velocity_max20140618233225048</id>
    455452            <datasource>./cairns_fixed_wave_velocity_max.tif</datasource>
    … …  
    459456                <value></value>
    460457            </keywordList>
    461             <layername>cairns_fixed_wave_velocity_max</layername>
     458            <layername>max_velocity</layername>
    462459            <srs>
    463460                <spatialrefsys>
    … …  
    477474            <provider>gdal</provider>
    478475            <pipe>
    479                 <rasterrenderer opacity="1" alphaBand="0" classificationMax="14.4888" classificationMinMaxOrigin="CumulativeCutFullExtentEstimated" band="1" classificationMin="0" type="singlebandpseudocolor">
     476                <rasterrenderer opacity="1" alphaBand="-1" classificationMax="14.4888" classificationMinMaxOrigin="CumulativeCutFullExtentEstimated" band="1" classificationMin="0" type="singlebandpseudocolor">
    480477                    <rasterTransparency/>
    481478                    <rastershader>
    … …  
    497494        <maplayer minimumScale="-4.65661e-10" maximumScale="1e+08" minLabelScale="0" maxLabelScale="1e+08" geometry="Line" type="vector" hasScaleBasedVisibilityFlag="0" scaleBasedLabelVisibilityFlag="0">
    498495            <id>coast_contour20140619121652371</id>
    499             <datasource>coast_contour</datasource>
     496            <datasource>./coast_contour</datasource>
    500497            <title></title>
    501498            <abstract></abstract>
    … …  
    686683                <edittype labelontop="0" editable="1" type="0" name="ID"/>
    687684            </edittypes>
    688             <editform></editform>
     685            <editform>.</editform>
    689686            <editforminit></editforminit>
    690             <annotationform></annotationform>
     687            <annotationform>.</annotationform>
    691688            <editorlayout>generatedlayout</editorlayout>
    692689            <excludeAttributesWMS/>

  • trunk/anuga_core/demos/cairns/export_results_to_qgis.py

    r9198 r9206  
    1010print 'output dir:', name
    1111
    12 #for which_var in [0,2,3,4,5]:
    13 for which_var in []:
     12"""
     13Produce .tif files extracting results of cairns simulation.
    1414
    15 
    16 
    17     if which_var == 0:    # Stage
    18         outname = name + '_stage'
    19         quantityname = 'stage'
    20         reduction = max
    21 
    22     if which_var == 5:    # Stage
    23         outname = name + '_stage_0'
    24         quantityname = 'stage'
    25         reduction = 0   
    26 
    27     if which_var == 1:    # Absolute Momentum
    28         outname = name + '_momentum'
    29         quantityname = '(xmomentum**2 + ymomentum**2)**0.5'    #Absolute momentum
    30         reduction = max
    31 
    32     if which_var == 2:    # Depth
    33         outname = name + '_depth'
    34         quantityname = 'stage-elevation'  #Depth
    35         reduction = max
    36 
    37     if which_var == 3:    # Speed
    38         outname = name + '_speed'
    39         quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-3)'  #Speed
    40         reduction = max
    41 
    42     if which_var == 4:    # Elevation
    43         outname = name + '_elevation'
    44         quantityname = 'elevation'  #Elevation
    45         reduction = max
    46 
    47     print 72*'='
    48     print 'sww2dem '+outname
    49     print 72*'='   
    50 
    51     anuga.sww2dem(name+'.sww',
    52             outname+'.asc',
    53             quantity=quantityname,
    54             cellsize=100,     
    55             easting_min=project.eastingmin,
    56             easting_max=project.eastingmax,
    57             northing_min=project.northingmin,
    58             northing_max=project.northingmax,       
    59             reduction=reduction,
    60             verbose=True)
     15Can be viewed with qgis (a project is available cairns_fixed_wave.qgs)
     16"""
    6117   
    6218from anuga.utilities.plot_utils import Make_Geotif
Note: See TracChangeset for help on using the changeset viewer.