Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  658 1611 (N51E09 L=264)   338     62        1 ( 638/1574- 673/1642)  59.0 E10-E06 N47-N53
      2 1328 1598 (N28W38 L=311)   355     36        0 (1300/1580-1359/1615) -33.7 W35-W40 N26-N30
      3  547 1438 (N41E22 L=251)   298     85        0 ( 519/1403- 580/1474) -56.3 E24-E19 N38-N44
      4  426 1155 (N26E35 L=238)   266     16        0 ( 407/1135- 445/1175) -45.0 E36-E33 N25-N27
      5  493 1079 (N20E31 L=242)   388    100        0 ( 452/1049- 547/1100)  -8.1 E34-E27 N19-N21
      6 1865  788 (S34W58 L=331)  1373    105        4 (1839/ 732-1891/ 853)  90.0 W54-W62 S37-S29
      7 1704  645 (S37W38 L=311)   404     51        1 (1682/ 621-1730/ 674)  56.3 W35-W40 S39-S35
      8 1533  624 (S33W22 L=295)   251     64        0 (1497/ 608-1562/ 646)   0.0 W18-W24 S33-S31
      9 1204  547 (S28E02 L=271)   391    155        0 (1147/ 500-1260/ 588)  20.6 E06-W02 S29-S25
     10 1341  226 (S60E05 L=269)   634     59        2 (1309/ 205-1377/ 250)  -6.3 E09-W00 S62-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