Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1125 1566 (N33W23 L=147)   594    121        1 (1065/1532-1174/1591) -33.7 W19-W25 N29-N35
      2 1347 1499 (N22W34 L=158)   385    127        1 (1297/1459-1400/1533)   5.7 W28-W38 N21-N23
      3  598 1187 (N24E22 L=103)   228     77        0 ( 582/1149- 615/1229) -90.0 E22-E20 N21-N26
      4  206 1125 (N29E58 L= 66)  1100    111        3 ( 167/1085- 248/1178)  45.0 E63-E52 N26-N32
      5  818 1052 (N11E11 L=113)   193     90        0 ( 804/1006- 834/1085)  53.1 E12-E09 N08-N13
      6 1875  997 (S21W59 L=184)   419     70        1 (1860/ 961-1893/1032) -90.0 W57-W62 S22-S18
      7 1832  948 (S22W53 L=177)   746     95        0 (1810/ 908-1857/ 997) -80.5 W50-W56 S25-S19
      8 1679  650 (S36W35 L=160)  2517    641        4 (1455/ 503-1890/ 857)   8.5 W14-W61 S40-S30
      9  731  650 (S09E27 L= 98)   358     77        1 ( 713/ 599- 758/ 680)  45.0 E29-E24 S12-S07
     10 1147  403 (S36E11 L=114)   243     65        1 (1113/ 382-1174/ 427)  21.8 E13-E08 S36-S34
     11  804  387 (S27E33 L= 91)   796    115        1 ( 746/ 367- 868/ 415) -35.5 E36-E28 S30-S24
     12  593  338 (S25E52 L= 73)   791    132        3 ( 533/ 299- 663/ 366)  -4.4 E57-E44 S27-S23
     13  995  142 (S53E43 L= 81)   491     42        1 ( 976/ 122-1020/ 159)  16.7 E48-E38 S55-S51
  

    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