Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  747 1705 (N39E27 L=298)  1107    175        3 ( 685/1612- 792/1764) -42.7 E34-E21 N31-N43
      2  611 1710 (N39E39 L=286)   434     62        1 ( 589/1669- 632/1736) -71.6 E41-E37 N35-N41
      3  462 1575 (N28E46 L=279)   528     88        1 ( 416/1550- 508/1594)  -5.7 E51-E41 N26-N29
      4  313 1152 (N01E48 L=277)   503     88        0 ( 274/1119- 350/1185)  18.4 E51-E45 S01-N02
      5  492  842 (S18E34 L=291)   333     95        0 ( 479/ 791- 511/ 890)  63.4 E35-E32 S20-S14
      6  650  808 (S19E23 L=302)   225     84        0 ( 614/ 774- 687/ 837)  38.7 E25-E20 S21-S17
      7  309  750 (S23E51 L=274)   968    153        1 ( 264/ 685- 363/ 791)  28.3 E58-E45 S27-S20
      8  649  741 (S23E23 L=302)   284    103        1 ( 599/ 724- 696/ 766)   0.0 E26-E19 S24-S21
      9  663  503 (S39E25 L=300)   937    118        3 ( 618/ 468- 719/ 544)  26.6 E29-E19 S41-S36
     10 1689  471 (S33W62 L= 27)  1948    243        2 (1585/ 395-1771/ 535)  45.0 W52-W74 S39-S27
     11 1259  414 (S42W24 L=349)   560    141        1 (1225/ 388-1320/ 469) -21.8 W19-W29 S43-S37
  

    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