Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1854 (N40E47 L=206)   567     92        1 ( 841/1804- 865/1898) -36.9 E54-E41 N35-N44
      2  712 1816 (N33E53 L=199)   949    126        3 ( 683/1762- 746/1866) -48.4 E60-E48 N28-N37
      3 1608 1663 (N56W26 L=279)   779     48        0 (1587/1629-1626/1701)  51.3 W22-W30 N51-N61
      4  657 1582 (N15E39 L=213)   434     70        0 ( 620/1550- 686/1611)  -9.5 E42-E36 N13-N17
      5  508 1255 (S07E37 L=215)   618    138        1 ( 444/1222- 567/1285)  35.5 E40-E32 S10-S04
      6 1603 1210 (N21W29 L=282)   806    213        0 (1525/1125-1651/1267) -20.0 W22-W34 N17-N23
      7  277 1083 (S22E52 L=200)  1163    274        1 ( 202/ 972- 363/1174)  49.1 E59-E45 S29-S14
      8  184 1110 (S22E64 L=189)   517    104        0 ( 153/1070- 210/1159)  45.0 E68-E61 S24-S18
      9  415  976 (S25E37 L=215)   277     50        2 ( 392/ 956- 446/ 998)  53.1 E39-E35 S26-S22
     10  609  810 (S30E19 L=233)  1278    356        0 ( 463/ 741- 770/ 860)   7.1 E31-E07 S31-S27
     11 1330  464 (S28W39 L=291)   394    109        1 (1286/ 426-1379/ 501)  71.6 W36-W40 S30-S24
     12  746  242 (S67W17 L=269)   459     75        0 ( 732/ 206- 759/ 286) -31.6 W10-W23 S71-S62
  

    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