~ubuntu-branches/ubuntu/hardy/libterralib/hardy

result = operation->compute (stos.begin(input_attrname), stos.end(input_attrname), output_columnName); // if property not found, stos iterator will return zero to operation. To check this, have to be less general and use getProperty inside the compute method.

799

800

TeFillCellInitSTO ((*cell_it), t, result, cellObjSet);

801

cell_it++;

802

803

}

804

805

// Update DB

806

if (!TeUpdateDBFromSet (&cellObjSet, cell_tablename))

807

return false;

808

return true;

809

}

810

811

812

813

814

815

bool TeFillCellNonSpatialRasterOperation (TeDatabase* db,

816

const string& input_raster_layername,

817

TeTimeInterval t,

818

const string& cell_layername,

819

const string& cell_tablename,

820

const string& output_columnName,

821

TeComputeAttrStrategy<TeRaster::iteratorPoly>* operation)

822

{

823

if (!db) return false;

824

825

// Load input layers

826

TeLayer* input_layer = new TeLayer (input_raster_layername);

827

if (!db->loadLayer (input_layer))

828

{

829

cout << "\tLayer de entrada inexistente: " << input_raster_layername << endl;

830

db->close();

831

return false;

832

}

833

834

TeRaster* raster = input_layer->raster();

835

if (!raster) return false;

836

837

838

// Load output cells layer with geometry previously created

839

TeLayer* cell_layer = new TeLayer (cell_layername);

840

if (!db->loadLayer (cell_layer))

841

{

842

cout << "\tLayer de entrada inexistente: " << cell_layername << endl;

843

db->close();

844

return false;

845

}

846

847

848

// Initialize cell set

849

TeCellSet cells;

850

if (!TeFillCellInitLoad (cell_layer, cell_tablename, cells)) return false;

851

852

853

// Initialize object set to store cell properties

854

TeSTElementSet cellObjSet (cell_layer);

855

856

857

// Process

858

TePropertyVector result;

859

TeCellSet::iterator cell_it = cells.begin();

860

while (cell_it != cells.end())

861

{

862

TePolygon p = TeBox2Polygon((*cell_it).box());

863

TeRaster::iteratorPoly rasterItBegin(raster);

864

rasterItBegin = raster->begin(p, TeBBoxPixelInters, 0);

865

TeRaster::iteratorPoly rasterItEnd(raster);

866

rasterItEnd = raster->end(p, TeBBoxPixelInters, 0);

867

868

result = operation->compute (rasterItBegin, rasterItEnd, output_columnName);

869

TeFillCellInitSTO ((*cell_it), t, result, cellObjSet);

870

cell_it++;

871

}

872

873

// Update DB

874

if (!TeUpdateDBFromSet (&cellObjSet, cell_tablename))

875

return false;

876

return true;

877

}

878

879

880

881

bool TeFillCellAggregateOperation ( TeDatabase* db,

882

const string& input_layername,

883

const string& input_tablename,

884

TeTimeInterval t,

885

const string& cell_layername,

886

const string& cell_tablename,

887

vector<string>& attrNames,

888

TeComputeAttrStrategy<TeSTElementSet::propertyIterator>* operation)

889

{

890

if (!db) return false;

891

892

893

// Load input layers

894

TeLayer* input_layer = new TeLayer (input_layername);

895

if (!db->loadLayer (input_layer))

896

{

897

cout << "\tLayer de entrada inexistente: " << input_layername << endl;

898

db->close();

899

return false;

900

}

901

902

903

// Load output cells layer with geometry previously created

904

TeLayer* cell_layer = new TeLayer (cell_layername);

905

906

if (!db->loadLayer (cell_layer))

907

{

908

cout << "\tLayer de entrada inexistente: " << cell_layername << endl;

909

db->close();

910

return false;

911

}

912

913

// Initialize cell set

914

TeCellSet cells;

915

if (!TeFillCellInitLoad (cell_layer, cell_tablename, cells)) return false;

916

917

918

// Initialize object set to store cell properties

919

TeSTElementSet cellObjSet (cell_layer);

920

921

// Initialize theme

922

923

TeTheme* theme = new TeTheme ("", input_layer);

924

vector<string> attrTableNames;

925

attrTableNames.push_back (input_tablename);

926

TeAttrTableVector atts;

927

if (!input_layer->getAttrTablesByName(attrTableNames, atts)) return false;

928

theme->setAttTables (atts);

929

930

931

// Process

932

TeCellSet::iterator cell_it = cells.begin();

933

while (cell_it != cells.end())

934

{

935

TeSTInstance cellObj;

936

cellObj.objectId ((*cell_it).objectId());

937

string uniqueId = (*cell_it).objectId();

938

cellObj.addUniqueId (uniqueId);

939

uniqueId += t.getInitialDate() + t.getInitialTime();

940

uniqueId += t.getFinalDate() + t.getFinalTime() ;

941

cellObj.addUniqueId (uniqueId);

942

cellObj.timeInterval (t);

943

944

// Set restrictions on a theme and create stoset

945

theme->setSpatialRest((*cell_it).box(), TeCELLS); //????? TeWITHIN ok? ANAP

946

TeSTElementSet stos (theme);

947

TeSTOSetBuildDB (&stos, true, false, attrNames);

948

949

// Create stoset based on input layer

950

TePropertyVector result;

951

vector<string>::iterator attIt = attrNames.begin();

952

while (attIt != attrNames.end())

953

{

954

result = operation->compute (stos.begin(*attIt), stos.end(*attIt), (*attIt)); // if property not found, stos iterator will return zero to operation. To check this, have to be less general and use getProperty inside the compute method.

955

TePropertyVector::iterator itProp = result.begin();

956

while (itProp != result.end())

957

{

958

cellObj.addProperty (*itProp);

959

itProp++;

960

}

961

attIt++;

962

}

963

cellObjSet.insertSTInstance (cellObj);

964

965

cell_it++;

966

}

967

968

// Update DB

969

if (!TeUpdateDBFromSet (&cellObjSet, cell_tablename))

970

return false;

971

972

return true;

973

}

974

975

976

977

978

979

980

bool TeFillCellConnectionOperation ( TeDatabase* db,

981

TeSTElementSet objects,

982

TeGeomRep rep,

983

TeGraphNetwork* net,

984

const string& input_attrName,

985

986

TeTimeInterval t,

987

const string& cell_layername,

988

const string& cell_tablename,

989

const string& output_columnName,

990

991

992

TeComputeAttrStrategy<TeSTElementSet::propertyIterator>* operation,

993

double local_distance_factor,

994

double net_distance_factor,

995

double net_conn_factor,

996

double mult_factor)

997

{

998

999

if (!operation) return false;

1000

if (!net) return false;

1001

if (!db) return false;

1002

1003

// Load output cells layer with geometry previously created

1004

1005

TeLayer* cell_layer = new TeLayer (cell_layername);

1006

1007

if (!db->loadLayer (cell_layer))

1008

{

1009

cout << "\tLayer de entrada inexistente: " << cell_layername << endl;

1010

db->close();

1011

return false;

1012

}

1013

1014

// Initialize cell set

1015

TeCellSet cells;

1016

if (!TeFillCellInitLoad (cell_layer, cell_tablename, cells)) return false;

1017

1018

1019

// Initialize object set to store cell properties

1020

TeSTElementSet cellObjSet (cell_layer);

1021

1022

1023

// Process

1024

TePropertyVector result;

1025

TeCellSet::iterator cell_it = cells.begin();

1026

while (cell_it != cells.end())

1027

{

1028

TeSTInstance cellObj;

1029

cellObj.objectId ((*cell_it).objectId());

1030

string uniqueId = (*cell_it).objectId();

1031

cellObj.addUniqueId (uniqueId);

1032

uniqueId += t.getInitialDate() + t.getInitialTime();

1033

uniqueId += t.getFinalDate() + t.getFinalTime() ;

1034

cellObj.addUniqueId (uniqueId);

1035

cellObj.timeInterval (t);

1036

cellObjSet.insertSTInstance (cellObj);

1037

cell_it++;

1038

}

1039

1040

TeProxMatrixOpenNetworkStrategy2 sc_net (&cellObjSet, TeCELLS, &objects, rep, 0, TeMAXFLOAT, TeMAXFLOAT, net);

1041

// TeProxMatrixConnectionStrenghtStrategy sw(4, 0.0, false);

1042

TeProxMatrixInverseDistanceStrategy sw (local_distance_factor, net_distance_factor, net_conn_factor, mult_factor, false);

1043

TeProxMatrixNoSlicingStrategy ss_no;

1044

TeGeneralizedProxMatrix mat (&sc_net, &sw, &ss_no);

1045

1046

// initialize object set to store cell properties

1047

TeSTElementSet cellObjSet2 (cell_layer);

1048

TeCellSet::iterator cell_it2 = cells.begin();

1049

while (cell_it2 != cells.end())

1050

{

1051

TeSTElementSet neigh = mat.getSTENeighbours((*cell_it2).objectId());

1052

1053

result = operation->compute (neigh.begin(input_attrName), neigh.end(input_attrName), output_columnName); // if property not found, stos iterator will return zero to operation. To check this, have to be less general and use getProperty inside the compute method.

1054

TeFillCellInitSTO ((*cell_it2), t, result, cellObjSet2);

1055

cell_it2++;

1056

}

1057

1058

1059

// Update DB

1060

if (!TeUpdateDBFromSet (&cellObjSet2, cell_tablename))

1061

return false;

1062

1063

return true;

1064

}

1065

1066

Older »