Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1639 1569 (N17W59 L=227)  2659    592        5 (1459/1531-1782/1596) -21.5 W45-W78 N11-N24
      2 1478 1341 (N11W36 L=203)   259    101        0 (1464/1287-1492/1390)  63.4 W34-W37 N07-N13
      3 1883 1195 (S11W65 L=232)  1154    201        3 (1855/1132-1904/1255)  54.5 W59-W69 S14-S06
      4  564 1169 (N26E24 L=143)   675    247        1 ( 478/1062- 634/1238) -63.4 E29-E21 N18-N32
      5 1820  952 (S22W50 L=218)  2752    408        5 (1766/ 798-1855/1111)  78.1 W44-W52 S31-S11
      6 1666  797 (S25W34 L=201)   768    209        1 (1646/ 689-1704/ 889)  63.4 W30-W38 S32-S20
      7 1037  610 (S17E11 L=157)   629    219        1 ( 961/ 516-1092/ 699)  28.3 E17-E04 S20-S13
      8  550  606 (S05E40 L=127)  1000    170        2 ( 502/ 520- 605/ 680) -90.0 E41-E38 S11-N00
      9 1179  441 (S32E08 L=160)   390     73        3 (1153/ 410-1204/ 474)  31.0 E10-E05 S33-S29
     10  990  337 (S34E24 L=143)  1064    261        3 ( 918/ 260-1067/ 411)  19.4 E33-E16 S37-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