Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  983 1753 (N48W11 L=164)   352     31        0 ( 965/1727-1004/1773)  36.9 W09-W13 N46-N50
      2  865 1571 (N36E02 L=151)   264     45        0 ( 834/1558- 897/1587) -26.6 E04-W00 N35-N37
      3  752 1547 (N36E11 L=142)   383     92        0 ( 717/1508- 787/1582)  38.7 E13-E08 N34-N38
      4 1380  883 (S13W20 L=173)   219     43        0 (1363/ 853-1397/ 910)  56.3 W18-W20 S14-S11
      5 1377  753 (S21W18 L=171)   669    171        0 (1304/ 691-1440/ 818)  33.7 W13-W23 S23-S17
      6 1204  697 (S22W07 L=160)   407    124        0 (1142/ 670-1267/ 729) -35.5 W02-W10 S24-S19
      7  415  646 (S14E47 L=106)   669     77        1 ( 377/ 609- 443/ 688)  38.7 E50-E43 S16-S11
      8  845  590 (S23E18 L=135)   509    179        0 ( 815/ 499- 885/ 662)  45.0 E22-E13 S28-S19
      9  480  635 (S15E42 L=111)   243     43        0 ( 461/ 609- 494/ 657) -76.0 E42-E40 S17-S13
     10  213  561 (S16E76 L= 77)  2261    106        1 ( 181/ 504- 238/ 607)  80.5 E82-E70 S19-S13
     11  545  388 (S32E49 L=104)  1898    307        0 ( 468/ 335- 626/ 469) -68.2 E55-E39 S35-S25
  

    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