//Получение координат центра контура
Point getContourCentre(CONT& vec)
{
double tempx = 0.0, tempy = 0.0;
for (int i = 0; i < vec.size(); i++) {
tempx += vec[i].x;
tempy += vec[i].y;
}
return Point(tempx / (double)vec.size(), tempy / (double)vec.size());
}
//Проверка находится ли один контур в другом
bool isContourInsideContour(CONT& in, CONT& out)
{
//каждую точку внутреннего контура проверяем на нахождение во внешнем контуре
for(int i = 0; i < in.size(); i++){
if(pointPolygonTest(out, in[i], false) <= 0) return false; //проверка, находится ли точка внутри контура
}
return true;
}
TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy>
cbind(const CONT<TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy> >& cont,
const bool intersection)
{
size_t sz = cont.size();
if(sz == 0) {
return TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy>();
}
if(sz == 1) {
return *cont.begin();
}
std::vector<TDATE> ans_dates;
if(intersection) {
date_intersection(cont,std::inserter(ans_dates,ans_dates.begin()));
} else {
date_union(cont,std::inserter(ans_dates,ans_dates.begin()));
}
if(ans_dates.size() == 0) {
return TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy>();
}
TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy> ans(ans_dates.size(),cbind_total_cols(cont));
std::copy(ans_dates.begin(), ans_dates.end(), ans.getDates());
// FIXME: set colnames
// if union, then we must initialize the values of ans to NA
if(!intersection) {
for(TSDIM i = 0; i < ans.nrow() * ans.ncol(); i++) {
ans.getData()[i] = numeric_traits<TDATA>::NA();
}
}
typename CONT<TSeries<TDATE,TDATA,TSDIM,TSDATABACKEND,DatePolicy> >::const_iterator it = cont.begin();
// walk along iterator and map values into ans
TSDIM offset = 0;
while(it != cont.end()) {
cbind_map_values(ans, *it, offset);
offset += it->ncol();
++it;
}
return ans;
}
//empty
bool empty() {return container.size() == 0;}
void write_array(const CONT<T, K>& a, std::ostream& s)
{
write_size(a.size(), s);
write(a, s);
}