Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1069 1628 (N38W23 L=160)   948    215        3 (1019/1548-1121/1699)  23.2 W16-W30 N33-N41
      2 1479 1604 (N24W48 L=185)  1096    189        3 (1412/1561-1558/1634) -56.3 W42-W52 N18-N27
      3  545 1251 (N30E24 L=113)  1020    394        1 ( 390/1222- 730/1279) -19.1 E36-E10 N26-N35
      4 1539 1166 (S01W33 L=170)  1070    259        3 (1429/1120-1629/1201) -12.1 W25-W39 S05-N01
      5 1497 1092 (S04W29 L=166)   247    104        0 (1457/1058-1538/1110)   0.0 W25-W31 S05-S02
      6 1702  988 (S16W40 L=177)   584     99        2 (1656/ 962-1737/1018)  -8.1 W36-W43 S17-S14
      7 1431  837 (S16W18 L=155)   495    171        1 (1414/ 767-1455/ 933)  73.3 W16-W20 S20-S10
      8 1098  694 (S15E05 L=132)   788    209        1 ( 974/ 655-1188/ 721) -39.3 E11-W00 S19-S10
      9 1199  557 (S26E03 L=134)   264     75        0 (1185/ 510-1216/ 589)  53.1 E05-E01 S28-S23
     10  669  438 (S19E39 L= 98)   404    117        2 ( 635/ 409- 697/ 477)  78.7 E41-E37 S21-S16
     11 1041  218 (S46E29 L=108)  1086    165        1 (1003/ 172-1068/ 268)  26.6 E35-E23 S49-S41
  

    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