Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  997 1870 (N60W31 L=224)  1730    244        1 ( 920/1802-1098/1933)  10.4 W17-W55 N55-N62
      2 1211 1737 (N42W35 L=228)   302     43        0 (1187/1712-1238/1759)   8.1 W31-W38 N40-N42
      3  238 1142 (N24E54 L=139)   527    123        0 ( 226/1087- 254/1213) -76.0 E55-E52 N20-N28
      4 1651 1017 (S11W39 L=232)   250     95        0 (1622/ 984-1679/1063) -63.4 W37-W40 S13-S07
      5 1873  720 (S35W64 L=258)  1330     86        1 (1841/ 692-1911/ 769) -45.0 W58-W76 S37-S31
      6  808  710 (S12E19 L=174)   232     97        0 ( 787/ 662- 831/ 758) -90.0 E19-E18 S14-S08
      7 1175  708 (S19W02 L=196)   244     62        1 (1142/ 695-1213/ 728) -14.0 W00-W04 S20-S17
      8 1045  615 (S22E07 L=186)   244     87        0 (1008/ 590-1082/ 639)   9.5 E10-E04 S23-S20
      9  874  347 (S36E28 L=165)  1261    312        1 ( 730/ 328-1028/ 377) -25.3 E37-E16 S40-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