Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1621 1663 (N50W35 L=195)   970    131        2 (1573/1636-1658/1687)  15.3 W29-W40 N48-N54
      2  886 1400 (N12E15 L=144)   799    255        1 ( 863/1285- 912/1522) -61.7 E19-E11 N05-N19
      3  240 1363 (N01E63 L= 97)   700    115        1 ( 211/1332- 265/1395) -20.6 E67-E59 N00-N03
      4 1569 1303 (N21W29 L=188)   369    114        1 (1553/1257-1584/1361)  90.0 W27-W29 N17-N24
      5 1096 1332 (N12E02 L=158)   256     89        0 (1059/1293-1139/1359)  59.0 E03-W00 N09-N14
      6 1153 1242 (N08W03 L=163)   478    166        0 (1113/1172-1184/1331) -84.3 W01-W04 N03-N13
      7  506 1308 (N01E37 L=122)   231     27        0 ( 484/1292- 525/1323)  26.6 E38-E35 N00-N02
      8  359 1195 (S08E46 L=113)   309     16        0 ( 336/1176- 378/1216) -18.4 E47-E44 S08-S06
      9 1593  726 (S11W42 L=201)   906    214        3 (1505/ 626-1657/ 816)  74.1 W39-W44 S19-S05
     10 1034  735 (S23W06 L=166)   617    177        1 ( 931/ 710-1113/ 770)  31.0 W00-W10 S26-S19
     11 1448  607 (S22W36 L=196)   265     32        1 (1422/ 590-1474/ 627)  18.4 W34-W38 S22-S20
     12 1249  600 (S27W24 L=183)   267     47        0 (1221/ 585-1282/ 618)  33.7 W21-W25 S28-S25
     13 1461  325 (S37W52 L=211)  2322    450        2 (1283/ 300-1636/ 362)  27.5 W39-W65 S43-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