Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1220 1784 (N47E15 L=107)   447     69        2 (1195/1755-1244/1812) -39.8 E18-E12 N44-N49
      2  655 1722 (N24E48 L= 74)   583     93        2 ( 633/1683- 671/1770) -49.4 E53-E44 N21-N28
      3  524 1660 (N17E54 L= 67)  1193    138        2 ( 473/1621- 573/1710)  74.1 E58-E50 N14-N21
      4 1542 1684 (N53W17 L=139)   434     57        0 (1512/1667-1573/1702)  45.0 W14-W21 N50-N57
      5  619 1514 (N10E38 L= 83)   345     56        1 ( 590/1482- 640/1542)  90.0 E39-E37 N07-N12
      6 1342 1516 (N31W05 L=126)   245     37        0 (1323/1490-1365/1536)  76.0 W03-W05 N28-N32
      7 1826 1207 (N29W50 L=171)   426     99        0 (1785/1170-1863/1238)  -5.7 W45-W55 N28-N30
      8 1613  481 (S17W58 L=179)   368     54        0 (1584/ 452-1641/ 506)  76.0 W56-W59 S19-S15
      9 1191  433 (S35W32 L=153)   625    144        1 (1159/ 366-1222/ 501) -28.6 W26-W38 S37-S30
     10 1362  398 (S31W45 L=167)   308     70        0 (1334/ 374-1394/ 419)  76.0 W43-W46 S32-S28
  

    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