Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  395 1090 (N24E40 L=317)   630    159        0 ( 336/1052- 457/1119) -29.1 E44-E35 N21-N26
      2  322  977 (N18E48 L=309)   634    133        0 ( 278/ 938- 376/1021) -40.6 E51-E43 N15-N21
      3 1404  925 (S08W20 L= 17)   400     96        0 (1357/ 893-1452/ 951)   0.0 W16-W22 S08-S06
      4 1507  910 (S11W26 L= 23)   212     51        1 (1488/ 884-1525/ 939)  71.6 W24-W27 S12-S09
      5  685  806 (N02E25 L=332)   849    190        1 ( 660/ 718- 713/ 893)  71.6 E28-E22 S02-N07
      6 1229  780 (S12W06 L=  3)   683    173        1 (1192/ 691-1257/ 869)  59.0 W01-W08 S16-S06
      7 1562  749 (S23W27 L= 24)   673    235        0 (1509/ 632-1587/ 867)  55.3 W21-W30 S28-S15
      8 1622  627 (S33W30 L= 27)   332     56        0 (1588/ 596-1647/ 654)  11.3 W27-W32 S34-S31
      9 1200  517 (S27E02 L=355)   918    287        0 (1061/ 415-1321/ 616)  15.9 E13-W08 S31-S23
     10  865  441 (S24E25 L=332)   343     75        0 ( 835/ 410- 901/ 473)  26.6 E27-E21 S25-S22
     11  850  238 (S39E38 L=320)   813    254        1 ( 827/ 172- 874/ 307)  40.2 E44-E31 S45-S33
  

    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