Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  843 1823 (N62E25 L= 13)  1243    167        1 ( 781/1801- 921/1849)   3.0 E35-E13 N59-N65
      2 1163 1729 (N53W14 L= 52)   607     84        1 (1124/1706-1209/1748)  21.8 W09-W19 N50-N54
      3  548 1461 (N31E37 L=  2)   326     75        0 ( 524/1429- 575/1491) -45.0 E39-E33 N28-N32
      4 1181 1385 (N26W11 L= 49)   693    164        0 (1108/1329-1227/1422) -32.0 W05-W13 N22-N28
      5 1758  854 (S08W53 L= 91)   859    123        1 (1721/ 829-1820/ 880)  11.3 W49-W59 S09-S06
      6 1382  757 (S13W23 L= 62)   590    199        0 (1350/ 695-1431/ 835) -71.6 W21-W27 S17-S08
      7 1676  709 (S17W47 L= 86)   302     90        0 (1656/ 672-1691/ 758)  78.7 W46-W48 S19-S14
      8  817  688 (S18E13 L= 25)   778    174        1 ( 744/ 660- 901/ 718)  15.3 E18-E07 S19-S15
      9  380  518 (S31E54 L=345)  2494    319        6 ( 300/ 438- 473/ 623) -50.2 E60-E45 S37-S23
     10 1386  567 (S26W26 L= 64)  1069    280        1 (1275/ 534-1518/ 600)   3.0 W17-W36 S28-S23
  

    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