Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  606 1720 (N32E44 L=284)   386     49        0 ( 579/1705- 635/1739)   0.0 E46-E41 N30-N33
      2  307 1469 (N12E58 L=270)   836    111        2 ( 278/1417- 330/1514) -56.3 E62-E54 N08-N14
      3 1166 1399 (N18W03 L=331)   495    191        2 (1075/1382-1267/1420)   5.2 E02-W09 N16-N19
      4  479 1372 (N07E40 L=288)   363     81        1 ( 439/1350- 521/1390)  18.4 E42-E36 N05-N08
      5 1547 1280 (N17W29 L=357)  1771    268        5 (1439/1189-1628/1363) -24.8 W20-W34 N12-N22
      6  349 1299 (N02E48 L=280)   555     89        2 ( 326/1253- 376/1338)  90.0 E50-E46 S01-N04
      7  224 1302 (N01E61 L=267)   512     58        2 ( 193/1276- 264/1322)   0.0 E64-E56 N00-N02
      8 1174 1294 (N12W05 L=333)   200     98        0 (1146/1253-1218/1319)  51.3 W03-W07 N08-N13
      9 1454 1251 (N14W23 L=351)   283     38        0 (1427/1225-1476/1276) -71.6 W20-W24 N12-N15
     10 1445 1200 (N11W22 L=350)   250     66        0 (1411/1180-1475/1219)  26.6 W20-W24 N09-N12
     11 1529 1079 (N05W29 L=358)   242     89        0 (1480/1067-1574/1092)  18.4 W26-W32 N03-N06
     12 1749  890 (S01W49 L= 17)   439     89        0 (1705/ 873-1798/ 910)  14.0 W45-W53 S02-N00
     13  459  817 (S26E35 L=293)   928    239        0 ( 366/ 758- 574/ 855)  21.3 E43-E25 S30-S22
     14 1267  642 (S25W22 L=350)  3671    898        4 ( 981/ 563-1668/ 804)   8.1 W02-W51 S29-S17
     15  499  655 (S36E32 L=296)   347     76        1 ( 464/ 638- 531/ 680)  31.0 E35-E29 S36-S33
     16 1825  595 (S15W71 L= 39)   634     87        1 (1790/ 578-1869/ 619)  24.0 W66-W77 S16-S12
     17  446  536 (S45E40 L=288)  1312    331        5 ( 290/ 507- 584/ 561)   1.7 E60-E25 S46-S42
     18  719  507 (S43E14 L=314)   466    112        1 ( 662/ 478- 779/ 533)  -5.7 E18-E08 S43-S41
     19  955  451 (S42W06 L=334)   218     56        1 ( 925/ 439- 989/ 469) -11.3 W03-W08 S43-S41
  

    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