Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  760 1707 (N39E26 L=298)  1124    176        3 ( 701/1619- 801/1769) -45.0 E33-E20 N32-N44
      2  620 1708 (N39E38 L=286)   382     34        1 ( 601/1679- 640/1729) -63.4 E40-E36 N36-N40
      3  948 1572 (N30E09 L=315)   274    114        0 ( 935/1523- 962/1639) -76.0 E10-E07 N26-N34
      4  470 1576 (N28E45 L=279)   616     93        1 ( 422/1550- 517/1596)  -5.7 E50-E40 N26-N29
      5  323 1153 (N01E48 L=277)   557    120        0 ( 271/1120- 372/1193)  18.4 E52-E43 S01-N03
      6  509  854 (S17E33 L=291)   365    117        0 ( 493/ 794- 531/ 910)  66.8 E34-E30 S20-S13
      7  660  803 (S19E22 L=302)   230     97        0 ( 620/ 763- 702/ 835)  33.7 E25-E19 S21-S17
      8  321  752 (S23E50 L=274)   926    163        0 ( 269/ 689- 377/ 804)  30.3 E57-E43 S27-S20
      9  652  734 (S24E23 L=301)   216     77        0 ( 617/ 721- 692/ 754) -11.3 E25-E20 S24-S22
     10 1711  635 (S23W54 L= 18)   352     21        0 (1690/ 620-1735/ 651)   0.0 W51-W56 S23-S21
     11  673  504 (S39E24 L=300)   836    120        3 ( 626/ 463- 731/ 547)  16.7 E28-E18 S42-S35
     12 1754  499 (S30W69 L= 33)  1485    126        2 (1717/ 455-1793/ 541)  56.3 W65-W75 S33-S27
     13 1610  430 (S37W56 L= 20)   464     43        0 (1591/ 402-1630/ 455) -76.0 W53-W58 S38-S34
  

    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