Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1416 1808 (N39W63 L=288)  2699    250        5 (1319/1775-1527/1840) -13.0 W50-W76 N34-N42
      2 1660 1630 (N21W67 L=292)  2765    329        1 (1567/1528-1791/1729) -70.6 W63-W76 N11-N28
      3 1300 1624 (N31W36 L=262)   640    164        2 (1220/1606-1389/1641) -31.0 W31-W41 N27-N34
      4  776 1442 (N34E05 L=220)  1461    665        0 ( 493/1378-1057/1487) -15.8 E27-W12 N26-N38
      5 1294 1292 (N12W22 L=247)   546    170        0 (1255/1219-1338/1358) -90.0 W21-W23 N07-N16
      6  862 1257 (N21E04 L=221)   193    111        0 ( 814/1226- 907/1290)   8.1 E07-W00 N19-N21
      7 1787 1119 (S10W52 L=277)   439     96        0 (1761/1074-1810/1161) -78.7 W50-W54 S12-S07
      8  214  899 (N14E59 L=166)  1080    125        2 ( 153/ 857- 266/ 943) -28.6 E66-E53 N10-N17
      9  935  825 (S05E10 L=216)   467    136        1 ( 899/ 755- 963/ 871) -90.0 E11-E08 S09-S02
     10 1842  777 (S33W55 L=280)   799    110        2 (1822/ 738-1876/ 837) -42.0 W52-W62 S36-S27
     11 1198  640 (S23W01 L=226)   841    204        0 (1143/ 565-1235/ 734)  41.6 E04-W05 S26-S17
     12 1685  694 (S33W37 L=262)   387     64        0 (1652/ 677-1718/ 712) -18.4 W34-W40 S34-S31
     13 1598  530 (S42W29 L=254)  1353    197        1 (1523/ 470-1655/ 610)  24.8 W21-W34 S44-S37
     14  407  427 (S18E62 L=163)   473     46        1 ( 391/ 397- 421/ 456)  33.7 E65-E58 S19-S15
     15  551  380 (S23E52 L=174)   499    116        0 ( 538/ 327- 572/ 433)  37.9 E55-E46 S27-S20
     16  678  329 (S30E45 L=180)   437     76        0 ( 657/ 290- 698/ 368)  32.0 E49-E40 S32-S27
     17 1158  254 (S50E16 L=209)  1710    288        0 (1031/ 202-1289/ 310)   0.0 E29-E02 S51-S47
  

    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