Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  435 1650 (N23E64 L= 41)   483     71        0 ( 401/1627- 471/1664)  -4.8 E69-E57 N22-N24
      2 1085 1528 (N30E08 L= 96)   300     96        0 (1041/1500-1128/1547)   0.0 E11-E04 N28-N30
      3  889 1388 (N16E16 L= 88)   578    163        3 ( 833/1329- 929/1423) -29.7 E19-E12 N13-N18
      4 1121 1385 (N21E02 L=102)   459    103        2 (1083/1355-1154/1419)  26.6 E04-W00 N19-N23
      5 1147 1313 (N17W01 L=106)   686    181        0 (1073/1283-1227/1340)  32.0 E03-W05 N14-N20
      6  186 1065 (S17E63 L= 41)  1152    137        3 ( 158/1002- 214/1131)  69.4 E67-E60 S20-S12
      7 1772 1013 (N14W51 L=156)   463    118        0 (1758/ 955-1789/1067) -56.3 W49-W53 N10-N17
      8  842  692 (S26E04 L=101)  1093    293        4 ( 722/ 642- 946/ 782)  52.4 E11-W01 S31-S18
      9 1622  423 (S24W64 L=168)  1113     71        2 (1587/ 397-1652/ 455)  76.0 W59-W69 S25-S21
     10  699  406 (S50E09 L= 96)   430     86        1 ( 668/ 385- 745/ 435)  35.5 E11-E04 S52-S47
     11  973  339 (S47W17 L=122)  1439    247        2 ( 859/ 319-1098/ 360)  22.8 W07-W26 S51-S41
  

    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