Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1916 (N60W52 L=258)  1246    124        3 (1058/1873-1156/1958)   7.7 W38-W75 N57-N62
      2  995 1881 (N61W33 L=240)   490     65        0 ( 973/1848-1015/1906)  35.0 W29-W39 N56-N63
      3  875 1816 (N59W14 L=221)   867    129        1 ( 829/1769- 929/1858)  17.4 W07-W23 N55-N61
      4 1661 1068 (S08W41 L=247)   544    105        2 (1629/1017-1698/1109)  23.2 W37-W44 S10-S06
      5 1781  687 (S35W50 L=257)   891    134        0 (1739/ 654-1835/ 736) -23.2 W44-W58 S37-S30
      6  976  638 (S19E11 L=195)   262     41        2 ( 948/ 619-1003/ 659) -33.7 E12-E09 S20-S17
      7  876  557 (S22E19 L=187)   346    117        0 ( 828/ 527- 932/ 579)   7.1 E23-E15 S23-S20
      8  693  304 (S35E45 L=162)  1118    212        0 ( 603/ 282- 798/ 341) -28.3 E50-E36 S38-S30
      9 1406  281 (S58W13 L=219)   542     67        2 (1364/ 266-1446/ 302) -10.3 W07-W18 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