Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  498 1274 (N25E30 L=207)   542    146        2 ( 442/1219- 544/1317) -54.5 E34-E26 N21-N28
      2  177 1210 (N27E62 L=174)   716     53        1 ( 151/1186- 207/1232) -26.6 E67-E58 N24-N28
      3  230 1138 (N21E55 L=182)   391     47        0 ( 202/1119- 263/1158)   0.0 E58-E51 N20-N22
      4  239 1082 (N18E54 L=183)   397     57        0 ( 214/1061- 270/1099)   9.5 E56-E50 N16-N19
      5 1900  958 (S19W65 L=301)   467     61        0 (1885/ 927-1917/ 985)  45.0 W62-W67 S20-S17
      6 1760  705 (S32W49 L=286)  2077    268        2 (1661/ 629-1836/ 807)  22.4 W38-W57 S36-S26
      7 1516  708 (S26W26 L=263)   447    125        0 (1471/ 658-1563/ 749)  26.6 W22-W30 S28-S24
      8  928  469 (S31E18 L=219)   388    118        0 ( 879/ 445- 991/ 489)   0.0 E21-E12 S31-S29
      9 1621  363 (S56W46 L=283)   852     84        0 (1593/ 342-1658/ 404) -16.4 W40-W57 S59-S51
     10  611  375 (S31E45 L=192)   586     82        0 ( 567/ 354- 658/ 399) -23.2 E47-E40 S32-S28
     11 1011  131 (S62E36 L=200)   667     78        0 ( 998/  89-1026/ 171)  29.7 E52-E29 S68-S56
  

    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