//
// Make a polyline for this vtLine by using the points from each node.
//
void vtLine::MakePolyline(DLine2 &polyline)
{
const uint num = m_poles.size();
polyline.SetSize(num);
for (uint i = 0; i < num; i++)
polyline[i] = m_poles[i]->m_p;
}
void vtVegLayer::AddElementsFromLULC(vtLULCFile *pLULC)
{
LULCSection *section;
LULCPoly *poly;
SetVegType(VLT_Density);
//set projections
vtProjection proj_new;
proj_new.SetProjectionSimple(0, -1, EPSG_DATUM_WGS84);
SetProjection(proj_new);
// figure out the number of polygons in file
uint size = 0;
for (uint sec = 0; sec < pLULC->NumSections(); sec++)
{
section = pLULC->GetSection(sec);
size = size + section->m_iNumPolys;
}
// Create density field
m_field_density = m_pSet->AddField("Density", FT_Float);
m_pSet->SetNumEntities(size);
// get each poly from LULC file
uint i, s, p, count = 0;
float density=0;
for (s = 0; s < pLULC->NumSections(); s++)
{
section = pLULC->GetSection(s);
for (p = 0; p < section->m_iNumPolys; p++)
{
poly = section->m_pPoly + p;
bool wild = false;
switch (poly->Attribute)
{
case 42: // forest
wild = true;
density = 1.0f;
break;
case 32:
case 33:
wild = true;
density = 0.5;
break;
case 22: // orchards
wild = false;
// no crops for now
break;
default:
density = 0.0f;
break;
}
DLine2 dline;
dline.SetSize(poly->m_iCoords);
// get Coords of LULCpoly and store as latlon, then save in VPoly
for (i = 0; i < dline.GetSize(); i++)
dline.SetAt(i, poly->m_p[i]);
DPolygon2 dpoly;
dpoly.push_back(dline);
GetPS()->SetPolygon(count, dpoly);
m_pSet->SetValue(count, m_field_density, density);
count++;
}
}
}
