Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  713 1585 (N21E32 L=205)  1055    206        1 ( 691/1493- 739/1679) -69.0 E36-E28 N14-N27
      2 1355 1557 (N33W10 L=247)  1255    313        3 (1288/1438-1418/1671) -49.4 W03-W15 N25-N39
      3 1868 1103 (N18W58 L=295)   487     36        0 (1838/1082-1895/1118)  15.9 W54-W61 N16-N19
      4  204  944 (S25E60 L=178)  1683    253        0 ( 164/ 862- 288/1052) -39.8 E65-E50 S28-S18
      5  336  768 (S34E46 L=192)   310     85        0 ( 295/ 754- 370/ 785)   0.0 E50-E42 S34-S32
      6  348  485 (S54E56 L=181)  1386    118        1 ( 294/ 442- 408/ 520)   8.6 E77-E44 S57-S51
      7 1196  345 (S44W34 L=271)   657    141        0 (1164/ 292-1242/ 388) -14.0 W29-W41 S46-S40
      8 1524  337 (S34W61 L=299)  1296    207        3 (1459/ 309-1620/ 377)  36.0 W55-W68 S37-S28
  

    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