Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1667 1603 (N22W69 L= 35)  1844    216        3 (1629/1510-1698/1694)  24.2 W59-W79 N17-N27
      2 1502 1602 (N27W52 L= 19)  1306    164        1 (1432/1545-1566/1654) -73.3 W49-W54 N22-N32
      3 1342 1437 (N23W31 L=358)   337     84        1 (1315/1395-1369/1467)  31.0 W29-W34 N20-N24
      4 1199 1430 (N26W22 L=348)   336     84        0 (1165/1389-1226/1465) -80.5 W20-W22 N23-N29
      5 1503 1411 (N17W41 L=  8)   803    237        0 (1408/1341-1581/1456) -57.5 W36-W44 N10-N21
      6  664 1372 (N36E16 L=311)   315     52        1 ( 636/1354- 692/1391)  11.3 E18-E13 N34-N37
      7 1137 1315 (N21W14 L=341)   226     52        1 (1113/1291-1156/1341)  71.6 W12-W15 N19-N22
      8 1601 1080 (S04W37 L=  3)   330    126        0 (1563/1035-1632/1128)  23.2 W32-W39 S05-S01
      9  192 1050 (N24E63 L=263)  1325    128        2 ( 139/1007- 251/1090) -13.2 E72-E55 N20-N26
     10  128  997 (N21E75 L=251)  1039    104        0 ( 106/ 947- 147/1037)  31.0 E82-E70 N17-N23
     11  906  841 (S01E11 L=316)   291     60        1 ( 884/ 809- 933/ 867) -90.0 E12-E09 S02-N01
     12 1515  785 (S19W25 L=351)   978    331        2 (1366/ 742-1590/ 864)  -7.1 W14-W30 S21-S15
     13 1056  808 (S06E03 L=324)   225    138        0 (1001/ 765-1114/ 845)  14.0 E07-W01 S07-S04
     14 1707  772 (S26W40 L=  7)  1336    228        2 (1622/ 718-1808/ 835)   0.0 W32-W51 S27-S23
     15  828  713 (S07E19 L=308)   365     98        2 ( 800/ 675- 857/ 761)  51.3 E21-E16 S08-S03
     16  866  289 (S36E33 L=293)  2912    460        4 ( 799/ 125- 960/ 506)  30.1 E66-E16 S51-S21
  

    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