Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1446 1600 (N25W45 L=212)  3989    702        8 (1133/1539-1704/1635) -25.4 W25-W65 N15-N34
      2  990 1431 (N28W09 L=177)   373    101        1 ( 951/1391-1017/1473) -10.7 W07-W10 N25-N31
      3 1929 1155 (S14W69 L=237)  2804    241        2 (1890/1091-1964/1233) -84.3 W65-W74 S19-S09
      4 1749  795 (S29W41 L=209)  2811    484        3 (1601/ 619-1815/1032)  39.6 W27-W50 S37-S16
      5  717  914 (N06E20 L=147)   487    161        2 ( 660/ 866- 805/ 944)   0.0 E24-E14 N03-N07
      6 1624  807 (S24W31 L=199)   663    178        2 (1567/ 721-1677/ 891)  32.5 W25-W36 S27-S20
      7  934  740 (S09E12 L=155)   208     72        0 ( 895/ 727- 974/ 753) -11.3 E14-E09 S09-S07
      8 1310  608 (S27W06 L=173)   433    170        0 (1261/ 542-1345/ 699)  37.9 W01-W10 S30-S22
      9  642  630 (S08E33 L=135)   293     80        0 ( 627/ 573- 657/ 662)  59.0 E34-E31 S11-S06
     10 1045  594 (S20E10 L=157)   250     82        0 (1005/ 558-1074/ 637)  26.6 E13-E07 S21-S18
     11  884  433 (S26E26 L=142)   812    278        1 ( 828/ 339- 932/ 538)  33.7 E33-E18 S31-S20
     12  392  416 (S16E63 L=105)  2136    392        8 ( 352/ 358- 429/ 494)  45.0 E69-E54 S19-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