area=0;
for col=1:464
    temp=find(ds(:,col)==1);
    if temp
        area= area + (temp(end)-temp(1)+1);
    end
end
area;