Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1489 1583 (N37W24 L=133)   305     27        0 (1469/1563-1509/1601)  56.3 W22-W26 N35-N38
      2  233 1348 (N05E62 L= 47)   446     93        0 ( 202/1334- 263/1366)   8.1 E65-E58 N03-N06
      3 1567 1345 (N22W30 L=139)   292     88        0 (1529/1333-1612/1358)  18.4 W27-W33 N20-N22
      4 1110 1212 (N06W02 L=111)   195     99        0 (1059/1193-1166/1229)  31.0 W00-W05 N04-N07
      5 1472 1074 (N04W26 L=135)   581    151        0 (1434/ 998-1507/1144)  77.5 W24-W27 S01-N08
      6  115 1068 (S12E70 L= 38)   854    174        1 (  92/1029- 137/1110)  51.3 E73-E67 S14-S09
      7 1649 1084 (N07W38 L=147)   302     70        1 (1622/1063-1679/1105)  26.6 W36-W40 N06-N09
      8  239 1052 (S13E55 L= 54)   542    130        0 ( 211/ 991- 275/1096)  50.2 E57-E51 S16-S10
      9  315 1014 (S15E47 L= 61)   508     80        1 ( 283/ 981- 348/1051) -45.0 E50-E44 S16-S12
     10  873  969 (S11E08 L=100)   731    293        2 ( 759/ 885- 994/1025)  34.7 E14-E01 S17-S06
     11 1006  922 (S12W00 L=109)   632    173        4 ( 953/ 871-1049/ 958) -21.8 E03-W03 S15-S10
     12  352  871 (S23E44 L= 65)   656    213        0 ( 285/ 815- 414/ 928)  32.5 E49-E38 S26-S19
     13 1092  823 (S17W07 L=116)   187     66        0 (1068/ 792-1115/ 854)  78.7 W05-W08 S19-S14
     14  970  564 (S35W04 L=113)   292     61        0 ( 935/ 540-1007/ 587)  38.7 W01-W06 S36-S32
     15 1451  507 (S30W39 L=148)  1429    251        2 (1340/ 450-1569/ 579)  45.0 W32-W45 S35-S23
     16 1111  357 (S47W21 L=130)   282     25        0 (1094/ 336-1128/ 382) -76.0 W19-W22 S48-S44
  

    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