Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1503 1589 (N26W50 L= 84)  1221    278        3 (1391/1513-1591/1659) -65.2 W44-W54 N20-N33
      2 1789 1352 (N05W62 L= 96)  2170    343        3 (1754/1195-1822/1489)  51.3 W56-W74 S04-N11
      3 1343 1316 (N16W26 L= 60)   254    104        0 (1295/1295-1395/1334)  -8.1 W22-W29 N14-N16
      4  953 1163 (N17E01 L= 33)   615    174        2 ( 869/1127-1032/1204) -45.0 E05-W02 N12-N20
      5  674 1148 (N22E19 L= 16)   420    152        2 ( 602/1122- 741/1181) -29.1 E23-E14 N19-N24
      6 1633  889 (S16W34 L= 69)  1458    461        1 (1493/ 783-1731/1064)  25.5 W23-W45 S19-S07
      7  926 1029 (N10E06 L= 29)   408    158        1 ( 892/ 986- 992/1058) -26.6 E07-E01 N07-N11
      8  539  859 (N08E33 L=  1)   538    219        0 ( 503/ 762- 562/ 960) -90.0 E34-E31 N01-N13
      9 1718  573 (S40W41 L= 76)  3821    485        3 (1532/ 433-1891/ 701)   7.6 W22-W68 S44-S36
     10  932  657 (S12E14 L= 21)   306     72        0 ( 890/ 639- 971/ 678) -45.0 E15-E11 S13-S09
     11  473  502 (S13E50 L=345)  1609    221        2 ( 451/ 409- 516/ 622)  57.3 E55-E44 S19-S05
     12 1316  321 (S46W01 L= 35)   603    109        2 (1256/ 296-1370/ 348)   5.7 E04-W06 S47-S44
  

    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