Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  424 1392 (N18E47 L=242)   897    238        6 ( 399/1301- 447/1485) -47.7 E52-E42 N12-N23
      2  710 1263 (N13E22 L=267)   286    108        2 ( 671/1241- 756/1282)  -9.5 E25-E19 N11-N14
      3  827 1032 (N00E12 L=277)   274     57        1 ( 801/1013- 854/1052)  18.4 E13-E10 S01-N01
      4  855  862 (S10E09 L=280)   366     61        2 ( 829/ 836- 882/ 887)  56.3 E10-E07 S11-S08
      5  893  750 (S17E06 L=283)   518    174        2 ( 810/ 729- 958/ 790)   0.0 E11-E01 S17-S14
      6 1554  718 (S13W39 L=328)   796    177        0 (1523/ 640-1604/ 788)  78.7 W37-W41 S17-S07
      7  218  739 (S25E64 L=225)   577     78        1 ( 178/ 727- 248/ 764)  13.0 E71-E58 S26-S23
      8 1320  694 (S16W23 L=312)   437    172        0 (1259/ 656-1371/ 763)  54.5 W18-W24 S19-S11
      9 1706  648 (S16W55 L=344)   400     65        0 (1681/ 628-1725/ 681)  63.4 W52-W56 S17-S13
     10 1569  556 (S23W45 L=335)   352    110        1 (1542/ 510-1602/ 598)  80.5 W43-W46 S26-S19
     11  741  543 (S32E15 L=274)   239     86        1 ( 689/ 529- 776/ 559)  -8.1 E19-E12 S33-S31
     12 1493  360 (S38W50 L=340)  1002    121        1 (1453/ 326-1537/ 411)  82.9 W45-W54 S41-S33
  

    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