Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1644 1553 (N20W61 L=188)   591    112        0 (1598/1530-1700/1575) -32.0 W56-W65 N17-N22
      2 1272 1500 (N28W29 L=156)   382     89        2 (1232/1468-1302/1538)   8.1 W25-W32 N26-N30
      3 1645 1480 (N17W56 L=182)   399     59        1 (1621/1452-1666/1508)  76.0 W53-W57 N14-N18
      4  182 1075 (N25E65 L= 62)  1050     75        0 ( 159/1018- 210/1117) -71.6 E68-E60 N21-N27
      5  232 1028 (N22E58 L= 69)   948    134        1 ( 213/ 959- 254/1087) -82.9 E60-E54 N17-N25
      6 1868  966 (S19W60 L=186)  1130    132        3 (1842/ 903-1899/1026) -60.9 W56-W64 S23-S14
      7 1675  903 (S16W39 L=165)   361     68        1 (1649/ 880-1712/ 928)   0.0 W36-W42 S17-S14
      8  864  503 (S20E22 L=104)  1819    429        1 ( 728/ 372- 983/ 663)  27.4 E37-E10 S26-S12
      9 1414  594 (S28W15 L=141)   333    144        0 (1345/ 574-1480/ 616)  -5.7 W09-W19 S29-S27
     10  991  465 (S26E16 L=111)   396    131        1 ( 941/ 416-1027/ 504) -66.8 E18-E14 S29-S22
     11  576  470 (S17E44 L= 83)   280     49        0 ( 554/ 446- 602/ 498)  26.6 E46-E40 S18-S15
     12 1104  446 (S30E10 L=117)   247     75        0 (1066/ 430-1142/ 470)   0.0 E12-E06 S30-S29
     13  630  348 (S27E47 L= 80)  1582    284        1 ( 557/ 303- 728/ 440) -67.4 E51-E38 S30-S18
  

    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