Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1076 1741 (N46W31 L=139)   371     93        1 (1054/1702-1108/1774)  18.4 W26-W35 N43-N47
      2 1246 1676 (N36W37 L=146)   778    206        1 (1154/1633-1329/1717) -61.4 W32-W41 N30-N41
      3  837 1636 (N46W07 L=115)   282     59        0 ( 806/1620- 876/1652)  -9.5 W03-W09 N45-N47
      4 1500 1584 (N22W49 L=157)   308     77        0 (1473/1550-1523/1622)  23.2 W45-W52 N21-N24
      5 1357 1579 (N27W38 L=147)   311     50        2 (1334/1550-1378/1610)  36.9 W36-W40 N24-N28
      6  446 1303 (N35E32 L= 76)   424    141        0 ( 408/1243- 496/1352) -52.1 E35-E28 N29-N38
      7  808 1285 (N25E06 L=102)   315     97        1 ( 785/1244- 834/1334) -90.0 E07-E05 N21-N27
      8  352 1172 (N29E42 L= 67)   269    113        0 ( 309/1150- 399/1189) -14.0 E45-E37 N27-N29
      9  251 1079 (N25E53 L= 56)   932    150        2 ( 209/1046- 317/1121) -22.6 E57-E45 N23-N28
     10  859 1068 (N11E08 L=100)   197     52        0 ( 845/1039- 871/1107)  76.0 E09-E07 N09-N13
     11  143 1020 (N24E68 L= 41)  1462    163        5 ( 112/ 974- 190/1101) -31.6 E74-E60 N20-N28
     12  974  785 (S08E09 L=100)   627    190        1 ( 934/ 700-1029/ 861)  41.2 E13-E03 S11-S03
     13  224  713 (N04E65 L= 44)   458     65        1 ( 203/ 679- 248/ 744) -53.1 E66-E62 N02-N06
     14 1000  468 (S27E17 L= 91)   450    121        2 ( 958/ 450-1066/ 502) -26.6 E19-E13 S29-S24
     15  753  405 (S25E35 L= 73)   372    119        0 ( 696/ 387- 813/ 418)  -6.3 E40-E30 S25-S23
     16 1147  199 (S54E23 L= 85)   520     60        0 (1114/ 181-1179/ 215)  -9.5 E26-E19 S55-S51
     17 1056  131 (S57E40 L= 68)   600     41        1 (1030/ 108-1076/ 156)  56.3 E46-E35 S60-S54
  

    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