Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  777 1650 (N27E29 L=195)   738    174        1 ( 735/1574- 809/1727) -36.9 E36-E24 N22-N32
      2  521 1583 (N19E45 L=180)   281     92        0 ( 479/1564- 548/1600)  21.8 E47-E42 N17-N19
      3  468 1494 (N12E45 L=179)   356     86        0 ( 432/1474- 502/1516)   9.5 E48-E42 N11-N14
      4  248 1339 (N01E60 L=165)  1153    247        2 ( 231/1237- 279/1464) -69.0 E63-E55 S05-N08
      5 1411 1309 (N17W18 L=242)   692    195        1 (1311/1287-1495/1328)  22.6 W11-W23 N14-N19
      6  274 1133 (S11E52 L=172)   304     38        1 ( 254/1114- 296/1162) -45.0 E53-E49 S12-S09
      7  378 1039 (S16E42 L=182)   904    283        0 ( 280/ 952- 452/1148)  47.3 E50-E37 S21-S08
      8  221 1097 (S14E57 L=167)   620    107        3 ( 191/1054- 249/1145) -51.3 E60-E53 S16-S10
      9  479 1123 (S09E35 L=189)   213     10        0 ( 462/1104- 497/1143) -63.4 E36-E33 S10-S08
     10  899 1028 (S09E07 L=217)   252     56        1 ( 872/1006- 923/1055)  45.0 E08-E05 S09-S07
     11 1200  860 (S13W13 L=238)   351    122        0 (1152/ 840-1254/ 894)   9.5 W10-W16 S14-S10
     12  319  727 (S35E48 L=176)   953    133        1 ( 275/ 678- 383/ 774)  23.2 E54-E40 S38-S32
  

    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