Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1528 (N36W12 L=144)  1053    350        1 ( 803/1497-1125/1548) -26.6 W01-W21 N29-N40
      2 1563  900 (S16W28 L=160)  1002    366        0 (1432/ 798-1672/1032)  15.3 W17-W39 S18-S12
      3 1185  891 (S06W05 L=138)  1034    263        2 (1120/ 765-1246/ 996)  86.2 W03-W06 S14-N01
      4  244  891 (N17E57 L= 76)   411     75        1 ( 215/ 863- 267/ 918) -31.0 E60-E54 N14-N18
      5 1362  843 (S13W14 L=146)   248    110        0 (1335/ 803-1383/ 875)  78.7 W12-W15 S15-S10
      6 1590  801 (S23W28 L=160)   321    119        0 (1558/ 753-1626/ 836)  18.4 W24-W31 S25-S20
      7 1606  710 (S29W28 L=160)   298     62        0 (1580/ 690-1633/ 737) -59.0 W26-W30 S31-S26
      8 1404  581 (S31W11 L=143)   282     50        0 (1378/ 562-1432/ 601)   0.0 W08-W12 S31-S29
      9 1276  478 (S33E01 L=132)  1239    252        2 (1220/ 386-1328/ 570)  69.0 E07-W03 S39-S25
     10  572  450 (S16E46 L= 86)   346     63        1 ( 555/ 415- 587/ 486)  51.3 E47-E43 S18-S13
  

    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