Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1409 1779 (N36W60 L= 71)   625     76        0 (1367/1757-1449/1802)  -5.2 W54-W66 N35-N37
      2  641 1632 (N53E10 L=  1)  1461    252        3 ( 551/1566- 754/1684)  -5.4 E20-W01 N50-N55
      3  448 1555 (N53E33 L=339)   585     89        2 ( 409/1521- 476/1589)  38.7 E37-E29 N51-N56
      4 1114 1429 (N25W17 L= 28)   462    136        0 (1078/1365-1138/1476)  66.8 W14-W18 N21-N28
      5  429 1328 (N37E33 L=338)   651    130        3 ( 371/1284- 473/1361) -31.0 E39-E29 N34-N40
      6 1173 1251 (N14W14 L= 26)   448    111        2 (1133/1223-1223/1284) -59.0 W12-W16 N11-N16
      7  239  999 (N21E55 L=316)  1461    332        1 ( 140/ 896- 341/1087) -32.5 E68-E46 N13-N28
      8 1890  988 (S22W61 L= 73)   690     93        1 (1850/ 966-1928/1013) -14.0 W56-W68 S24-S19
      9  191  848 (N13E64 L=307)  3181    401        9 ( 112/ 698- 274/ 976) -67.6 E75-E54 N04-N21
     10 1115  723 (S14E03 L=  9)   992    331        0 (1034/ 593-1206/ 858)  36.3 E10-W05 S20-S08
     11 1340  662 (S24W09 L= 20)   858    196        2 (1256/ 599-1395/ 743)  18.4 W02-W14 S26-S19
     12 1594  646 (S33W27 L= 38)   272     48        1 (1564/ 627-1618/ 663) -21.8 W24-W29 S35-S31
     13  868  607 (S15E20 L=351)   267     92        0 ( 829/ 582- 912/ 640) -63.4 E21-E17 S17-S11
     14  450  489 (S12E53 L=318)   730    161        0 ( 388/ 455- 517/ 543) -49.4 E56-E49 S15-S07
     15 1532  518 (S40W20 L= 32)   298     47        1 (1512/ 497-1554/ 540) -56.3 W18-W22 S41-S38
     16  567  295 (S28E59 L=313)   651     90        1 ( 547/ 247- 598/ 324)  14.9 E67-E52 S30-S25
  

    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