Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  974 1467 (N32E02 L= 29)   322     56        0 ( 951/1430-1002/1498) -59.0 E03-W00 N29-N34
      2  262 1421 (N30E65 L=326)  1593    133        3 ( 229/1361- 299/1481) -19.3 E76-E56 N26-N33
      3  868 1416 (N29E09 L= 21)   334    140        1 ( 828/1358- 898/1444) -45.0 E12-E07 N25-N31
      4  764 1302 (N22E16 L= 14)   266     76        1 ( 740/1278- 790/1324) -33.7 E17-E14 N20-N23
      5  500 1269 (N20E35 L=355)  1020    233        0 ( 403/1241- 602/1298)   0.0 E43-E27 N18-N22
      6 1095 1249 (N18W05 L= 36)   205     85        0 (1055/1238-1132/1265) -11.3 W02-W07 N17-N18
      7 1611  966 (S02W38 L= 69)   841    236        0 (1590/ 857-1633/1074) -90.0 W36-W39 S08-N04
      8  614  972 (N02E26 L=  5)   225    112        0 ( 561/ 953- 668/ 994) -15.9 E29-E22 N00-N03
      9 1243  952 (S01W13 L= 44)   192     40        0 (1227/ 924-1262/ 976)  63.4 W12-W14 S02-N00
     10  789  808 (S08E15 L= 15)   438    136        0 ( 735/ 764- 846/ 841)  26.6 E19-E11 S10-S06
     11 1528  802 (S12W32 L= 63)   698    174        1 (1449/ 754-1606/ 836)  20.0 W26-W38 S14-S10
  

    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