Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1010 1790 (N41E25 L=152)   329     44        1 ( 980/1774-1037/1809)  45.0 E26-E22 N39-N43
      2  626 1591 (N18E40 L=137)   488    149        1 ( 571/1559- 677/1628)  45.0 E43-E36 N15-N20
      3  729 1569 (N18E32 L=145)   241     65        0 ( 711/1534- 746/1603) -53.1 E34-E29 N16-N20
      4 1754 1369 (N33W43 L=220)   415     39        2 (1731/1342-1777/1390) -45.0 W40-W46 N31-N34
      5 1229 1206 (N08W07 L=184)   828    256        0 (1132/1167-1356/1233)  28.3 W02-W15 N03-N10
      6 1741 1140 (N18W43 L=220)   478    132        3 (1695/1094-1786/1172)  -5.7 W38-W48 N16-N20
      7  630 1014 (S16E23 L=153)   313    118        1 ( 605/ 958- 655/1069)  90.0 E24-E22 S19-S12
      8  233  953 (S27E56 L=120)  1588    232        1 ( 187/ 863- 281/1055)  52.4 E62-E51 S32-S19
      9 1756  980 (N09W48 L=224)   400     61        2 (1734/ 958-1785/1010) -14.0 W45-W50 N07-N10
     10 1476  579 (S20W41 L=217)   285     52        1 (1447/ 562-1508/ 593)  11.3 W38-W43 S21-S19
  

    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