Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1607 1673 (N51W34 L=194)   737    137        1 (1566/1655-1634/1702)  28.6 W28-W39 N47-N55
      2  867 1400 (N12E16 L=144)   853    335        2 ( 838/1290- 892/1525) -63.4 E20-E12 N05-N19
      3  229 1364 (N01E64 L= 96)   675    126        0 ( 199/1334- 256/1399) -15.3 E70-E59 N00-N03
      4 1551 1261 (N18W28 L=188)   516    174        1 (1527/1176-1571/1341)  90.0 W26-W28 N12-N23
      5 1137 1248 (N08W02 L=162)   537    161        1 (1099/1176-1166/1335)  90.0 W00-W03 N03-N13
      6  490 1311 (N01E39 L=122)   271     43        0 ( 462/1292- 510/1327)   0.0 E40-E36 N00-N02
      7  347 1196 (S08E47 L=113)   329     54        1 ( 320/1178- 367/1216) -14.0 E49-E45 S08-S06
      8 1610  780 (S08W42 L=202)   482     93        1 (1577/ 733-1642/ 824)  90.0 W39-W43 S11-S05
      9  976  726 (S25W03 L=163)   885    279        0 ( 812/ 668-1093/ 781)  38.7 E06-W09 S31-S19
     10 1439  618 (S21W35 L=195)   325     81        3 (1412/ 595-1478/ 639)  45.0 W33-W37 S22-S19
     11 1237  605 (S27W23 L=183)   257     31        0 (1212/ 590-1265/ 625)   0.0 W20-W24 S27-S25
     12 1461  330 (S37W52 L=212)  2551    460        5 (1278/ 297-1633/ 368)  27.5 W39-W64 S43-S30
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y