Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1797 (N73E24 L= 74)   535     53        0 ( 624/1768- 665/1821)  56.3 E29-E19 N69-N76
      2  821 1612 (N49E00 L= 98)   916    145        2 ( 755/1560- 888/1645) -26.6 E05-W05 N44-N51
      3  832 1326 (N29E06 L= 92)   414    144        0 ( 800/1276- 890/1379)  26.6 E08-W00 N26-N31
      4  392 1162 (N26E41 L= 57)  1448    375        2 ( 284/1080- 504/1221) -29.1 E52-E32 N20-N30
      5  219 1063 (N22E60 L= 38)   902    133        1 ( 185/ 997- 255/1124)  41.2 E64-E55 N18-N25
      6  978  854 (S02E06 L= 92)   204     53        0 ( 962/ 831-1002/ 883) -71.6 E07-E04 S03-N00
      7 1232  796 (S11W08 L=106)   208     72        0 (1207/ 763-1258/ 817)  36.9 W05-W09 S12-S09
      8 1041  471 (S27E11 L= 87)   305     52        1 (1016/ 444-1069/ 502)  36.9 E13-E08 S29-S25
      9  506  427 (S21E52 L= 47)   623    107        1 ( 475/ 394- 555/ 454) -68.2 E55-E48 S23-S18
     10  631  368 (S27E44 L= 54)  1299    213        2 ( 580/ 279- 702/ 446)  78.7 E46-E41 S35-S20
  

    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