Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1542 1402 (N13W43 L=342)   395     89       -1 (1511/1383-1569/1421) -45.0 W40-W45 N11-N14
      2 1477 1061 (S03W27 L=326)   532    158        3 (1405/1011-1533/1092)   6.3 W21-W30 S05-S01
      3  986 1001 (N07E03 L=296)   225    121        0 ( 942/ 955-1042/1043)  15.9 E06-W01 N05-N07
      4  615  987 (N15E25 L=274)   343     79        0 ( 577/ 959- 652/1010)   0.0 E27-E22 N13-N15
      5  802  985 (N10E14 L=286)   565    194        1 ( 717/ 965- 902/1006) -35.0 E18-E08 N06-N13
      6  175  824 (N13E67 L=232)   866    168        4 ( 140/ 799- 207/ 862)  16.7 E72-E62 N11-N15
      7  244  764 (N09E61 L=238)  1174    267        2 ( 183/ 740- 321/ 793) -10.6 E68-E52 N07-N10
      8 1307  739 (S17W09 L=308)   626    118        3 (1261/ 705-1359/ 784) -18.4 W05-W11 S19-S13
      9 1174  550 (S24E04 L=295)   463    183        0 (1136/ 477-1198/ 646)  51.3 E08-W00 S28-S18
     10  394  477 (S11E59 L=240)  2326    464        7 ( 340/ 405- 441/ 552) -77.5 E65-E52 S16-S05
     11  481  495 (S11E50 L=249)   637    168        5 ( 445/ 438- 509/ 524)  35.5 E54-E47 S14-S08
  

    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