Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  757 1618 (N21E33 L=310)   269     21        0 ( 735/1600- 780/1635)  18.4 E34-E31 N19-N22
      2  517 1536 (N11E46 L=297)   938    220        0 ( 470/1457- 550/1607) -30.3 E52-E40 N05-N14
      3  266 1109 (S19E52 L=290)  1442    445        0 ( 147/1069- 391/1164)   4.8 E66-E42 S21-S15
      4  906 1040 (S09E07 L=336)   381    112        0 ( 874/ 988- 963/1074)  63.4 E08-E04 S12-S06
      5  425 1043 (S20E38 L=305)   698    165        1 ( 355/1026- 500/1062)  14.0 E43-E31 S21-S18
      6 1234  982 (S04W13 L=355)   229     95        0 (1192/ 954-1284/1007)  31.0 W10-W15 S06-S01
      7 1897  896 (N12W64 L= 47)  2936    494        5 (1852/ 820-1939/ 977) -73.3 W57-W72 N07-N17
      8  303  818 (S36E45 L=297)   517    182        0 ( 243/ 795- 362/ 838)   9.5 E51-E39 S37-S34
      9 1570  602 (S15W46 L= 29)   931    181        0 (1549/ 514-1597/ 684) -55.0 W43-W50 S19-S09
     10  325  592 (S51E46 L=297)   347     86        0 ( 300/ 571- 351/ 612)  23.2 E48-E41 S52-S48
     11 1267  487 (S30W32 L= 14)  2480    618        1 (1002/ 435-1547/ 547)  17.5 W13-W51 S35-S22
     12 1655  477 (S18W60 L= 42)   426     71        2 (1635/ 452-1676/ 496)  71.6 W57-W62 S19-S16
  

    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