Hi

you have to use tcl for that ( or maybe phyton)  there is no standard 
transformer in WB,

attatched is a mapping file transformer and a tcl script doing this - and a 
lot more, maybe you can get some idearts from it

what I do is  I save the geometry as a ogc well known text attribute attract 
the z values off all points, the I use the SurfaceModelFactory to calculate 
the ground level of the building, repeat step one, now i can calculate the 
height etc on the buildings

If you just need teh std dev. you can use the FME_Coordinates function to 
get all z values


NOTE from a statistical point off view it only give sens to calculate the 
std. dev. to the buildings iff the z values are on the same 'level' - they 
are representing a normal distribution

NB maybe it could be don using the statitstis calculater - split a copy of 
the buildings into individual points with a uniqu id on all buildings 
calculate std and the join it back to the original building with a 
ReferanceFactory



have fun
Peter Laulund


>From: "cityoftoronto_mapping" <[EMAIL PROTECTED]>
>Reply-To: [email protected]
>To: [email protected]
>Subject: [fme] Standard deviation of z values (elevation) for individual 
>shapes in a  DGN
>Date: Mon, 17 Apr 2006 01:30:51 -0000
>
>Wondering if anyone has come across a workflow for this... I have a
>Microstation (DGN) file that is comprised of 3D shapes (building
>outlines).  Each 3D building outline (shape) can consist of unique
>coordinates (in particular the z value).  Im looking to calculate the
>std. dev. of the z value for each individual shape, and was wondering
>which engine would work I ... I've tried 'Statistical' with no luck.
>In total, i'm looking at approximately 500,000 shapes.  Any ideas?
>Btw, the output file format is an ESRI shapefile.
>
>
>
>
>
>
>
>
>
>Get the maximum benefit from your FME, FME Objects, or SpatialDirect via 
>our Professional Services team.  Visit www.safe.com/services for details.
>Yahoo! Groups Links
>
>
>
>
>
>

_________________________________________________________________
Opret en personlig blog og del dine billeder på MSN Spaces:  
http://spaces.msn.com


Get the maximum benefit from your FME, FME Objects, or SpatialDirect via our 
Professional Services team.  Visit www.safe.com/services for details. 
Yahoo! Groups Links

<*> To visit your group on the web, go to:
    http://groups.yahoo.com/group/fme/

<*> To unsubscribe from this group, send an email to:
    [EMAIL PROTECTED]

<*> Your use of Yahoo! Groups is subject to:
    http://docs.yahoo.com/info/terms/
 

Attachment: GetStat.TCL
Description: Binary data

/* ============================================================================

  Navn          : Peter Laulund, KMS

  Oprettet dato : 30 januar 2004 - 10:34

  Beskrivelse   :

  Problemer     : Hvor stort et område ekstra skal med for at være sikker
                  på at der ikke er fejl i kanten + 1500 m




============================================================================*/

# ============================================================================= # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - FACTORY_DEF * TestFactory \ FACTORY_NAME "GENERATE 3D: Soter input data " \ INPUT FEATURE_TYPE * \ @SupplyAttributes(datatype, @Lookup(SortData, &OBJEKTKODE)) \ @FeatureType(&datatype) \ @Count(@Concatenate(&OBJEKTTYPE, " - ", &OBJEKTKODE), 1) \ TEST &datatype != no \ BOOLEAN_OPERATOR AND \ OUTPUT PASSED FEATURE_TYPE * \
     @KeepAttributes(fme_type, fme_geometry, OBJEKTKODE, OBJEKTTYPE, ID_NR)

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - FACTORY_DEF * TeeFactory \ FACTORY_NAME "GENERATE 3D: Tilføj attributter til Bygninger" \ INPUT FEATURE_TYPE byg \ @RoundOffCoords(X,2) @RoundOffCoords(Y,2) \ OUTPUT FEATURE_TYPE * \ @OGCGeometry(to_attribute, wkt, BygGeom ) \ @SupplyAttributes(ID, @Count(id_tæller, 1)) \ @SupplyAttributes(POLY, @Evaluate( @NumHoles() + 1)) \ @SupplyAttributes(ANTAL_PKT, @Evaluate( @NumCoords() - &POLY)) \
     @Log(bygning, 2, 1)

# -----------------------------------------------------------------------------
INCLUDE GET_GRID.FMI

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# meget afhænger af de indstillinger der er her i farbrikken
FACTORY_DEF * SurfaceModelFactory \ FACTORY_NAME "GENERATE 3D: Tilføj grundniveau som Z" \ INPUT DRAPE_FEATURES FEATURE_TYPE byg \ INPUT POINTS FEATURE_TYPE ref \ SURFACE_BOUNDING_BOX $(X1) $(Y1) $(X2) $(Y2) \ INTERPOLATION_TYPE PLANAR \ DRAPE_TYPE VERTEX \
   OUTPUT DRAPED_FEATURES FEATURE_TYPE byg


# TOLERANCE 1 \

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - FACTORY_DEF * TeeFactory \ FACTORY_NAME "GENERATE 3D: Tilpas og beregn sokkel" \ INPUT FEATURE_TYPE byg \ @OGCGeometry(to_attribute, wkt, Sokkel ) \ @TCL( CalculateStat ) \ OUTPUT FEATURE_TYPE * \ @OGCGeometry(from_attribute, wkt, BygGeom ) \ @SupplyAttributes(FinalFeature, TDCBYGN) \
    @SupplyAttributes(fme_template_feature_type, TDC)

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - FACTORY_DEF * TeeFactory \ FACTORY_NAME "GENERATE 3D: Drop Z-Værdier" \ INPUT FEATURE_TYPE * \ OUTPUT FEATURE_TYPE * \
     @Force2D()

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - FACTORY_DEF * CreationFactory \ FACTORY_NAME "GENERATE 3D: Dan Boundingbox" \ 2D_GEOMETRY $(X1) $(Y1) $(X2) $(Y1) $(X2) $(Y2) $(X1) $(Y2) $(X1) $(Y1) \ COORDINATE_SYSTEM $(DestCoordSys) \ CREATE_AT_END \ NUMBER_TO_CREATE 1 \ OUTPUT FEATURE_TYPE * \ @SupplyAttributes(FinalFeature, BBOX, OBJEKTKODE, 100, OBJEKTTYPE, "Bounding box for område" )\
     @SupplyAttributes(LB_NR, 1)

# -----------------------------------------------------------------------------
Lookup SortData   \
  2655   byg      \
  3210   byg      \
  4625   byg      \
  8111   dhm      \
  3113   ref      \
  3116   ref      \
  3117   ref      \
  3118   ref      \
  2700   ref      \
  3319   ref      \
  2237   ref      \
  5500   ref      \
  2439   ref      \
  2710   ref      \
  7710   ref      \
  5158   ref      \
  ""     no

# er 2710 en fejl for 7710 dvs. kyst?
# skal vi alternativt bruge veje / dhm  eller hvad?
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # =============================================================================

Reply via email to