Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  647 1648 (N22E43 L=274)  4331    827       10 ( 594/1438- 694/1874) -45.0 E62-E32 N09-N39
      2  467 1588 (N14E54 L=263)   873    273        1 ( 443/1493- 486/1684) -53.1 E60-E49 N08-N20
      3  364 1540 (N09E61 L=256)  1852    347        5 ( 339/1454- 410/1637) -36.3 E68-E52 N04-N15
      4 1695 1526 (N48W37 L=354)   701    106        2 (1660/1484-1721/1551)  33.7 W32-W41 N44-N51
      5 1593 1017 (N11W34 L=351)   842    232        1 (1536/ 898-1641/1118) -36.9 W27-W40 N05-N14
      6  476 1076 (S15E35 L=283)   520    123        1 ( 423/1053- 533/1095)  26.6 E38-E30 S16-S12
      7  164 1069 (S22E63 L=255)   950    253        1 ( 112/1056- 216/1087)   8.7 E69-E56 S23-S20
      8 1690  593 (S09W57 L= 15)  1135    204        5 (1656/ 522-1724/ 649) -63.4 W52-W61 S13-S05
      9 1592  579 (S13W49 L=  7)   804    171        2 (1515/ 538-1656/ 642)  59.0 W45-W52 S17-S07
  

    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