Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  714 1655 (N27E38 L=244)   928    316        2 ( 626/1573- 784/1742) -20.2 E49-E30 N23-N31
      2  427 1630 (N19E60 L=222)   485    134        1 ( 412/1592- 439/1673) -50.2 E62-E57 N16-N22
      3  496 1590 (N18E51 L=231)   464     79        1 ( 472/1551- 527/1626)  90.0 E52-E49 N15-N20
      4  372 1469 (N08E56 L=227)  1062    116        3 ( 340/1411- 402/1518) -60.3 E59-E51 N04-N11
      5  341 1122 (S13E47 L=236)   406     61        2 ( 318/1089- 363/1153) -76.0 E48-E44 S14-S10
      6  279 1064 (S17E51 L=232)   673    173        3 ( 241/1006- 312/1121)  76.0 E53-E48 S21-S13
      7  208 1037 (S20E58 L=225)   908    192        0 ( 168/ 950- 232/1109)  79.7 E60-E55 S26-S15
      8  832  913 (S14E09 L=274)   315     79        1 ( 796/ 883- 866/ 938)  63.4 E10-E06 S16-S12
      9  167  769 (S38E61 L=222)   695    165        2 ( 148/ 741- 198/ 800)  63.4 E63-E56 S40-S36
     10  787  687 (S29E06 L=276)   312     70        0 ( 751/ 656- 818/ 712)  59.0 E08-E04 S31-S26
     11  627  571 (S41E16 L=267)   292     36        2 ( 605/ 548- 653/ 594)  26.6 E17-E13 S42-S38
  

    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