Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1681 (N25E52 L=107)  1468    317        2 ( 502/1620- 596/1757) -21.8 E61-E46 N20-N29
      2  379 1588 (N15E62 L= 97)  3038    489       10 ( 305/1479- 442/1666) -53.1 E70-E55 N08-N20
      3  563 1398 (N07E38 L=122)   420    125        1 ( 521/1345- 589/1444) -71.6 E40-E36 N04-N10
      4  929 1170 (N03E09 L=150)   314     95        1 ( 893/1118- 960/1221)  80.5 E10-E07 N00-N06
      5 1393  957 (N02W23 L=182)   286    112        1 (1343/ 905-1432/1008)  74.1 W21-W24 S02-N05
      6 1168  982 (S02W09 L=169)   190     83        1 (1127/ 970-1219/1001)   9.5 W06-W12 S03-N00
      7 1396  645 (S15W33 L=192)   267     89        1 (1373/ 602-1415/ 676) -36.9 W30-W34 S17-S13
      8 1602  559 (S15W51 L=210)   661     99        0 (1558/ 521-1641/ 609)  74.1 W49-W52 S17-S10
      9 1689  560 (S12W59 L=218)   482     91        2 (1657/ 520-1724/ 594)  90.0 W57-W61 S14-S09
     10 1055  423 (S38W21 L=180)  1172    318        1 ( 919/ 399-1236/ 455)  26.6 W10-W31 S42-S32
     11 1435  308 (S34W54 L=213)   467    130        0 (1408/ 290-1471/ 340)  -8.1 W49-W56 S35-S32
  

    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