Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1682 (N27E33 L= 89)   373     78        0 ( 820/1637- 855/1720) -39.8 E36-E30 N24-N29
      2  711 1475 (N11E31 L= 91)  1344    406        0 ( 508/1412- 826/1510)  18.4 E44-E23 N05-N14
      3 1045 1380 (N14E09 L=114)   331    120        0 (1021/1328-1077/1421) -71.6 E09-E06 N11-N17
      4 1707 1307 (N31W37 L=159)   327     70        0 (1682/1273-1734/1343) -21.8 W34-W39 N30-N33
      5 1225 1159 (N07W07 L=130)   503    161        0 (1211/1080-1241/1241) -60.9 W04-W09 N02-N11
      6 1397 1091 (N08W19 L=141)   294     73        2 (1379/1050-1413/1132) -68.2 W17-W20 N05-N10
      7  783  791 (S26E07 L=115)   286     81        0 ( 743/ 772- 826/ 813)  45.0 E09-E05 S27-S23
      8  922  526 (S37W10 L=132)   303    115        0 ( 875/ 486- 967/ 563)  69.4 W08-W11 S41-S33
      9 1563  497 (S18W53 L=175)   331     30        1 (1541/ 478-1587/ 513)  33.7 W50-W54 S18-S16
     10 1214  454 (S32W32 L=154)   266     59        1 (1185/ 437-1239/ 478)  33.7 W29-W33 S33-S30
     11 1594  401 (S21W64 L=186)  1039    171        0 (1501/ 368-1669/ 437)  39.3 W57-W69 S26-S16
  

    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