Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  305 1538 (N17E67 L=349)   721     86        1 ( 282/1507- 325/1574)  90.0 E69-E64 N14-N19
      2 1197 1531 (N34W01 L= 57)   500    136        0 (1124/1499-1263/1558)  -5.7 E04-W06 N32-N34
      3 1782  982 (N11W53 L=108)   340     89        1 (1751/ 966-1808/1000)   0.0 W49-W55 N09-N12
      4  262  890 (S23E51 L=  5)   339     79        0 ( 236/ 877- 295/ 904)   0.0 E53-E48 S24-S22
      5 1201  801 (S11W15 L= 71)   416    100        1 (1185/ 746-1219/ 853) -56.3 W12-W16 S14-S08
      6 1014  779 (S16W04 L= 60)   203    104        1 ( 975/ 764-1071/ 803)   9.5 W01-W07 S16-S14
      7  898  665 (S25E01 L= 55)   256    128        0 ( 833/ 650- 964/ 676)  14.0 E05-W03 S27-S23
      8 1128  632 (S23W14 L= 70)   210     47        0 (1105/ 605-1144/ 665)  90.0 W13-W15 S24-S20
      9 1738  512 (S17W66 L=122)  1609    219        3 (1687/ 444-1789/ 563)  60.9 W61-W71 S22-S13
     10 1370  468 (S28W36 L= 92)  2829    588        2 (1139/ 376-1645/ 530)  28.6 W23-W58 S39-S21
  

    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