Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1552 1573 (N20W52 L=  5)   528    202        0 (1497/1541-1579/1632)  14.0 W45-W57 N18-N22
      2  180 1131 (N32E57 L=256)  1823    361        0 ( 120/1027- 251/1220) -49.4 E63-E51 N24-N38
      3  821  983 (N09E12 L=301)   449     93        1 ( 763/ 963- 877/1002) -33.7 E15-E09 N06-N10
      4  946  960 (N04E06 L=308)   329     89        0 ( 907/ 932-1000/ 995)   9.5 E08-E02 N03-N05
      5  526  940 (N13E32 L=282)   838    281        2 ( 419/ 886- 654/ 975)  -3.2 E40-E22 N11-N15
      6 1323  892 (S10W13 L=327)   276     74        1 (1278/ 864-1355/ 924)   0.0 W09-W15 S10-S08
      7  310  872 (N14E50 L=264)   477    165        1 ( 248/ 847- 356/ 916) -32.0 E54-E45 N11-N17
      8 1527  844 (S18W24 L=338)   249     63        2 (1509/ 813-1546/ 875)  56.3 W22-W25 S19-S16
      9  200  791 (N11E64 L=250)   536     69        0 ( 182/ 759- 225/ 836) -78.7 E65-E61 N09-N14
     10 1026  519 (S23E14 L=299)  1467    365        6 ( 994/ 382-1079/ 652)  49.1 E20-E07 S31-S15
     11  373  490 (S09E60 L=253)  3552    524       12 ( 280/ 426- 445/ 565) -16.4 E73-E53 S14-S04
     12  491  455 (S14E51 L=262)   374    123        2 ( 476/ 413- 515/ 500)  39.8 E54-E48 S16-S11
  

    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