Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1576 1692 (N49W36 L=  7)   420     84        0 (1561/1661-1596/1715)  68.2 W33-W39 N46-N51
      2  569 1592 (N23E40 L=291)   765    189        3 ( 515/1539- 616/1625) -29.1 E44-E35 N19-N25
      3  372 1535 (N18E54 L=277)  2274    307        6 ( 303/1419- 473/1612)  58.4 E57-E47 N09-N23
      4  711 1519 (N20E27 L=304)   482     75        1 ( 691/1486- 735/1558) -36.9 E29-E25 N17-N22
      5  579 1404 (N11E33 L=298)  1079    267        3 ( 507/1327- 697/1460)  42.0 E37-E26 N05-N14
      6  494 1402 (N10E39 L=292)   266    103        1 ( 475/1353- 524/1454)  90.0 E40-E38 N06-N13
      7 1179 1395 (N18W04 L=335)   347     98        0 (1151/1343-1204/1444) -56.3 W01-W06 N15-N21
      8  236 1276 (N01E59 L=272)   390     71        0 ( 223/1239- 250/1317) -63.4 E60-E57 S01-N03
      9 1108 1270 (N09W01 L=333)   325     74        1 (1059/1246-1139/1299) -26.6 E01-W03 N08-N10
     10  153 1086 (S11E65 L=266)  1745    398        1 (  72/1053- 225/1127)  10.8 E78-E56 S12-S08
     11 1339  979 (S04W19 L=350)   270     89        0 (1296/ 945-1385/1004) -18.4 W16-W22 S05-S03
     12  549  857 (S22E28 L=303)   773    198        1 ( 467/ 820- 647/ 901)  21.0 E34-E21 S24-S18
     13 1310  800 (S15W21 L=352)   240     56        0 (1280/ 776-1348/ 825)  53.1 W19-W22 S16-S12
     14 1789  690 (S12W60 L= 31)  1394    189        3 (1741/ 618-1831/ 785)  79.7 W56-W62 S16-S05
     15  322  678 (S35E48 L=284)   302     50        1 ( 305/ 648- 339/ 705) -90.0 E49-E46 S37-S33
     16  992  597 (S32W05 L=336)   413     79        2 ( 956/ 565-1024/ 631) -18.4 W01-W07 S34-S30
  

    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