Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  835 1824 (N62E27 L= 12)  1229    147        2 ( 777/1800- 912/1850)   6.3 E37-E14 N59-N65
      2 1160 1727 (N53W14 L= 52)   543     98        1 (1116/1706-1205/1744)  11.3 W08-W18 N50-N54
      3 1793 1414 (N26W67 L=106)   649     60        2 (1766/1388-1828/1432)  -5.7 W63-W74 N24-N27
      4 1176 1379 (N26W10 L= 49)   685    162        1 (1104/1319-1219/1418) -35.5 W05-W12 N22-N28
      5 1755  856 (S08W53 L= 92)   862    109        2 (1720/ 829-1815/ 884)  16.7 W49-W59 S09-S06
      6 1879  829 (S11W70 L=109)   551     77        1 (1857/ 799-1896/ 867)  31.0 W66-W72 S12-S08
      7 1375  759 (S13W23 L= 62)   660    226        2 (1342/ 696-1425/ 834) -66.0 W20-W26 S17-S08
      8 1642  672 (S20W45 L= 84)   569    179        0 (1592/ 618-1685/ 759)  60.9 W41-W47 S23-S13
      9  808  687 (S18E14 L= 25)   812    187        1 ( 719/ 657- 895/ 714)  14.0 E20-E08 S19-S15
     10  374  525 (S31E54 L=345)  2419    329        3 ( 295/ 438- 466/ 631) -47.5 E60-E46 S37-S22
     11 1383  567 (S26W26 L= 65)  1030    282        0 (1265/ 536-1519/ 600)   2.9 W16-W36 S28-S23
  

    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