Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1099 1856 (N56W25 L=113)   466     55        1 (1061/1842-1134/1867) -23.2 W21-W29 N54-N57
      2 1469 1685 (N36W46 L=133)  1554    251        2 (1365/1645-1577/1739)   4.1 W35-W63 N33-N38
      3  543 1496 (N35E29 L= 58)   410     69        0 ( 499/1468- 582/1518) -26.6 E33-E25 N32-N36
      4  754 1410 (N26E13 L= 74)   325    144        0 ( 680/1380- 807/1430) -26.6 E17-E09 N23-N28
      5 1659 1397 (N14W48 L=135)   398    103        0 (1610/1363-1694/1426) -33.7 W44-W50 N12-N16
      6  384 1349 (N27E41 L= 46)   718    204        0 ( 300/1296- 451/1382) -25.0 E50-E35 N22-N30
      7  391 1084 (N10E40 L= 47)   542    141        1 ( 361/1029- 415/1146)  71.6 E42-E38 N07-N13
      8  256 1095 (N13E52 L= 35)   629     99        3 ( 224/1055- 292/1132)  53.1 E55-E48 N10-N14
      9  288 1042 (N09E49 L= 38)   848    163        3 ( 242/ 980- 337/1102)  54.5 E53-E44 N05-N12
     10 1848  962 (S13W59 L=146)   611     95        1 (1834/ 908-1865/1006)  90.0 W56-W60 S16-S10
     11  349  916 (N01E45 L= 42)   846    220        3 ( 284/ 865- 419/ 961)   0.0 E50-E39 S02-N04
     12 1390  872 (S14W21 L=108)   369    138        3 (1343/ 829-1436/ 929)  39.8 W17-W24 S15-S10
     13  736  647 (S20E23 L= 64)   312    136        0 ( 700/ 599- 788/ 676)  33.7 E26-E19 S22-S18
     14 1455  610 (S31W25 L=112)   773    183        0 (1370/ 564-1518/ 666)  20.0 W18-W29 S33-S28
     15 1339  621 (S29W16 L=103)   209     43        0 (1316/ 593-1362/ 647)  45.0 W14-W17 S30-S27
     16  756  556 (S26E24 L= 63)   347    114        2 ( 710/ 508- 788/ 594) -74.1 E26-E22 S29-S22
  

    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