Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  852 1750 (N53E20 L=158)   512     72        1 ( 819/1724- 887/1771)  26.6 E23-E16 N50-N55
      2  347 1535 (N33E62 L=116)  1800    167        4 ( 297/1475- 414/1606) -13.1 E82-E51 N29-N37
      3 1324 1553 (N38W22 L=200)   447     75        2 (1288/1531-1363/1571)  -9.5 W19-W25 N36-N39
      4 1425 1428 (N29W28 L=206)  1064    301        2 (1370/1277-1459/1536)  86.6 W24-W30 N19-N36
      5  612 1461 (N29E31 L=147)  1089    233        2 ( 515/1428- 727/1494)  -6.3 E40-E22 N27-N31
      6  434 1132 (N08E40 L=138)   578    191        1 ( 387/1061- 472/1211) -48.4 E44-E36 N03-N12
      7 1776 1154 (N11W54 L=232)   624    101        0 (1730/1110-1809/1188) -32.0 W49-W57 N08-N13
      8  777 1054 (N04E15 L=163)   188     50        0 ( 757/1027- 804/1074)  33.7 E16-E13 N02-N05
      9 1829  970 (S00W59 L=237)   610     67        0 (1809/ 931-1846/1005)  45.0 W56-W61 S02-N01
     10  390  847 (S10E43 L=135)   795    321        0 ( 336/ 731- 426/ 966) -74.1 E46-E39 S17-S02
     11 1732  759 (S13W51 L=229)   882     96        1 (1692/ 719-1768/ 803) -26.6 W47-W55 S15-S10
     12 1055  567 (S26W03 L=181)  1121    262        2 ( 979/ 526-1165/ 649) -32.5 E02-W10 S28-S20
     13  851  573 (S26E11 L=167)   370     98        0 ( 805/ 548- 894/ 611)  29.7 E14-E07 S27-S23
     14 1660  582 (S25W49 L=227)   365     56        0 (1634/ 557-1691/ 603)  21.8 W47-W52 S26-S23
     15  628  243 (S55E44 L=134)   506     54        0 ( 602/ 226- 659/ 265) -36.9 E47-E40 S57-S53
     16 1363  218 (S56W43 L=221)   569     28        0 (1350/ 194-1382/ 248) -28.6 W38-W49 S58-S52
  

    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