Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1735 (N36E60 L=247)   836    106        2 ( 519/1690- 578/1782) -18.4 E70-E52 N32-N39
      2  710 1606 (N30E35 L=272)   315    118        0 ( 673/1560- 751/1640) -21.8 E39-E29 N27-N31
      3 1139 1540 (N34E02 L=304)   294     59        0 (1119/1512-1162/1563) -76.0 E04-W00 N31-N35
      4 1211 1490 (N32W04 L=310)   368    115        0 (1171/1438-1258/1521) -29.7 W00-W07 N28-N33
      5  872 1426 (N21E17 L=290)   359    140        1 ( 823/1382- 921/1477) -24.0 E21-E12 N18-N23
      6 1683  702 (S08W51 L=358)   347     68        0 (1649/ 687-1713/ 717)  31.0 W48-W54 S09-S06
      7  752  624 (S30E11 L=296)   785    259        1 ( 646/ 586- 876/ 653)  23.6 E18-E02 S34-S26
      8 1217  473 (S31W25 L=332)  1703    292        2 (1147/ 373-1291/ 571) -30.5 W17-W34 S37-S25
      9 1430  382 (S33W47 L=353)   887    112        2 (1379/ 354-1481/ 414)   5.2 W41-W52 S35-S31
  

    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