Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1819 (N47W51 L=358)   546     79        0 (1368/1794-1435/1842)  12.1 W44-W58 N45-N48
      2 1128 1621 (N35W15 L=322)   306     60        0 (1090/1607-1161/1635) -21.8 W12-W17 N33-N36
      3  461 1437 (N31E35 L=272)   812    164        1 ( 410/1378- 522/1508) -39.8 E41-E29 N26-N36
      4 1682 1379 (N13W49 L=356)  1106    186        2 (1629/1317-1749/1445) -49.4 W45-W53 N09-N17
      5 1426 1380 (N16W29 L=337)   258     97        0 (1411/1337-1448/1423)  78.7 W28-W30 N13-N18
      6  723 1399 (N25E15 L=292)   228     75        0 ( 684/1378- 757/1420)  11.3 E17-E12 N23-N25
      7 1817 1098 (S04W55 L=  2)   311     44        1 (1800/1071-1830/1125) -90.0 W53-W56 S05-S02
      8  533 1053 (N06E30 L=278)  1894    419        7 ( 349/1000- 721/1103)  -6.6 E43-E17 N01-N09
      9 1881 1001 (S11W62 L=  9)  2070    155        8 (1840/ 929-1926/1071)  63.4 W56-W68 S14-S06
     10 1542  692 (S26W31 L=339)  4772    878        6 (1141/ 582-1913/ 838)   6.9 W02-W68 S30-S20
     11 1179  773 (S17W07 L=314)   216    120        0 (1133/ 742-1238/ 788)   8.1 W03-W10 S18-S16
     12  619  716 (S15E29 L=279)   579    216        1 ( 519/ 669- 698/ 740) -20.0 E35-E23 S18-S12
     13  939  631 (S23E10 L=298)   422    118        1 ( 906/ 596- 972/ 675)  45.0 E12-E07 S25-S20
     14  284  597 (S18E61 L=247)   816     93        2 ( 261/ 545- 304/ 645)  60.3 E64-E57 S21-S14
     15 1795  526 (S41W66 L= 13)   631     45        2 (1767/ 509-1825/ 543)   5.7 W60-W72 S41-S39
     16 1647  477 (S43W46 L=354)   537     79        1 (1610/ 459-1689/ 498)  -5.7 W41-W51 S43-S40
     17  427  414 (S31E57 L=250)  1618    208        1 ( 357/ 361- 521/ 440)  -5.7 E70-E47 S34-S28
     18  751  351 (S40E31 L=276)   439     90        0 ( 703/ 326- 799/ 374)   6.3 E36-E27 S41-S38
  

    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