Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1220 1569 (N30W25 L=138)   266     73        0 (1186/1541-1253/1597) -63.4 W23-W26 N27-N32
      2 1322 1467 (N22W28 L=141)   342     63        1 (1292/1436-1355/1487) -45.0 W25-W29 N19-N23
      3  438 1357 (N32E35 L= 79)   309     95        0 ( 398/1330- 478/1373)   8.1 E38-E31 N30-N33
      4 1893 1019 (S15W64 L=178)  1160    154        1 (1844/ 960-1928/1079)  18.4 W57-W70 S17-S11
      5 1344  953 (S08W17 L=131)   324    102        1 (1324/ 910-1362/1007) -90.0 W15-W18 S10-S04
      6  179  820 (N04E66 L= 48)  1119    126        1 ( 122/ 789- 219/ 867)  10.0 E76-E59 N02-N06
      7  864  747 (S12E14 L= 99)   939    256        1 ( 805/ 627- 930/ 846)  39.8 E20-E08 S17-S06
      8 1830  761 (S30W57 L=170)   460     55        2 (1805/ 741-1853/ 790)  21.8 W53-W60 S31-S28
      9  338  697 (S06E52 L= 62)   601     92        1 ( 302/ 657- 375/ 734)  26.6 E57-E48 S07-S03
     10 1748  647 (S36W49 L=162)  1570    241        0 (1678/ 580-1848/ 729) -25.8 W39-W71 S43-S28
     11  362  570 (S14E56 L= 58)  1208    128        1 ( 323/ 526- 402/ 637)  26.6 E61-E50 S16-S10
  

    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