Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1426 1633 (N40W18 L=235)  1510    330        4 (1294/1544-1573/1725) -16.3 W06-W30 N36-N45
      2  513 1622 (N23E48 L=169)   808    138        4 ( 475/1576- 554/1674) -81.9 E50-E44 N19-N26
      3 1219 1507 (N26W04 L=221)   227    124        1 (1183/1478-1270/1546) -18.4 W01-W07 N25-N28
      4 1490 1386 (N24W24 L=241)   281    113        1 (1443/1349-1535/1428) -23.2 W19-W27 N22-N25
      5  723  838 (S22E16 L=202)   321    130        1 ( 683/ 802- 759/ 904)  66.8 E18-E13 S24-S17
      6  233  822 (S29E58 L=160)  1301    159        4 ( 159/ 780- 309/ 857)   7.8 E70-E48 S31-S26
      7 1507  734 (S14W36 L=253)   401    142        0 (1452/ 688-1546/ 784)  66.8 W33-W37 S17-S10
      8 1748  715 (S10W56 L=273)   752    132        2 (1695/ 682-1809/ 749)  26.6 W51-W62 S12-S07
      9  921  705 (S27E01 L=216)   309    124        2 ( 867/ 680- 979/ 730)  18.4 E04-W02 S28-S24
     10 1720  641 (S14W57 L=274)   433     44        1 (1700/ 618-1745/ 665)  76.0 W54-W59 S16-S12
     11 1363  558 (S27W32 L=250)  2440    552        6 (1134/ 500-1611/ 613)  25.8 W17-W48 S35-S18
     12 1760  560 (S18W66 L=284)   728     87        1 (1718/ 522-1790/ 582)  18.4 W60-W70 S19-S15
  

    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